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FOREWORD 


This  final  technical  report  was  prepared  under  contract 
F33615  67C  1870,  University  of  Miami,  Florida.  The 
symposium  covered  herein  was  accomplished  on  Project 
7114  under  the  technical  cognizance  of  Dr.  Kenneth  E.  Kissell, 
Director  of  the  General  Physics  Research  Laboratory, 
Aerospace  Research  Laboratories. 
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TEN  YEARS  OF  OPTICAL  SATELLITE  OBSERVATIONS 

Douglas  Duke 
University  of  Miami 

It  is  not  generally  known  that  within  a  matter  of  weeks  after  the  launching 
of  the  first  two  earth  satellite  vehicles  by  Russia  late  in  1957,  large  scale 
photographs  were  made  of  these  objects  at  the  Air  Force  Missile  Test  Center 
at  Cape  Kennedy,  Florida.  Even  the  very  earliest  photographs  allowed  a 
reasonably  accurate  measurement  of  the  size  of  these  objects  and  their 
orientation  changes  (tumbling  rate) .  Much  experimentation  was  carried 
out  during  the  first  six  months  in  an  effort  to  obtain  the  clearest  possible 
photographs.  Various  combinations  of  exposure  times,  focal  lengths  and 
film  speeds  were  tried.  Methods  for  rapid  acquisition  and  the  smoothest 
possible  tracking  were  developed. 

Despite  all  these  efforts ,  even  the  clearest  photographs  made  in  1958 
did  not  have  sufficient  resolution  to  allow  the  determination  of  structural 
details  below  approximately  5  feet  in  extent.  While  some  improvement 
might  be  expected  from  increasing  the  telescope  aperture  to  more  than 
the  24  inch  instruments  then  available,  such  improvement  would  be  nominal 
since  the  best  observed  resolution  was  several  times  the  diffraction  limit 
of  the  24  inch  instruments.  As  is  the  case  in  ail  astronomical  observations 
where  high  resolution  is  desirable,  the  problem  lay  in  the  distortion  of  the 
optical  wavefronts  produced  by  the  refractive  index  turbulence  in  the  earth's 
atmosphere . 

Two  independent  and  complementary  actions  were  called  for  to  improve  this 
problem  of  atmospheric  seeing.  The  first  was  to  design  tracking  telescopes 
with  apertures  and  driving  mechaAiams  optimized  for  satellite  observations  and 
to  locate  these  at  high  elevation  sites  so  as  to  minimize  the  atmospheric 
transmission  path  and  the  resulting  distortion.  This  consideration  resulted  in 
the  construction  of  the  Cloudcroft  Electro-Optical  Facility  by  the  Air  Force 
and  the  Mount  Haleakala  Facility  by  the  Advanced  Research  Projects  Agency. 


v 


4 


♦ 

i 

BLANK  PAGE 


The  second  indicated  action  was  a  comprehensive  program  of  study  and 
experimentation’ for  the  purpose  of  improving  our  knowledge  of  the  exact  nature 
of  the  atmospheric  distortions  so  as  to  attempt  to  devise  new  devices  and 
techniques  reducing  their  effect.  Unfortunately  the  pleas  and  urgings  of  the 
author  fell  on  deaf  ears  for  almost  five  years.  Finally  in  1963  the  increased 
realization  in  the  possibility  of  optical  methods  due  to  the  new  attention 
being  paid  to  the  laser  caused  a  revival  of  interest. 

For  the  past  18  months  a  group  under  the  direction  of  the  author  and  with  the 
support  of  the  Air  Force  Avionics  Laboratory  and  the  Air  Force  Eastern 
Test  Range  has  been  working  on  these  problems  of  atmospheric  distortion 
to  optical  images.  To  date  the  results  have  been  modest,  principally  due 
to  the  difficult  nature  of  the  phenomena  being  investigated  and  the  need  to 
design  and  construct  entirely  new  registration  and  measurement  devices. 

It  is  felt  that  only  a  beginning  has  been  made  and  that  real  results  will  take 
several  additional  years. 
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OPTICAL  AND  INFRARED  CHARACTERISTICS  OF 
SATELLITE  SURFACE  MATERIALS 

S.  A.  Greenberg,  Research  Scientist 

LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 
PALO  ALTO,  CALIFORNIA 


The  optical  properties  of  typical  spacecraft  surface  materials  are  discussed. 
The  four  basic  classes  of  materials  (solar  absorbers,  solar  reflectors,  flat 
absorbers  and  flat  reflectors)  are  described  in  terms  of  directional  spectral 
reflectance  in  the  wavelength  region  0.2  -  30  microns.  The  concept  of  remote 
spectral  scanning  of  satellites  for  identification  of  surface  materials  is 
examined  in  relation  to  identification  of  vehicle  origin,  temperature , 
mission  and  lifetime.  Environmentally  induced  degradation  of  optical  proper¬ 
ties  of  spacecraft  coatings  is  discussed  in  reference  to  signature  characteris¬ 
tics.  Methods  of  remote  sensing  of  surface  temperatures  are  proposed,  based 
on  the  temperature  variation  of  optical  absorption  edges  of  semiconductor 
materials.  The  limitations  of  remote  optical  characterization  of  orbiting 
spacecraft  are  outlined  and  requirements  for  closer  coordination  of  designers 
and  observers  are  emphasized. 
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INTRODUCTION 


Late  in  1967,  the  National  Aeronautics  and  Space  Administration  will 
launch  the  third  satellite  in  the  ATS  series,  the  spin-stabilized  syn¬ 
chronous  orbit  ATS-C.  Among  the  several  experiments  carried  by  this 
spacecraft  will  be  the  Ref lectometer  Experiment,  an  experiment  designed 
primarily  for  obtaining  engineering  data  necessary  for  the  design  of 
spacecraft  solar  power  systems  and  other  reflective  optical  systems. 

Power  systems  for  advanced  space  missions  are  recognized  to  be  one  of 
the  most  critical  problem  areas  in  space  technology.  The  two  basic 
energy  sources  apparently  showing  the  most  promise  of  meeting  space 
power  requirements  are  solar  energy  and  nuclear  energy.  One  critical 
component  of  the  solar  energy  system  is  the  solar  concentrator,  which 
depends  upon  geometrical  accuracy  and  spectral  and  specular  reflectance 
of  the  concentrator  surface.  The  efficiency  is  directly  related  to  the 
specular  reflectivity  of  its  mirrored  surface. 

Considerable  effort  has  already  been  expended  in  developing  fabrication 
techniques  for  obtaining  the  essential  geometric  accuracies,  and  in 
investigating  materials  which  demonstrate  in  a  laboratory  environment 
the  high  specular  reflectivity  required  for  space  solar  systems.  How¬ 
ever,  almost  without  exception,  data  regarding  the  performance  and  the 
useful  life  of  highly  reflective  surfaces  in  a  space  environment  is  not 
available.  As  technological  advances  in  power  systems  are  made,  knowl¬ 
edge  about  the  durability  of  materials  used  in  all  types  of  optical 
systems  is  not  only  desirable  but  indispen.'ible . 

Essential  to  the  design  of  solar  space  power  systems  and  other  systems 
requiring  the  use  of  specularly  reflective  surfaces  in  space  is  the 
determination  of  how  reflective  surfaces  are  affected  when  exposed  for 
prolonged  periods  to  a  space  environment. 
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The  difficulty  in  estimating  specular  reflectance  loss  in  a  space  envi¬ 
ronment  results  from  the  following  conditions.  Firsts  the  melanism  and 
environmental  influences  which  cause  reflectance  are  understood  only 
hypothetically.  Second,  flux-intensity  levels  of  possible  degrading 
factors  in  space  generally  have  not  been  determined  directly  but  have 
only  been  estimated  from  extrapolation  of  other  measurements.  Third, 
because  of  the  uncertainties  regarding  the  degrading  mechanisms  them¬ 
selves  and  the  character  of  the  space  environment,  it  simply  has  not 
been  practical  to  simulate  reflective-surface  degradation  in  the 
laboratory. 

The  hostile  aspects  of  the  space  environment  which  have  significant  ef¬ 
fects  on  reflective  surfaces  have  been  determined  to  be  micrometeoroids , 
charged  particles,  ultraviolet  radiation,  temperature,  and  the  inherent 
vacuum  of  space.  In  space,  measurement  of  the  total  effect  of  these 
environmental  factors  is  essential  to  the  development  of  various  systems 
and  instruments  which  depend  upon  specularly  reflecting  surfaces  to 
insure  adequate  performance  throughout  the  planned  lifetime  of  the  system 
or  instrument. 

One  example  of  a  degrading  mechanism  is  micrometeorites.  Most  space 
measurements  made  to  date  have  been  concerned  with  the  penetration  of 
pressure  vessels  (such  as  space  power-system  radiators  and  pressured 
space  vehicles)  by  large  particles.  As  a  result,  most  mic rometcor ite 
instrumentation  has  been  designed  to  detect  encounters  with  these  larger 
particles  in  space.  With  reflective  surfaces,  however,  the  problem  is 
not  penetration  but  gradual  erosion  and  roughening  of  the  surface  layer 
by  the  continual  impact  of  smaller  particles  which  cause  no  major  struc¬ 
tural  damage  such  as  penetration.  This  effect  can  be  especially  pro¬ 
nounced  in  samples  having  interference-type  surface  coatings  which  are 
especially  sensitive  to  changes  in  thickness.  In  extreme  cases,  the 
entire  reflective  coating  may  be  worn  away,  exposing  the  reflective 
substrate  material. 
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Although  many  attempts  have  been  made  to  simulate  various  space-environment 
factors,  there  is  not  only  uncertainty  as  to  what  factor  to  simulate,  but 
to  what  degree.  The  micrometeorite  conditions  suspected  to  exist  in 
space  involve  large  flux  levels  of  micron-sized  particles  moving  at  20 
to  70  kilometers  per  second  relative  to  the  surface  of  the  spacecraft. 

It  has  been  impossible  to  simulate  small  particles  at  high  speeds  simul¬ 
taneously  in  the  laboratory.  Slow  moving  large  particles  with  the  equiva¬ 
lent  kinetic  energy  do  not  have  the  same  effect  and  the  validity  of  the 
extrapolation  of  these  parameters  to  simulate  space  effects  is  the  subject 
of  much  controversy. 

Much  of  the  same  sort  of  argument  applies  to  other  degrading  influences 
assumed  to  exist  in  space.  An  area  of  uncertainty  even  greater  "han  that 
associated  with  the  individual  effects  of  erosion  mechanisms  is  their  com¬ 
bined  effect.  All  these  uncertainties  taken  together  weigh  heavily  in 
favor  of  the  direct  measure  of  reflective  degradation  in  space. 

The  data  obtained  from  the  ref lectometer  experiment  will  be  directly 
and  immediately  applicable  to  solar-concentrators  and  related  power  sys¬ 
tem  designs  and  to  all  optical  systems  depend  t  upon  reflective  surfaces 
in  space.  Examples  of  types  of  systems  and  instruments  which  will  benefit 
greatly  from  the  data  gathered  by  the  reflectometer  experiment  are:  solar- 
energy  converters,  space  telescopes,  optical  communication  systems — 
including  laser  systems--and  specularly  reflective  thermal  control  surfaces. 

Although  the  minimum  desired  lifetime  of  the  reflectometer  experiment  is 
one  year,  every  effort  will  be  made  to  assure  an  instrument  life  of  three 
years,  the  expected  life  of  ATS-C,  the  space  platform  for  the  reflectom¬ 
eter  experiment. 
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EXPERIMENT  OBJECTIVES 


The  objective  of  the  reflectometer  experiment  is  to  test  in  a  space 
environment,  for  an  extended  period  (one  to  three  years)  the  durability 
of  specularly  reflective  surface  materials.  The  experiment  is  designed 
primarily  to  obtain  engineering  data  necessary  for  spacecraft  optical 
systems . 

A  secondary  aim  is  to  obtain  data  of  scientific  value  relating  to  the 
degradation  process.  However,  this  objective  is  pursued  only  to  the 
extent  that  the  primary  objective  is  not  compromised.  For  this  reason, 
no  direct  systematic  measurement  of  the  environment  is  planned. 

Ideally  it  would  be  desirable  to  make  the  acquisition  of  scientific  data 
the  primary  aim  in  order  to  gain  a  thorough  scientific  understanding  of  the 
space  environment.  From  this  more  complete  set  of  test  data,  it  would 
be  possible  to  deduce  engineering  data  for  any  application.  But  fulfill¬ 
ment  of  this  objective  would  require  a  large  number  of  space  flights  and 
many  more  experiments.  Therefore,  the  experiment  is  primarily  aimed  at 
obtaining  much  needed  engineering  data  only. 

The  reflectometer  experiment  will  provide  data  consisting  of  measure¬ 
ments  of  both  the  magnitude  of  reflectivity  and  spectral  characteristics 
of  reflectivity;  analysis  of  data  will  provide  insight  into  the  cause  and 
extent  of  indicated  damage  to  sample  surfaces.  To  further  pinpoint  the 
cause  of  surface  erosion  an  effort  will  be  made  to  separate  the  effects 
on  various  surfaces  by  selectively  shielding  particular  samples. 

The  surface  samples  will  be  mounted  in  pairs  and  one  member  of  each  pair 
will  be  partially  protected  from  the  space  environment  by  a  clear  quartz 
crystal.  The  quartz  shield  will  effectively  block  most  of  the  charged 
particles  and  micrometeorites  except  for  a  very  short  interval  during 
which  data  is  recorded.  The  quartz  shield  will  permit  passage  of  ultra¬ 
violet  radiation,  however.  The  other  member  of  the  sample  pair  will  be 
continually  exposed  to  the  entire  space  environment. 
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The  reflectivity  of  each  of  the  samples  will  be  determined  by  sequen¬ 
tially  projecting  a  beam  of  light  from  a  tungsten  filament  light  bulb 
onto  the  sample  and  measuring  the  amount  reflected.  The  intensity  of 
the  return  beam  will  be  measured  by  a  "white"  light  photodetector,  and 
the  spectral  composition  of  the  reflected  beam  will  be  analyzed  by  a 
low  resolution  spectrometer,  measuring  the  reflectivity  in  each  of  the 
four  spectral  bands. 

A  more  complete  ref lectome ter  experiment  might  include  other  environ¬ 
mental  measurements  in  order  to  correlate  cause  and  effect  data.  How¬ 
ever,  the  complexity  of  the  space  environment  would  make  the  complete 
measurement  and  cataloging  of  all  environmental  parameters  an  expensive 
undertaking.  Thus,  the  success  of  the  ref lectometer  experiment  will 
depend  upon  correlation  with  the  data  from  other  sources  such  as  the 
previously  flown  ATS-1  and  other  satellite  and  space  probes.  In  addition, 
evaluation  of  other  data  collected  during  the  active  lifetime  of  the  ATS-C 
ref lectometer  experiment--such  as  solar  storm  activity  and  critical  on¬ 
board  temperature  measurements — will  be  extremely  valuable. 

Since  most  optical  systems  are  primarily  concerned  with  the  use  of 
specular  rather  than  diffuse  reflecting  surfaces,  no  attempt  will  be  made 
to  evaluate  diffuse  reflectivity. 

REFLECTOMETER  DESIGN 

Figure  1  is  an  overall  view  of  the  ref lectometer  and  its  location  in  the 
ATS-C  spacecraft.  The  essential  elements  of  the  ref lectometer  experiment 
are:  the  optics,  internal  light  source,  light  detectors,  test  sample  con¬ 
figuration,  the  electric  and  mechanical  subsystems,  and  the  telemetry  and 
command  interfaces. 

a .  Optics 

Figure  2  illustrates  the  ref lectometer  optical  schematic.  Radiant  energy 
from  one  of  the  two  incandescent  tungsten  lamps  is  passed  through  (or 
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Figure  2.  Reflectometer  Optical  Schematic 
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reflected  from,  depending  upon  the  lamp  in  use)  a  mechanical  shopper 
driven  at  2400  rpm  producing  a  pulse  rate  of  480  Hz.  An  aperture  stop 
limits  the  light  beam  passing  to  the  beam  splitter  lens  to  0.2  inch. 

The  beam  splitter  projects  30  percent  of  the  source  light  into  the 
telecentric  optical  system.  The  remaining  light  is  reflected  out  of 
the  system. 

The  optical  path  to  the  ref lectometer  sample  consists  of  two  0.3  inch 
diameter  lens  elements  and  two  mirrors  which  are  an  integral  part  of 
the  rotating  arm.  The  two-element  mirror  arrangement  focuses  the  light 
onto  the  sample  surface.  The  light  is  subsequently  reflected  off  the 
sample  surface  and  returned  along  the  same  mirror  and  lens  system  to  the 
beam  splitter.  Seventy  percent  of  the  light  is  dissipated  and  the  re¬ 
maining  30  percent  is  reflected  via  another  mirror  into  an  integrating 
sphere  spectrometer. 

Two  lead  sulphide  photoconductors,  one  silicon  photoconductor,  and  two 
photomultiplier  tubes  (P-M  tubes)  are  mounted  at  apertures  on  the  inte¬ 
grating  sphere.  Mounted  in  front  of  three  of  the  five  detectors  are 
different  spectral  filters  forming  a  system  capable  of  providing  a  nomi¬ 
nal  breakdown  of  the  entire  tungsten  lamp  spectrum.  The  spectrum  from 
0.3  to  0.65  microns  is  detected  by  the  two  P-M  tubes  ("seeing"  0.3  to 
0.4,  and  0.3  to  0.65  microns,  respectively).  The  0.65  to  1.2  micron 
spectral  band  is  detected  by  an  "nfiltered  silicon  detector.  The  first 
of  the  lead  sulfide  detectors  senses  the  1.2  to  3  micron  range  of  re¬ 
flected  light,  and  the  other  detector  examines  the  entire  (unfiltered) 
visible  spectrum. 

The  effective  optical  efficiencies  in  the  five  spectral  detection 
domains  of  the  ref lectometer  are  determined  by  a  large  number  of  factors, 
including  the  operating  temperature  of  the  tungsten  lamp,  the  beam 
splitter  efficiency,  the  various  lens  transmissivities,  the  mirror  chro¬ 
matic  vignetting,  and  finally  the  bandpass  characteristics  of  the  various 
filter-detector  combinations. 
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A  simple  tungsten  filament  incandescent  bulb  has  been  chosen  as  the  light 
source  for  the  ref lectometer ,  being  the  best  trade-off  between  relia¬ 
bility  and  lamp  lifetime,  and  the  desired  lamp  response  in  the  ultra¬ 
violet  range.  The  lamp  will  be  operated  at  a  voltage  as  close  as  pos¬ 
sible  to  that  of  sunlight  in  space.  Figure  3  shows  the  electromagnetic 
spectrum  of  both  the  tungsten  lamp  and  the  sun.  From  this,  it  can  be 
seen  that  the  lamp  does  provide  a  satisfactory  range,  with  the  possible 
exception  of  too  little  radiation  in  the  ultraviolet  region. 

The  usable  output,  which  extends  from  less  than  0.3  micron  into  the  3 
micron  region,  compares  favorably  with  the  solar  spectrum  in  which  98 
percent  of  the  solar  energy  is  found.  The  solar  spectrum  extends  some¬ 
what  beyond  that  of  a  lamp  operating  at  2800°K.  (The  tungsten  lamp  has 
little  output  below  0.35  micron  or  above  2.0  microns.) 

At  the  lower  end  of  the  spectrum  9  percent  of  the  solar  energy  is  below 
0.4  micron  and  1.2  percent  is  below  0.3  micron.  The  tungsten  lamp  ex¬ 
hibits  less  than  0.08  percent  of  its  total  energy  at  wavelengths  less 
than  0.4  micron  when  operated  at  2800°K.  Assuming  80  percent  of  the  2- 
watt  input  to  the  lamp  is  radiated,  about  4  milliwatts  of  light  energy 
with  wavelengths  below  0.45  micron  will  be  available  to  the  optical 
system. 

The  problem  of  tungsten  lamp  reliability  was  solved  by  using  two  espe¬ 
cially  rugged  lamps  and  by  operating  them  below  rated  voltage.  The  lamps 
are  low-voltage,  high-intensity  devices  which  operate  with  0.65  amp  at 
3.2  volts.  Should  the  primary  lamp  fail,  the  backup  lamp  can  be  commanded 
into  action  from  the  ground. 

The  partial  aluminum  coating  on  a  Vycor  substrate  of  the  beam  splitter 
gives  it  nearly  neutral  reflectance-transmittance  characteristics. 

Two  sets  of  ghosts  are  formed  by  the  beam  splitter.  One  set  is  formed 
by  the  outgoing  bundle  and  the  other  set  by  the  return  bundle.  The 
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Plgure  3.  Spectrum  of  the  Tungsten  Lamp 
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ghosts  formed  by  the  outgoing  rays  are  vignetted  mostly  by  the  optical 
system  and  the  ghosts  formed  by  the  return  rays  are  formed  near  the 
primary  image  at  the  integrating  sphere  entrance  window. 

To  reduce  the  sensitivity  to  sample  misalignment  the  optics  are  designed 
to  focus  the  filament  image  on  the  sample.  This  image  is  then  reimaged 
by  the  lens  onto  the  detector  plane  through  the  beam  divider.  Angular 
movement  of  the  sample  does  not  change  the  position  of  the  image.  In 
order  to  avoid  an  energy  change  as  the  sample  is  tilted  the  lens  diam¬ 
eter  is  substantially  larger  than  the  input  beam  diameter  of  0.2  inch. 

The  slit  through  which  the  reflected  rays  enter  the  spectrometer  is 
1.1  mm  x  0.284  mm,  large  enough  to  permit  the  beam  to  enter  regardless 
of  the  optical  aberrations  and  geometrical  uncertainties.  The  lens 
diameter  is  0.5  inch  with  a  focal  length  of  2.73  inches.  Although  most 
spectromatic  designs  are  extremely  sensitive  to  image  movement,  the  use 
of  an  integrating  sphere  allows  the  optical  beam  to  wander  in  the  aper¬ 
ture  without  affecting  the  output. 

After  penetrating  the  sphere  aperture,  the  beam  strikes  a  diffusely  MgO 
reflecting  interior  surface.  This  results  in  an  infinite  number  of  re¬ 
flections  producing  a  constant  power  density  and  spectral  content  on  any 
internal  point.  The  photodetectors,  with  their  respective  bandpass 
filters  are  located  on  the  surface  of  the  integrating  sphere  spectrometer. 

The  major  shortcoming  of  the  spectrometer  is  its  power  efficiency  which 
is  proportional  to  the  reflectance  of  the  interior  coating  and  to  the 
ratio  of  the  detector  area  to  the  total  surface  area  of  the  coating. 

The  spectrometer's  low  power  efficiency  requires  high  light  power  and 
highly  sensitive  photodetectors. 

The  system  uses  three  types  of  light  detectors.  One  type  is  a  lead 
sulfide  photoconductive  device  whose  electrical  resistance  is  inversely 
proportional  to  the  impinging  light  level.  PbS  detectors  have  a  broad 
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response  spectrum  extending  from  0.3  micron  to  about  3.0  microns.  Their 
response  falls  off  rapidly  at  the  near-ultraviolet  wavelengths.  The  PbS 
detector  responsivity,  on  the  order  of  100  microvolts  per  microwatt,  is 
an  order  of  magnitude  better  than  other  solid-state  or  photoconductive 
devices.  They  also  operate  at  a  fairly  low  voltage  dc  bias. 

Lead  "sulfide  photodetectors  are  extremely  temperature-sensitive,  having 
a  dc  resistivity  temperature  change  on  the  order  of  A  percent  per  degree 
centigrade.  To  reduce  the  problems  caused  by  extreme  temperature  sensi¬ 
tivity,  a  combination  of  temperature  monitoring  and  gain  sensitivity 
compensation  is  used. 

Another  undesirable  characteristic  of  lead  sulfide  is  that  its  respon¬ 
sivity  diminishes  with  age.  To  compensate  for  this  and  changes  due  to 
temperature,  the  signal  conditioning  electronics  include  ground- 
controlled  range  switching. 

The  second  type  of  light  detector  is  a  silicon  photodetector  whose  out¬ 
put  current  is  proportional  to  the  impinging  light  level.  The  silicon 
photodiode  responsivity  is  on  the  order  of  0.5  microampere  per  microwatt. 
This  type  photodetector  has  a  broad  response  spectrum  from  approximately 
0.32  micron  to  1.2  microns.  By  combining  the  spectral  output  of  the 
tungsten  source  and  the  spectral  response  of  the  silicon  photodetector 
the  capability  to  detect  the  0.65  to  1.2  micron  spectrum  was  established. 

In  a  manner  similar  to  that  of  the  lead  sulfide  detector  the  silicon 
photodetector  operates  at  a  relatively  low  dc  bias  level.  The  silicon 
diode  is  not  as  temperature  sensitive  as  the  lead  sulfide  detectors; 
however,  ground-controlled  range  switching  is  also  provided  in  this 
channel . 

The  third  type  of  light  detector,  a  photomultiplier  tube,  RCA  8571,  is 
used  to  monitor  the  near  ultraviolet  end  of  the  reflected  spectrum. 
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Its  responsivity  is  in  the  order  of  10^  microamps  per  microwatt.  The  PMT 
overall  spectral  response  is  much  narrower  than  that  of  the  PbS  detectors, 
but  high  in  the  near  ultraviolet  region.  However,  the  tube  does  require 
a  high-voltage  bias  of  between  500  and  1,000  volts. 

b .  Kef lectometer  Sampling 

Figure  1  illustrates  the  arrangement  of  the  samples  in  the  ref lectometer 
instrument.  Exposure  of  the  samples  to  the  space  environment  is  made 
possible  by  allowing  the  ref lectometer  sample  disc  to  protrude  through  a 
7-inch  diameter  hole  in  the  ATS-C  solar  panel.  By  means  of  the  two 
reflecting  mirrors  in  the  rotating  arm,  the  beam  from  the  tungsten  lamp 
is  directed  onto  the  surface  of  the  sample,  returns  to  the  optical  path 
and  finally  enters  the  spectrometer  where  it  is  quantitatively  examined. 

The  experiment  includes  20  samples,  10  of  which  are  unshielded,  8  of 
which  are  shielded  with  quartz  inserts,  and  2  of  which  are  reference 
samples  (one  full  scale  and  one  zero) .  The  reference  surfaces  are  care¬ 
fully  and  continuously  shielded  at  all  times  in  order  to  maintain  their 
calibrated  reflectivity. 

Since  the  number  of  surfaces  to  be  tested  was  limited,  considerable  ef¬ 
fort  was  devoted  to  selection  of  materials  currently  of  interest  to  the 
space  environment  experimenters.  Joint  conferences  between  NASA  experi¬ 
menters,  private  industrial  researchers,  and  ATS  program  management  per¬ 
sonnel  resulted  in  selection  of  materials  which,  from  most  points  of 
view,  offered  the  best  promise  of  resisting  the  erosion  forces  encountered 
in  a  space  environment. 

Appendix  I  lists  the  number  and  kind  of  materials  selected  for  reflectom- 
eter  test  surfaces.  The  reflective  materials  have  as  substrates  either 
nickel,  aluminum,  epoxy,  or  quartz.  The  reflective  surfaces  are  either 
aluminum  or  silver.  The  various  over-and-under-coatings  are  derivatives 
of  the  silicon  oxide  family.  The  samples  themselves  are  quite  small, 
three-quarters  of  an  inch  in  diameter. 
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Certain  geometrical  reactions  can  be  expected  from  thermally  induced 
structural  stresses.  To  reduce  the  number  of  uncontrolled  variables  and 
increase  the  reliability  of  the  instruments  and  the  validity  of  the  data, 
every  effort  is  made  to  prevent  geometrical  changes  in  the  samples  being 
tested.  It  is  imperative  that  sample  mounts  be  insensitive  to  distortion 
and  temperature  gradients.  To  accomplish  the  latter,  a  material  with 
maximum  conductivity  is  used.  The  sample  mount  will  have  a  spherical 
recess  at  the  center,  such  that  the  sample  can  be  adjusted  in  any  direction. 
To  minimize  the  creep  problem,  metal  inserts  are  used  under  the  mounting 
bolt  heads  to  redistribute  the  loading.  The  plate  upon  which  the  test 
samples  are  mounted  is  designed  such  that  it  is  a  rigid  section. 

The  arm  containing  two  mirrors  to  direct  the  tungsten  beam  on  each  of  the 
20  samples  rotates  at  a  rate  of  about  an  eighth  of  a  revolution  per  min¬ 
ute  sampling  each  reflective  surface  at  a  very  slow  rate.  Arm  rotation 
is  not  synchronized  with  the  100  rpm  rotation  rate  of  the  spin  stabilized 
ATS-C  spacecraft.  This  permits  data  readout  independent  of  the  sending 
rate  of  the  telemetry  system  or  the  spin  rate  of  the  spacecraft. 

c .  Mechanical  Subsystem 

Figure  A  shows  the  mechanical  subsystem  which  includes  a  housing  for  the 
optical  system  and  electronics  and  provides  the  drive  system  which  ro¬ 
tates  the  optical  scanner  arm. 

The  housing,  machined  from  a  solid  block,  provides  structural  support 
for  the  various  components  and,  as  a  sealed  container,  shields  reflectom- 
eter  components  against  the  hostile  space  environment.  Most  of  the  mech¬ 
anical  elements  are  inside  the  hermetically  sealed  enclosure  with  the 
exception  of  two  gears  and  a  pair  of  bearings.  The  gears  and  bearings 
which  serve  to  mount  the  arm  assembly  and  which  are  exposed  to  the  space 
environment  are  expected  to  operate  without  difficulty.  The  housing  also 
acts  as  a  conductive  heat  path  between  scanner  arm,  samples,  and  space¬ 
craft  to  maintain  a  minimum  temperature  gradient. 
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The  source  of  power  for  the  drive  system  is  a  brushless  24  volt  dc  motor 
operating  at  2400  rpm.  Power  is  transferred  through  a  100:1  speed- 
reducer  to  a  harmonic  drive  assembly  which  provides  an  effective  dynamic 
coupling  through  the  hermetically  sealed  housing  enclosure.  The  optical 
chopper  wheel  is  driven  at  1:1  speed  ratio  through  a  spur  gear  train. 

The  in-line  branch  of  the  miter  gear  train  drives  the  output  shaft  of  a 
180:1  speed  reducer,  which  in  turn  drives  the  scanner  arm  position  en¬ 
coder  through  a  1:1  spur  gear  train. 

The  90°  branch  of  the  miter  gear  is  directly  coupled  to  the  72:1  harmonic 
reducer.  The  harmonic  drive  provides  the  rotary  interference  between  the 
vacuum  of  space  and  the  sealed  housing.  The  output  of  the  harmonic  drive 
with  a  speed  reduction  72:1  drives  the  ref lectometer  rotary  scanner  arm 
through  a  spur  gear  train. 

In  order  to  produce  chopping  frequency  as  high  as  practicable,  the  chopper 
wheel  is  driven  at  the  motor  speed  of  2400  rpm.  Spacecraft  inertial  prop¬ 
erties  and  gyroscopic  couple  specification  considerations  have  forced  a 
chopper  wheel  design  of  essentially  a  flat  disc  with  twelve  equally  spaced 
sharp-edged  slots. 

A  parallel  path  via  a  speed  reducer  to  the  encoder  translates  mechanical 
angular  position  of  the  scanner  arm  into  a  digital  electrical  signal. 

The  digital  encoder  reads  out  the  position  of  the  scanner  arm  with  8-bit 
resolution,  allowing  detection  of  minor  local  variations  in  sample  damage. 

The  hard  vacuum  of  space,  which  at  ATS-C  synchronous  altitude  is  about 
-12 

10  torr,  poses  a  problem  for  the  moving  parts  contained  in  the  instru¬ 
ment.  The  most  susceptible  elements  in  the  mechanical  system  are  those 
moving  at  high  speed  or  under  a  heavy  load.  Thus,  as  far  as  possible, 
the  critical  moving  parts  are  located  inside  the  hermetically  sealed 
enclosure.  This  includes  the  brushless  dc  motor,  the  gear  box  speed 
reducer,  the  digital  shaft-position  encoder,  the  bulk  of  the  optics,  and 
all  the  electronics.  Dry  film  lubricants  (MoS2_Tef Ion)  are  employed  for 
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the  few  bearings  and  gears  which  operate  outside  the  hermetic  housing 
in  the  space  vacuum. 

d .  Electronics  Subsystem 

Figure  5  is  the  block  diagram  of  the  electronic  subsystem.  It  is  pur¬ 
posely  a  simple  system  to  maintain  the  high  reliability  needed  to  assure 
a  three-year  lifetime  in  the  experiment. 

There  are  three  principal  electrical  interfaces  with  the  spacecraft: 
spacecraft  power,  the  command  system,  and  digital  and  analog  telemetry 
systems . 

Bias  voltages  for  the  photodetectors,  high  voltage  for  the  photomultiplier, 
transistor  bias  for  the  signal  conditioning,  and  the  highly  regulated 
voltage  required  for  the  tungsten  lamps  are  all  supplied  by  dc-to-dc 
converters  in  the  ref lectometer . 

The  high  voltage  converters  are  adjustable  over  the  range  of  500  to  1000 
volts  in  order  to  set  the  peak  output  of  the  associated  photomultiplier 
tube  at  2  microamps.  Limiting  the  P-M  tube  anode  current  to  2  microamps 
tends  to  minimize  photomultiplier  drift  and  instability. 

Electrical  power  for  the  lamp's  3. 2  volts  dc  at  650  milliamps  each  is 
developed  by  a  single  low  voltage  dc-to-dc  converter.  The  output  of  the 
lamp  interface  circuits  is  used  to  activate  highly  reliable,  subminiature 
relays  which  in  turn  direct  power  to  the  selected  lamps.  Since  measure¬ 
ment  of  all  bias  voltages  is  essential  for  data  reduction  and  correlation, 
analog  values  of  the  bias  voltages  are  supplied  to  the  spacecraft  telem¬ 
etry  system. 

Of  the  four  conmand  channels,  one  is  used  to  start  the  motors,  two  others 
select  and  apply  power  to  the  desired  lamp,  and  a  fourth  changes  the  gain 
of  the  signal  conditioning. 
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Figure  5.  Electronics  Subsystem 
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The  brushless  dc  motor  is  operated  from  the  spacecraft  regulator  bus  at 
-24V  ±2°L  through  the  motor  control  interface  circuit. 

The  power  conversion  subsystem  of  the  ref lectometer  consists  of  two 
adjustable  high  voltage  power  supplies  and  a  single,  multiple-output, 
low  voltage  dc-to-dc  converter. 

There  are  basically  two  signal  processing  circuit  designs,  one  associated 
with  the  lead  sulfide  and  silicon  detectors  and  one  associated  with  the 
P-M  tube  detectors.  These  circuits  consist  of  an  FET  buffer  amplifier, 
an  integrated-circuit  amplifier  as  a  bandpass  filter,  and  a  signal  de¬ 
modulator  and  impedance  matching  buffer  circuit.  A  multiple-range 
feature  is  included  in  the  photodetector  signal  conditioning  channels 
to  compensate  for  the  undesirable  features  of  the  lead  sulfide  and  silicon 
detectors,  viz.  temperature  sensitivity  and  unpredictable  performance 
degradation  with  age. 

The  primary  telemetry  data  consists  of: 

a.  Five  measurements  of  reflectivity  (one  broadband  channel  covering 
the  entire  tungsten  spectrum--nominally  from  0.3  to  3  microns--and 
four  spectral  reflectivity  channels  extending  from  0.30  to  0.40, 

0.30  to  0.65,  0.65  to  1.2,  and  1.2  to  3  microns) 

b.  Digital  encoder  readout  (of  the  scanner  arm  position)  with  8-bit 
resolution 

c.  Temperature  of  lead  sulfide  detector,  necessary  to  determine 
absolute  detector  sensitivity  at  any  given  time 

d.  Record  of  photomultiplier  high  voltage  to  determine  PMT  gain 

e.  Record  of  the  tungsten  lamp  voltage,  to  determine  stability  of  the 
light  output 

Figure  6  is  an  analog  representation  of  the  telemetry  data.  With  time 
as  the  ordinate,  the  voltage  level  shown  between  sample  spaces  is  an 
analog  of  the  reflectivity  of  the  instrument  from  between  the  samples. 

It  represents  the  reflected  light  from  the  intersample  space.  Nominally 
ahd  ideally,  the  zero  reference  is  zero  volts  because  the  "zero  reference" 
represents  a  black  cavity. 
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Figure  6.  Format,  Reflectivity  Data 
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The  full-scale  (FS)  sample  although  not  necessarily  100  percent  reflec¬ 
tive,  is  known  and  calibrated.  The  system  in  turn  is  scaled  so  that  the 
signal  from  the  full-scale  reference  sample  represents  full-scale  and 
all  other  samples  are  relative  to  it.  The  test  sample  voltage  levels 
represent  typical  real  numbers  associated  with  particular  samples.  The 
reflectivity  of  a  given  sample  is  determined  by  the  equation  shown. 


R 

s 


where  R 

s 

V 

s 


test  sample  reflectivity 

applicable  test  sample  reflectivity  signal 
voltage  (volts) 


V^g  *  applicable  full-scale  sample  reflectivity 
signal  voltage  (volts) 


V  ■  applicable  zero  samples  reflectivity  signal 
voltage  (volts) 


End-to-end  system  calibration  is  performed  by  measuring  the  zero  and 
full-scale  reflectivity  samples  each  time  reflectivity  data  is  recorded. 
This  alleviates  or  controls  the  parameters  likely  to  contribute  to  system 
decay.  In  this  way  the  reflectivity  measurement  of  a  given  test  sample 
does  not  depend  upon  absolute  reflectance  measurement  but  upon  a  relative 
measurement  and  comparison  with  the  reference  samples.  The  scanner  arm 
contains  the  shield  assembly  which  protects  the  reference  samples  since 
confident  use  of  reference  samples  requires  that  the  samples  be  shielded 
from  the  space  environment . 

SUMMARY  OF  OPERATION 

The  following  is  a  description  of  the  normal  operating  cycle  of  the 
ref lectometer . 


The  first  ground  command  turns  on  the  spacecraft  voltage  regulator.  This 
is  followed  by  a  warmup  period  that  permits  dc-to-dc  converters  to  sta¬ 
bilize.  Stable  voltages  are  essential  for  proper  operation  of  the  detec¬ 
tors  and  P-M  tubes.  The  telemetry  is  then  activated  to  verify  the  scanner 
arm  position  and  critical  parameters. 
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The  drive  motor  is  then  started  and  during  the  next  eight  minutes  the 
scanner  arm  is  pivoted  at  the  center  of  the  array,  scanning  each  of  the 
20  samples  in  turn.  The  arm  projects  the  beam  of  light  from  the  tungsten 
lamp  onto  the  sample  being  tested.  The  portion  of  the  light  that  is  re¬ 
flected  returns  through  the  arm  down  into  the  body  of  the  ref lectometer 
to  the  detectors  and  electronics.  The  return  beam  is  analyzed  by  a  "white" 
light  detector  as  well  as  by  a  liw  resolution  spectrometer  which  analyzes 
the  return  radiation  in  each  of  four  spectral  bands.  The  scanner  arm 
stops  automatically  when  it  returns  to  the  initial  position.  Finally  the 
critical  parameters  are  again  checked  and  the  spacecraft  voltage  regula¬ 
tor  is  de-energized. 

The  frequency  of  operation  will  depend  upon  the  quality  of  data  obtained. 
For  the  first  month  or  two,  the  ref  lectometer  will  be  interrogated  fre¬ 
quently  to  procure  initial  trend  data  on  seise  of  the  dielectric  coatings 
which  may  be  particularly  susceptible  to  ultraviolet  damage.  Once  ini¬ 
tial  trends  have  been  established,  data  acquisition  will  be  less  frequent, 
possibly  only  once  a  month. 

The  total  actual  operating  time  of  the  ref lectometer--ovcr  a  three-year 
period,  the  lifetime  of  ATS-C--is  anticipated  to  be  no  more  than  50 
hours  total. 
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NOTES  (a)  Where  a  quantity  of  two  samples  per  Instrument  are  indicated,  one  is 
provided  with  a  quarts  shield. 

(b)  EOS  supplies  all  sample  mounts. 

(c)  Each  sample  is  one-half  overcoated  with  no  overcoating  on  the  other  half. 
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Introduction 


The  nroblem  to  be  treated  here  is  the  nature  of  the  light  curve  of  a 
tumbling,  diffusely  reflecting  cylindrical  satellite  under  various  orienta¬ 
tions  of  the  spin  axis.  Several  assumptions  are  made:  the  reflecting 
surface  obevs  the  Lambert  law,  the  cylinder  tumbles  end  over  end,  and  the 
annular  velocity  of  the  satellite  with  respect  to  the  star  field  is  small 
compared  to  the  annular  rate  of  snin.  This  latter  assumption  restricts  the 
problem  to  one  where  the  channinp  geometry  of  the  sun,  observer,  and 
satellite  has  negligible  effect  on  one  cycle  of  the  light  curve.  It  will  be 
indicated  how  the  method  of  solution  can  also  handle  the  case  of  a  nreces- 
sinp  cylinder  and  the  case  of  cylinders  with  end  nieces  of  various  simple 
shapes. 

Discussion 


The  expression  for  the  apparent  magnitude  of  a  satellite  may  be  conveniently 
written  as1 


M  -  -26.50  -  2.5  log(AYF)  ♦  5  log(R) 


where  y  is  the  albedo,  A  is  the  cross-sectional  area  of  the  cylinder  as  seen 
from  the  side,  and  R  is  the  slant  range  of  the  satellite.  The  brightness 
function  F  for  a  diffusely  reflecting  cylinder  obeying  Lambert's  law  of 
reflection  is  given  bv1  »- 


sin  v7!  sin 

(1)  F  =  - 1 - i[(w  -  8)  cos  0  ♦  sin  6] . 

4ir 


In  Equation  (1)  spherical  coordinates  arc  used  which  are  fixed  in  the 
rotating  bodv  so  that  the  symmetry  axis  points  in  the  direction  of  the 
coordinate  role.  The  sun  and  the  observer  have  polar  angles  which  are  given 
by  Vj  and  ^  while  6  measures  the  difference  in  their  azimuthal  angles. 

Among  the  three  quantities  and  0  there  are  only  two  independent 

variables  since  the  three  quantities  are  related  by  means  of  the  equation 
for  the  phase  ancle  p: 
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cor  0  -r  cos  cor  *  sin  sin  (-°o  cos  Q 


f  2) 


A  evince  of  coordinatos  from  the  soinninr  svstem  to  an  inertial  system  is 
now  made.  The  new  svstem  of  molar  coordinate*"  fPimire  11  is  chosen  to  taV« 
advantage  of  the  svmmotries  in  the  hrichtncss  function  which  vi  1 1  become 
annaront.  The  crcnt  circle  ioinine  the  sun  and  observer  directions  defines 
the  zero  azimuth  nn^le,  and  the  noint  midway  between  the  sun  and  observer 
directions  is  the  eouator.  Tn  this  new  svstem  the  nolar  anele  will  be 
called  A  and  the  azimuthal  anele  is  r.  “  • 


FTCIJRF  1 


Solvinr  for  the  hriehtness  function  in  terms  of  these  new  coordinates  vields 


F,  =  sin  0  |  s in  r|  sin  A 


cos^  -  -sin^  A(sin^  —  ♦  cos^  —  cos^  rl 
2  2  2 


(3) 


F(A,  r1  s  77l(1»“  tan*l(Fi/F2l)F2  ♦  Fj], 


4* 


Figure  2  is  a  sketch  of  the  hriphtness  function  in  the  A,  r  syster.  When  the 
cylinder's  axis  points  in  the  direction  of  the  pole  the  obiect  has  maximum 
brightness.  When  pointing  toward  the  arc  joining  the  sun  and  observer  no  light 
is  received  from  the  object.  Saddle  ooints  are  located  on  the  equator  at 
azimuths  *-n/2. 


This  problem  is  so  highly  symmetrical  that  the  brightness  function  for  one 
octant  is  repeated  for  each  of  the  other  seven  octants.  These  symmetries 
are  best  thought  of  as  two  mutually  perpendicular  planes  of  mirror  symmetry 
plus  a  point  symmetry.  A  brief  discussion  of  the  origins  of  these  symmetries 
will  shed  some  light  on  what  is  to  be  expected  in  less  simple  problems.  The 
mirror  symmetry  F(A,  T)  *  F(A,  -T)  will  result  for  any  object  with  cylindrical 
symmetry.  The  point  symmetry  F(A,  r)  »  F(n  •  A,  n  ♦  r)  results  from  the  fact 
that  a  cylinder  has  a  plane  of  mirror  symmetry  perpendicular  to  its  symmetry 
axis.  Thus  a  cylinder  with  dissimilar  end  caps  would  break  this  symmetry. 

The  other  mirror  symmetry  F(A,  r)  ■  F(w  -  A,  T)  comes  from  the  use  of  Lambert's 
law  of  reflection.  It  would  result  for  any  reflection  law  where  the  angles  of 
incidence  and  reflection  are  treated  identically  . 

For  a  tumbling  object  we  now  wish  to  know  the  location  and  brightness  of  the 
maxima  and  minima  and  their  relation  to  the  location  of  the  spin  axis.  The 
straightforward  approach  of  setting  up  the  equation  of  the  brightness  function 
in  terms  of  the  coordinates  of  the  spin  axis  and  a  parameter  measuring  the 
rotation,  setting  the  derivative  equal  to  zero,  and  solving  for  the  location 
of  the  extrema  is  disappointing.  The  transce*. dental  equation  which  must  be 
solved  is  extremely  difficult,  if  not  impossible,  to  solve.  However,  the 
reverse  problem,  that  of  finding  the  location  of  the  spin  axis  once  the 
location  of  the  maxima  are  know,  is  soluable. 

The  solution  is  effected  by  noting  that  the  great  circle  swept  out  by  the 
rotating  cylinder  is  tangent  to  the  isophotes  at  the  extrema.  There  is  one 
exception  to  this  rule,  a  minimum  which  has  zero  intensity.  First  the 
equation  for  the  brightness  function  will  be  set  up  in  terms  of  the  slope  of 
the  isophotes  at  the  point  tangency,  then  the  slope  of  the  isophotes  will  be 
calculated  directly  from  the  expression  for  the  brightness  function.  The 
resulting  expression  for  F  is  in  terms  of  a  measurable  parameter,  the  frac- 
tional  rotation  from  a  maximum  or  a  minimum. 
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The  location  of  the  extremum  is  designated  Ao,  rp.  The  preat  circle  swept 
out  by  the  tumbling  cylinder  has  a  slope  of  a  with  respect  to  the  line  of 
constant  colatitude.  The  ancle  X  measures  the  angle  of  rotation  away  from 
the  extremum.  All  of  this  mav  be  visualized  with  the  aid  of  Figure  5. 
Solving  the  relevant  spherical  triangles  allows  Equations  (3)  to  be  written 
as 


Fj  =  sin  0  |  cos  Tfi  cos  o  sin  X  ♦  sin  rp(cos  X  sin  Aq  ♦  cos  Ap  sin  o  sin  X)  | 
F2  *  cos^  -|cos  X  sin  Aq  4  cns  Vi  s*n  x  si n  o|“|sin2  £  ♦  cos^  j  cos^  fpj 


?  7 

-sin  X  cos  0 


/  .  a  0  i»  ■  2  r  \ 

Ism-  j  ♦  cos  -  sin  r 0J 


7  0 

♦2  cos*  —  sin  Tp  cos  fp  sin  X  cos  n(cos  X  sin  Ap  ♦  cos  Ap  sin  X  sin  o) 


(4) 


F(X)  -  i— [  (it  -  tan_1(F1/F2))F2  ♦  FjJ. 
4  tr 


FTHfRF  X 


The  slope  of  the  isonhotes,  0,  is  calculated  from  liquations  (3)  by  means  of 


(S) 


tan  OIF,  A)  , 

’  sin  A  ar,  \  nl 
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Fvaluatinp  (5)  at  the  extrema  eives 


(h) 


sin  To  cos  rp 

tan  ofTo,  A0)  =  - - - j  8ttF 


cos  A 
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vn  (l  - 
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-|  sin3  r0|  sin3  0  sin  A0 


2  I*  2  r  ■  2 

L  —  cos  rn  ♦  stn‘  -r 


cos-  y  cos~  i o  ♦  sin-  y) 


7  0 

COS-  y 


-sin2  Ap(cos2  y  sin2  fp  -  lj 


The  slone  ancle  from  (6)  is  to  be  used  in  Fquations  (4).  This  eives  the 
result  which  is  soupht. 

The  nuantitics  r0,  Ap,  and  o(rp,  Ap)  easilv  related  to  the  coordinates  of  the 
spin  axis  Tj,  \\.  Solvinp  the  proper  spherical  triancle  (Figure  *1)  yields 


(7)  cos  A]  ■  sin  Ap  cos  o 

tan(F|  -  Fp)  =  sin  o/(-cos  Ap  cos  a) 

Rotation  of  these  coordinates  into  richt  ascension  and  declination  can  be 
done  by  conventional  means. 


The  extension  of  the  above  method  to  the  more  complex  situation  follows.  For 
end  caps  on  the  cylinder  the  appropriate  brightness  functions  for  the  end  caps 
would  have  to  be  added  onto  the  cylinder's  brightness  function.  There  would 
be  a  constant  multiplying  the  term  for  the  end  caps  relating  the  cylinder's 
diameter  to  its  length.  For  a  precessing  cylinder  there  would  be  an  angle  of 
precession  to  be  determined.  The  great  circle  path  of  the  cylinder  axis  used 
in  Equations  (4)  would  have  to  be  replaced  by  the  path  of  the  precessing  axis. 

The  important  point  to  note  is  that  at  each  point  of  the  foregoing  analysis 
these  additional  complexities  can  be  added  without  preventing  the  method  from 
yielding  an  expression  for  Che  brightness  function  in  terms  of  the  fraction 
of  the  rotation  between  che  maxima  or  minima. 

One  simple  case  with  end  caps  can  be  done.  The  case  with  two  diffuse  or 
specular,  hemispherical  end  caps  would  only  increase  the  brightness  function 
F  by  a  constant  amount;  the  slope  equation  (b)  remains  unchanged  (interpreting 
the  r(A0,  r q)  in  it  as  the  result  without  end  caps).  The  constant  added  amount 
would  depend  upon  the  phase  angle  4  for  the  diffuse  end  caps. 

For  the  simple  case  of  the  cylinder  only  two  unknowns  are  involved  in  getting 
the  spin  axis  location.  With  precession  and  end  caps  added  an  additional 
unknowns  must  be  solved  for.  In  addition  the  product  of  albedo  and  area  of 
the  body  can  be  derived  from  the  absolute  magnitude  of  the  maxima. 

A  sample  tabulation  is  shown  for  the  case  of  the  cylinder.  For  a  given  value 
of  phase  angle  and  various  values  of  Ty,  Ay,  there  is  tabulated  the  value  of 
the  brightness  function  F(Ty,  Ay),  the  difference  in  magnitude  for  three  equally 
spaced,  intermediate  points,  and  the  location  of  the  spin  axis  Tj,  Aj .  Due  to 
the  symmetries  of  F(Ty,  Ay)  inherent  in  the  simple  case,  the  determination  of 
the  proper  octant  for  the  spin  axis  would  require  light  curves  for  several 
different  points  across  the  sky. 

In  the  sample  table  a  100  or  >100  has  been  printed  whenever  the  magnitude  of 
a  quantity  of  zero  intensity  is  involved.  For  slopes  which  are  indeterminant, 
o  was  set  equal  to  n/2.  We  note  that  a  great  variety  of  light  curves  are 
represented,  for  example  at  Ay  ■  7U*,  ry  ■  80*  the  tumbling  cylinder  would 
appear  constant  which  means  that  the  isophote  F(Ay,  Ty)  is  nearly  a  great 
circle.  There  are  symmetrical  light  curves  such  as  Ay  >0,  ry  «  anything. 

There  are  strongly  unsymmetrical  curves  such  as  Ay  *  20*,  Ty  ■  80*  or  Ay  ■  70*, 

Ty  »  10*.  It  is  noted  that  for  4  •  00*  when  the  spin  axis  lies  half  way  between 

the  observer  and  the  sun  then  Ay  ■  0,  Fy  ■  0  and  the  variation  amounts  to  .5 
magnitude. 

In  practice  the  reduction  of  light  curves  by  means  of  the  above  model  or  the 

more  complex  analogues  would  be  best  done  by  means  of  a  least-squares 

procedure.  Since  the  unknowns  in  Equations  (4)  and  (b),  Ty  and  Ay,  are  expli¬ 
citly  given,  the  least-squares  procedure  could  be  carried  out  in  a  straight¬ 
forward  manner.  This  advantage  should  carry  over  into  the  more  complex  models. 
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TABLE  I:  SAMPLE  TABULATION  FOR  CYLINDER" 


4  V 


32tt0  0 . 079465  0.077235  0.076S02 
•  1 -0.007  -0.137  -0.202 
.474  -0.309  -0.390  -0.533 
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Time-Varying  Reflection  from  Rotating  Convex  Surfaces 


A  method  is  presented  for  determining  the  time  variation  of  the 
photometric  flux  density  received  in  the  far  field  by  diffuse  Lambert 
reflection  from  the  surfaces  of  rotating  objects  of  different  shapes. 

The  surfaces  discussed  are  restricted  to  the  class  of  convex  quadric- 
surfaces  which  can  be  identified  with  coordinate  surfaces  of  systems 
of  orthogonal,  curvilinear  coordinates  in  euclidean  3- space.  The 
relative  motion  between  the  object  and  observer,  and  object  and  source 
is  described  by  rotation  with  three  degrees  of  freedom.  Formulas 
permitting  calculation  of  the  received  flux  density  by  routine  substitution 
are  presented.  The  results  for  some  of  the  more  general  quadric  surfaces 
can  be  obtained  by  simple  extension  of  the  results  for  conical  surfaces. 
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ABSTRACT 


A  satellite  photometric  signature  is  a  record  of  the 
variations  in  intensity  of  sunlight  reflected  either 
diffusely  or  specularly  from  an  orbiting  satellite. 
The  variations  in  photometric  intensity  that  arise 
from  precessing  cylindrical  satellites  are  consid¬ 
ered  in  this  paper.  Procedures  are  described  for 
extracting  satellite  attitude,  as  a  function  of  time, 
and  satellite  size,  in  terms  of  a  length-times- 
diameter  product,  from  photometric  signatures. 
Results  from  the  analysis  of  representative  rocket- 
body  signatures  are  given. 
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1.  INTRODUCTION 

a.  Usefulness  of  Optical  Data 

The  measurable  characteristics  of  sunlight  reflected  from  a  satellite  can 
provide  significant  information  concerning  the  satellite.  Instantaneous 
attitude  and  position,  the  satellite's  size,  shape,  and  exterior  coating 
properties,  may  all  be  determined  by  an  analysis  of  the  reflected  sun¬ 
light.  Optical  monitoring  of  satellites  is  relatively  inexpensive,  and  can 
be  a  valuable  supplement  to  radar  monitoring. 

b.  Description  of  Optical  Data 

Optical  data  are  acquired  through  a  telescope  that  tracks  the  satellite  in 
its  transit  across  the  sky.  The  data  can  be  either  photometric  signatures 
or  sequences  of  photographic  images  on  film.  Tracking  is  done  before 
sunrise  or  after  sunset  so  that  the  brightness  of  the  sky  background  is 
well  below  that  of  the  satellite.  The  capability  to  track  optically  extends 
well  into  the  evening  hours  and  may  last  through  most  of  the  night  if  the 
satellite  is  sufficiently  high.  For  example,  Echo  1  satellite  at  about  900 
mi  can  be  seen  throughout  the  night  at  northern  latitudes  near  the  period 
of  the  summer  solstice. 

The  photometric  signature  is  the  recorded  analog  output  of  the  changes  in 
intensity  of  the  received,  reflected  sunlight  from  the  orbiting  satellite. 
The  satellite  (unless  it  is  a  uniform  sphere)  exhibits  changes  in  apparent 
brightness  as  it  tumbles,  or  if  stabilized,  the  satellite  exhibits  changes 
in  apparent  brightness  as  it  changes  aspect  angle  with  respect  to  the  ob¬ 
server.  The  photoelectric  photometer,  when  placed  at  the  focus  of  the 
telescope,  records  these  changes  in  brightness  as  a  change  in  electrical 
output,  which  in  turn  provides  a  trace  or  analog  of  light  intensity  versus 
time.  This  record  is  calibrated  by  recording  the  photometric  intensity 
of  stars  of  known  magnitude  during  a  given  observational  period. 
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If  an  image  orthicon  camera  or  photographic  film  is  placed  at  the  tele¬ 
scope's  focus,  an  image  of  the  telescope-tracked  satellite  is  obtained. 

USAF  Avionics  Laboratory  Observatory  at  Sulphur  Grove,  near  Dayton, 
Ohio,  first  acquired  photometric  data  in  mid-1962.  Two  additional  fa¬ 
cilities  at  Cloudcroft,  New  Mexico  (Air  Force),  and  at  Hawaii  (ARPA) 
are  now  being  made  available.  These  facilities  are  described  in  Ref¬ 
erence  1.  Goodyear  Aerospace  recently  completed  a  NASA- sponsored, 

truck-mounted  portable  24-in.  telescope  for  the  photometric  examina- 
2* 

tion  of  satellites. 

Although  GAC  has  analyzed  satellite  imagery  as  well  as  the  photometric 
signature,  the  emphasis  has  been  placed  on  photometry.  This  paper  will 
be  restricted  to  the  interpretation  of  photometric  signatures.  Techniques 
will  be  presented  for  the  determination  of  attitude  and  size  from  the  pho¬ 
tometric  signature  of  rocket  bodies;  results  of  the  application  of  these 
techniques  to  representative  signatures  will  also  be  given. 

c.  Description  of  Photometric  Signatures 

The  photometric  signature  of  an  orbiting  satellite  is  made  up  of  either 
specular  or  diffuse  reflection  components,  or  a  combination  of  the  two. 
Typical  specular  and  diffuse  signatures  for  rocket  bodies  are  shown  in 
Figure  1.  Specular  (mirror -like)  reflection  from  a  given  satellite  sur¬ 
face  is  observed  only  if  the  normal  to  the  surface  bisects  the  angle  be¬ 
tween  the  sun-pointing  (S)  and  observer -pointing  (O)  vectors.  The 
proper  specular  geometry  1  >r  a  rocket  body  (cylinder)  occurs  briefly, 
twice  per  tumbling  cycle,  to  yield  for  each  cycle  a  pair  of  specular 
flashes,  as  in  the  Agena  (upper)  signature  of  Figure  1 . 


* 

Superior  numbers  in  the  text  refer  to  items  in  the  List  of  References. 
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Figure  1  -  Typical  Photometric  Signatures 
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A  diffusely-reflecting  surface  point  may  be  viewed  anywhere  along  a 
hemispheric  solid  angle  that  is  forward  of  the  surface  point.  Outside 
this  solid  angle,  the  point  will  be  in  shadow.  Variations  in  the  diffuse 
signature  (the  lower  signature  of  Figure  1)  occur  as  portions  of  the  cyl¬ 
inder’s  illuminated  surface  enter  and  leave  shadow,  and  as  the  cylinder  - 
to -observer  aspect  angle  varies. 

d.  Interpretation  of  Photometric  Signatures 

The  variations  in  the  aspect  angle  or  attitude  of  an  orbiting  satellite  pro¬ 
duce  (unless  the  satellite  is  a  uniformly  reflecting  sphere)  corresponding 
variations  in  the  amplitude  of  the  received  photometric  signature.  The 
changes  in  satellite  attitude,  which  result  whether  the  satellite  is  stable 
or  is  tumbling,  are  identified  by  matching,  in  the  least -squares  sense 
and  over  the  entire  transit,  the  satellite's  photometric  signature  with  a 
properly  generated  theoretical  equivalent.  The  attitude  history,  from 
which  the  "best-matched"  theoretical  signature  has  been  generated, 
characterizes  the  attitude  history  that  is  associated  with  the  orbiting 
satellite. 

Mathematical  models  have  been  developed  by  Goodyear  Aerospace  that 
describe  the  diffuse  and  specular  reflection  from  such  simple  shapes  as 
cylinders,  cones,  and  spheres,  and  from  plates  that  are  flat  or  have  sym¬ 
metrical  curvature.  These  models  are  expressed  as  a  function  of  the  to- 
be-determined  satellite  attitude  and  amplitude  parameters,  and  are  used 
to  generate  the  theoretical  signatures  that  are  then  matched,  as  described 
in  the  previous  paragraph,  with  the  acquired  experimental  signature.  The 
best  matches,  in  the  least  squares  sense,  determine  values  for  the  desired 
satellite  parameters. 

The  theoretically  generated  diffuse  signature  predicts  at  each  instant  of 
time  the  intensity  of  the  received  diffusely  reflected  sunlight.  However, 
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the  specularly  reflected  light  contribution,  because  it  occurs  only  briefly 
at  irregularly  spaced  instants  of  time,  is  handled  differently.  In  the  spec¬ 
ular  case,  the  satellite  parameters  are  used  to  generate  the  different  atti¬ 
tudes  of  the  vehicle  that  exist  at  the  times  of  the  received  flashes.  The 
generated  attitude  for  each  flash  is  then  matched  with  the  corresponding 
vehicle  attitude  that  is  required  by  specular  geometry  to  generate  each 
flash.  The  desired  parameter  set  is  then  determined  when  the  best 
match  has  occurred. 

The  attitude  parameter  solution  in  the  specular  case,  since  it  requires 
only  the  measurement  of  time  and  is  independent  of  amplitude,  will  be 
more  accurate  than  the  corresponding  diffuse  reflection  solution.  How¬ 
ever,  since  the  diffuse  albedo  is  more  reliably  known,  and  since  the  dif¬ 
fuse  return  represents  an  averaged  reflection  from  a  large  portion  of  the 
vehicle,  the  diffuse  signature  provides  a  more  accurate  indication  of  ve¬ 
hicle  size. 

The  size  of  the  satellite  is  determined  from  the  now  known  satellite  aspect 
geometry  as  a  function  of  the  properly  calibrated  (in  terms  of  stellar  mag¬ 
nitudes)  amplitude  of  the  photometric  signature.  In  the  case  of  the  cylin¬ 
der,  the  computed  "size"  is  in  reality  an  albedo-times-length-times- 
diameter  product.  If  *he  albedo  is  known  or  is  assumed  from  experience 
obtained  with  typical  spacecraft  coatings,  then  the  length-times-diameter 
product  can  be  determined. 

2.  MATHEMATICAL  MODEL  FOR  DESCRIPTION  OF  SATELLITE  ATTITUDE 

The  instantaneous  attitudes  of  a  tumbling  or  precessing  satellite  are  defined 

* 

by  the  orientations  of  the  satellite’s  longitudinal  axis  Lt  and  are  determined 


* 

The  axis  vectors,  L^  and  f2,  are  assumed  to  be  of  unity  magnitude. 
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over  the  entire  transit  of  the  satellite.  The  assumption  is  made  that  the 
ends  of  the  cylinder  do  not  reflect  light  and  that  during  the  time  the  satel¬ 
lite  is  being  observed  the  precession  axis  and  precession  angle  remain 

constant.  Accordingly,  all  possible  orientations  of  L.  lie  on  the  surface 

*  * 

of  a  conic  (see  Figure  2),  the  axis  of  which  is  the  precession  axis 
(with  right  ascension  R  and  declination  D  ).  The  precession  angle  is 
obtained  from  the  scalar  product  (fi,  L^),  namely  cos  1  (ft,  L^).  Two 
additional  scalars,  the  rotation  per.od  T^  and  an  initial  phase  angle  ftjr, 
are  required  to  determine  within  this  locus  of  possible  positions  of  the 
exact  position  at  any  time  t. 

A  desired  attitude  or  position  Lt  is  determined  as  follows.  The  rotation 
of  a  vector  about  an  axis  has  been  extensively  described  in  the  literature 
(Reference  3,  see  page  9).  The  position  L^,  after  a  rotation  of  ftt  about 
an  axis  ft  and  irom  a  reference  position  Lq,  is  given  by 

Lt  =  Lq  cos  ftt  +  (ft  x  Lq)  sin  ftt  +  (ft)(ft,  Lq)(1  -  cos  ftt)  .  (1) 


The  rotation  angle  ft.  is  defined,  in  terms  of  the  rotation  period  T  and  an 

*  P 

Initial  phase  angle  ftj  ,  by 


ft 


t 


(2) 


The  reference  position  Lq  lies  on  the  acceptable  locus  of  positions  of 
as  defined  by  the  given  scalar  product  (ft,  L^). 


+ 

The  axis  vectors,  and  ft,  are  assumed  to  be  of  unity  magnitude. 


Satellite  attitude  model 


Figure  2  -  Satellite  Attitude 


3.  INTERPRETATION  OF  DIFFUSE  SIGNATURES  FOR  CYLINDERS 


a.  Mathematical  Model  of  Diffuse  Signature 

The  received  illuminance  E  .,  diffusely  reflected  from  a  cylinder,  is  de- 

4  a 

scribed  by 


/A  .H  R  \ 

dec 


r 

(cos  cos  6  gt)  |  (tt  -  X^)  (cos  Xt)  +  sin  Xt 


(3) 


The  terminology  entering  into  the  amplitude  (the  first  factor  of  Equation  3) 

of  the  received  illuminance  ratio  E./E  is: 

d  c 


A^  =  diffuse  albedo  of  cylinder  reflecting  surface, 

Ec  =  illuminance  (lumens  per  unit  area)  from  the  sun,  and 


H  R  =  height  and  radius,  respectively,  of  the  cylinder, 
c  c 

The  time-dependent  orientation  angles  (0Qt,  0gt,  and  Xt),  which  determine 
the  remainder  of  Equation  3,  are  related  to  the  instantaneous  positions  of 
the  cylinder  axis  vector  L^,  of  the  sun-pointing  vector  (S^),  and  of  the  ob¬ 
server-pointing  vector  (Ot).  Physically  the  angles  are  defined  as  follows. 
In  a  coordinate  system  (Figure  3)  where  the  axis  position  is  made  to  co¬ 
incide  with  the  Z  (vertical)  axis,  the  angles  0  t  and  0g^  define  the  elevation, 
from  the  coordinate  system’s  horizontal  plane,  of  Ot  and  S^,  respectively, 
while  the  azimuthal  absolute  difference  angle  between  Ot  and  is  given  by 

xt- 

The  instantaneous  axis  position  is  available  from  Equation  1;  the  instan¬ 
taneous  positions  of  the  sun  (S.)  and  observer  (0.)  vectors  are  determined 

5  1  1 

respectively  from  tables  and  from  the  orbital  elements  of  the  satellite. 
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SIGNATURE  MODEL 


Figure  3  -  Signature  Model  in  the  Interpretation  of  Cylinder  Diffuse 
Signatures 
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b.  Matching  of  Theoretical  and  Experimental  Diffuse  Signatures 

The  variations,  that  occur  in  the  received  photometric  diffuse  signature, 
result  from  the  contribution  of  G^,  defined  from  Equation  3  as 


r  1 

Gt  =  (cos  cos  0gt)  (tt  -  Xt)  (cos  Xt)  +  sin  Xt|  . 


(4) 


The  following  amplitude  matching  procedure  permits  extraction  of  the  cyl¬ 
inder's  attitude  and  size  from  the  variations  in  the  signature. 


The  received  photometric  signature,  covering  the  entire  observed  transit 
of  the  satellite,  is  suitably  sampled  at  from  50  to  150  points;  and  the  re¬ 
sulting  sample  amplitudes  are  normalized  with  respect  to  the  amplitude 

* 

and  satellite  range  that  exist  at  a  signature  reference  point. 


The  normalized  signature  amplitudes  Jt  are  matched  with  the  theoretical 
amplitudes  Gt,  each  of  which  is  computed  from  Equation  4  at  the  time  of 
the  corresponding  sample  point  Jt-  A  least  squares  matching  of  the  am¬ 
plitudes  G^  and  is  carried  out. 

An  amplitude  J^,  inserted  so  that  the  amplitude  scales  of  G^  and  can  be 
properly  compared,  represents  the  actual  contribution  of  the  "G^"  term  at 
the  signature  reference  point. 

The  various  computed  values  of  G^  are  functions  of  the  five  attitude  model 
parameters  described  in  Item  2.  These  parameters,  in  addition  to  the  am¬ 
plitude  scalar,  J^,  are  varied  until  a  minimum  difference  between  G^  and 
J.  is  achieved.  The  values  of  the  five  attitude  parameters  and  of  the  sig- 

nature  reference  amplitude  J  ,  that  correspond  to  this  minimum,  charac- 

P 

terize  the  reflection  conditions  of  the  orbiting  satellite  cylinder,  and 


The  selection  of  the  signature  reference  point  is  arbitrary;  usually  the  maxi 
mum-valued  signature  amplitude  is  made  the  reference  point. 


53 


permit  the  computation  of  its  attitude  and  size  (in  terms  of  a  length  times 
diameter  product). 


4.  INTERPRETATION  OF  SPECULAR  SIGNATURE  FOR  CYLINDERS 


ocular  Reflection  Requirements 


Specularly  reflected  illuminance  from  a  given  surface  is  received  if  the  bi¬ 
sector  Bt  of  the  sun  and  observer  vectors  St  and  Ot  coincides  with  the  nor¬ 
mal  to  the  given  reflecting  surface.  Specular  reflection  from  a  cylinder 
thus  requires 


(Lt,  Bt)  =  0  . 


(5) 


The  geometry  of  specular  reflection  is  illustrated  in  Figure  4.  The  inter¬ 
section  of  the  "(L,  B)  =  0tt  plane  and  the  locus  of  Lt  define  for  each  tum¬ 
bling  cycle  the  two  cylinder  attitudes  (La  and  L^  in  Figure  4)  from  which 
specular  reflection  is  observed  as  a  pair  of  specular  flashes.  In  addition, 
specular  reflection  requires  that  the  "(L^,  B^)  =  0"  plane  and  the  L^  locus 
be  sufficiently  close  to  one  another  to  intersect. 


In  a  similar  manner,  the  requirements  for  specular  reflection  from  a  conic 
reflecting  surface  are  determined  (see  Figure  4). 


Extraction  of  Attitude  Information,  General  Method 

The  times  of  occurrence  of  the  specular  flashes  provide  an  excellent  de¬ 
scription  of  the  attitude  history  of  the  reflecting  satellite.  For  example, 
as  a  satellite  cylinder  moves  across  the  sky,  the  orientation  of  the  bisec¬ 
tor  (of  St  and  Ot)  vector  Bt  and  therefore  the  required  cylinder  specular 
attitudes,  La  and  L^  in  Figure  4,  are  continually  changing.  The  times  of 
specular  flash  occurrence,  which  are  irregularly  spaced  to  accommodate 
the  continually  changing  specular  attitudes,  thus  determine  the  attitude 
history  of  the  reflecting  cylinder  and  therefore  permit  the  determination 
of  the  five  attitude  parameters.. 


CYLINDER 

REFLECTION 


A 


REQUIREMENTS: 
<L  B)=  0 


CONE 

REFLECTION 


REQUIREMENTS: 

( a  =  ANGLE  OF  CONE 

REFLECTION  SURFACE) 
<L.  fi.)  =  S  I  N  a 


Figure  4  -  Specular  Reflection  Requirements 
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c .  Angular  Separation  y^  between  Required  and  Generated  Cylinder  Attitudes 

The  matching  of  the  required  (L^)  and  generated  (L^)  cylinder  attitudes 
is  accomplished  by  minimizing  the  angles  yh  that  separate  the  two  cylinder 
attitudes  that  occur  at  each  of  the  specular -flash  times.  The  generated 
attitude  (L^)  is  obtained  at  time  t^  for  a  given  set  of  five  attitude  parame- 
ters  from-thc  attitude  Equations-L-and2^._The  reference  attitude  Lq  cor¬ 
responds  to  the  time  (t  =  0)  of  the  first  flash. 


The  required  specular  flash  attitude  Lah  is  also  obtained  from  Equations  1 
and  2  upon  insertion  of  the  constraint,  (L^,  B^)  =  0. 


The  angle  y^,  separating  and  L^,  lies  along  the  surface  of  the  conic 
locus  of  and  is  determined  from  Equation  1. 

The  separation  angle  y^  is  restricted  to  the  range 


-7T  <  yh  «  7T  . 

d.  Evaluation  of  Cylinder  Attitude 

The  proper  set  of  five  attitude  parameters,  which  defines  the  actual  atti¬ 
tude  history  of  the  specularly  reflecting  cylinder,  occurs  when  the  required 
attitude  and  the  generated  attitude  Lth  coincide  at  all  the  specular 
flashes;  or  yh  is  minimum. 


5.  MINIMIZATION  OF  SQUARE  DIFFERENCE  EXPRESSION 

Minimization  of  the  square  difference  (SD)  expressions  in  the  diffuse  and 
specular  reflection  cases  is  required  to  determine  the  proper  values  of  the 
five  attitude  parameters  and  of  the  amplitude  parameter  when  its  solution 
is  applicable.  The  adopted  minimization  procedure  is  general  and  applies 
not  only  to  the  two  SD  expressions  already  developed  for  cylinders  but  can 
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be  used  with  other  expressions  of  SD  to  be  developed  for  signatures  of 
other  shapes. 

The  minimization  of  SD  is  accomplished  in  an  iterated  procedure  through 
proper  incrementing  of  the  attitude  (and  amplitude)  parameters.  A  Tay¬ 
lor's  series  approximation  of  SD  is  used  to  facilitate  the  determination  of 
the  proper  parameter  increments  at  each  iteration.  Constraints  are  in- 
troduced  in  the  procedure  to  insure  that  the  increments  selected  are  suf¬ 
ficiently  small  so  that  the  Taylor's  series  approximation  of  SD  is  reason¬ 
ably  accurate  and  converges. 

6.  SIGNATURE  ANALYSIS  RESULTS 

a.  Computation  of  Specular  Signature 

The  specular  solution  described  in  Item  4  has  been  applied  to  eight  pairs 
of  specular  flashes  from  an  Agena  rocket  body  ( 62(3a2 )  signature  that  was 
obtained  during  a  22  July  1965  orbit. 

Computer  solutions  of  the  specular  signature  require  initial  estimates  of 
the  input  parameters,  of  which  all  but  the  two  orientation  angles  of  ft  can 
be  furnished.  Ten  arbitrary  values  of  ft  that  are  randomly  spaced  within 
the  equatorial  coordinate  system  are  inserted  as  initial  estimates  of  ft. 
One  half  of  the  above  initial  estimates  of  ft  converged  to  two  near-zero 
minimums  of  SD  (Item  n  of  Figure  5)  -  Solution  A  and  Solution  B;  the 
other  half  converged  to  local  minimums  and  were  discarded. 

Other  data  are  necessary  to  resolve  the  ambiguity  between  the  two  near¬ 
zero  Solutions  A  and  B.  In  the  present  case,  Solution  B  was  determined 
to  be  correct  through  use  of  specular  flashes  from  a  conic  reflecting  sur¬ 
face  portion  of  the  rocket  body.  The  spacing  between  these  flashes  pro¬ 
vided  an  approximate  but  independent  solution  of  ft  (see  Item  HI  of  Fig¬ 
ure  5). 
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Figure  5  -  Computation  of  Agena  Rocket  Body  (620ar2)  Specular  Signature 
(Input  Signature  Data,  USAF  Sulphur  Grove  Observatory) 
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Additional  sources  of  data  that  could  be  used  to  resolve  the  solution  am¬ 
biguity  are  available  from  the  accompanying  diffuse  reflection  portion  of 
the  vehicle's  signature  or  from  the  specular  data  from  other  orbits  of  the 
same  vehicle. 

The  last  column  in  Item  n  gives  an  estimate  of  the  average  angular  devia¬ 
tion  yh,  between  the  required  and  generated  attitudes  of  L^.  The  deviation 
of  0. 563  deg  found  for  Solution  B  is  equivalent  to  an  average  uncertainty  in 
the  times  of  the  flashes  of  about  0. 05  sec. 

b.  Diffuse  Reflection  Computer  Program  Results 

Diffuse  reflection  signatures  of  several  orbits  of  two  rocket  bodies  were 
analyzed;  simultaneous  matching  of  three  or  four  selected  cycles  from 
each  signature  was  obtained. 

The  solutions  for  the  attitude  parameters  (f 2,  L),  Rn,  D  ,  and  T  (see 

Mb  p 

Item  2)  arr  shown  in  Table  I  for  each  rocket  body.  The  resulting  value 
of  SD,  when  the  minimal  solution  is  obtained,  provides  an  estimate  of 
the  degree  of  mismatch  between  theory  and  experiment. 

The  actual  theoretical  and  experimental  sample  point  values  are  plotted 
in  Figures  6  and  7. 

c.  Error  Analysis  of  Diffuse-Reflection  Computations 

An  error  analysis  of  the  diffuse-reflection  computations  has  been  carried 
out. 

The  following  input  errors  were  considered  in  the  error  analysis:  AE^/E^, 
photometric  signature  amplitude;  ARr/Rr,  satellite  range;  At,  time  of  photo¬ 
metric  measurements;  and  AR  and  AD  the  right  ascension  and  declination 

o  o 

of  the  observer  vector.  Expressions  have  been  obtained  that  describe  the 
effects  of  these  errors  upon  the  computed  orientations  of  the  cylinder  axis, 
L^,  and  of  the  tumbling  axis  ft  and  upon  the  computed  cylinder  size  or 
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Table  I  -  Diffuse  Reflection  Computer  Program  Results 
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Figure  6  -  Comparison  of  Experimental  and  Theoretical  Diffuse  Cycles. 
Satellite  1092  (3  Dec  1966  Orbit)  (Input  Signature  Data, 
USAF  Sulphur  Grove  Observatory) 
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Figure  7  -  Comparison  of  Experimental  and  Theoretical  Diffuse  Cycles. 
Satellite  1950  (Input  Signature  Data,  USAF  Sulphur  Grove 
Observatory) 
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length  times  diameter  product.  The  orientation  inaccuracies,  AL^  and  AO, 
are  given  as  the  degrees  of  separation  between  the  corresponding  true  and 
"noisy"  vectors;  while  the  size  inaccuracy  appears  as  AS/S. 

The  above  three  errors  (AL^,  Afi,  and  AS/S)  have  been  determined  for  the 
four  diffuse- signatures  analyzed  in  Table  I,  and  are  tabulated  in  Table  n. 
The  assumed  distribution  of  input  errors  also  appears  in  Table  n. 

7.  CONCLUSIONS  AND  FUTURE  EFFORTS 

A  technique  for  analyzing  the  specularly  and  diffusely  reflecting  optical  sig¬ 
natures  of  cylindrical  satellites,  to  provide  satellite  attitude  and  size  (as  a 
length-times-diameter  product),  has  been  presented. 

Payloads  in  general  have  shapes  that  may  be  derived,  singly  or  in  combina¬ 
tion,  from  cones,  spheres,  and  plates  (curved  or  flat)  as  well  as  from  cyl¬ 
inders.  The  technique  of  matching  corresponding  experimental  and  theo¬ 
retical  data  samples  may  be  extended  to  signatures  received  from  such 
payloads.  Also,  the  analysis  of  stabilized  bodies,  in  which  the  observed 
aspect  angle  of  the  stabilized  body  is  varying  continually  because  of  the 
satellite’s  orbital  motion,  is  possible  with  this  technique.  Future  efforts 
will  be  directed  toward  obtaining  the  necessary  computer  programs,  based 
upon  the  derived  theoretical  models,  that  will  permit  the  extraction  of 
useful  information  from  the  signatures  of  such  payloads. 

Consideration  should  be  given  to  other  information  available  in  the  optical 
signature.  For  example,  the  optical  signature  can  be  used  to  monitor  the 
in-flight  variations  in  albedo  of  a  given  satellite  reflection  coating.  If  the 
dimensions  of  the  satellite  are  known,  the  albedo  of  its  reflection  coating 
may  be  determined  very  accurately  through  an  analysis  of  the  satellite's 
optical  signature.  Thus,  the  in-flight  performance  and  changes  with  ex¬ 
posure  time  of  various  satellite  coatings  may  be  determined  by  painting 
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Table  n  -  Accuracies  of  Diffuse  Computer  Program  Results 
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rocket  bodies  that  are  used  for  other  purposes  and  observing  these  bodies 
photometrically  over  periods  of  months  or  years. 
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TRANSLATION  EDITOR'S  SUMMARY 

A  method  is  described  for  determining  the  orientation  of  the  rotation 
axis  cf  a  spin-stabilized  satellite,  fitted  with  reflecting  surfaces  parallel 
to  the  axis,  from  observations  of  flashes  of  reflected  sunlight.  The  conditions 
in  which  such  flashes  may  be  observed  are  set  out,  and  some  typical  visibility 
patterns  have  been  worked  out  and  are  shown  in  the  form  of  maps . 

The  errors  to  be  expected  are  discussed  in  detail,  and  for  most  cases  the 
direction  of  the  spin  axis  may  be  determined  to  within  ±0.5°. 


67 


CONTENTS  «... 

-  Page 

1  INTOOL  -CTION  3 

1 .1  Techniques  for  determining  the  orientation  In  space  of 

artificial  satellites  3 

1 .2  Co-ordinate  systems,  symbols,  conversion  formulae  3 

2  DETERMINIM}  THE  ORIENTATION  IN  SPACE  8 

2.1  Determining  the  normal  to  the  reflecting  surface  6 

2.2  Determining  the  axis  of  rotation  9 

3  VISIBILITY  OF  LIGHT  FLASHES  10 

3.1  Reflection  cone  10 

3.2  A  sufficient  condition  for  the  avoidance  of  reflection 

flashes  1 1 

3.3  Conditions  for  the  appearance  of  light  flashes  13 

3.4  Prograsme  for  the  calculation  of  flash  curves  14 

3.5  Results  15 

4  ACCURACY  OF  THE  TECHNIQUE  17 

4.1  Errors  in  the  normal  to  the  reflecting  surface  17 

4.2  Errors  in  the  determination  of  the  axis  18 

4.3  The  normal  angle  20 

4.4  Effect  of  the  finite  diameter  of  the  sunrs  disc  23 

5  CONCLUSIONS  28 

Acknowledgments  20 

Tables  1-6  29-35 

References  3^ 

List  of  Illustrations  38 

Illustrations  Figures  1-12 

Detachable  abstract  cards 


68 


1 


THTRORTCTION 


1 . 1  Techniques  for  determining  the  orientation  In  apace  of  artificial 
sate ill tea 

For  numerous  problems  In  Bpace  research,  It  is  necessary  to  control 
actively,  or  at  least  to  establish  passively,  the  orientation  in  space  of  the 
satellites  used.  In  the  case  of  large  astronomical  satellites  such  as  the 
American  Orbiting  Astronomical  Observatory,  it  is  aimed  to  achieve,  where  opti¬ 
cal  orientation  in  relation  to  reference  stars  is  concerned,  an  active  control 
of  the  position  of  the  space  vehicle  and  a  stability  of  the  adjustment  of  the 
instruments  to  within  ±0.1  second  of  arc.  Unlike  these  heavy  satellites,  with 
their  complicated  position-control  systems,  small  satellites  are  spin-stabilized 
by  rotation  about  their  axis  with  the  greatest  moment  of  inertia.  As  a  result 
of  the  forces  exerted  on  the  satellites  by  the  earth *s  magnetic  field  and  the 
gradient  of  the  gravitational  field,  the  axis  of  rotation  precesses  slowly.  For 

the  satellites  Explorer  IV  and  Explorer  XI,  which,  like  propellers,  rotate  about 

11  12 

a  lateral  axis,  the  direction  of  the  axis  changes  by  about  10°  per  day  *  .  On 
the  other  hand,  in  the  case  of  satellites  Tiros  I  (Pef.1),  Explorer  VTII  (Ref. 13), 
Telstar  I  (Ref .9)  and  Explorer  VII  (Ref. 13),  the  average  precession  rate  of  the 
rotation  axis  is  only  about  3°  (Tiros)  to  0.2 v  (Explorer  VII)  per  day.  The 
satellites  last  mentioned  are  so  constructed  and  balanced  that  their  axis  of 
symmetry,  which  is  also  the  axis  of  the  greatest  moment  of  inertia,  affords  a 
stable,  nutation-free  (nutation  damping)  rotation.  The  spatial  orientation  of 
such  a  satellite  is  determined  by  the  direction  of  the  axis  of  rotation  end  by 
the  angle  of  rotation  in  relation  to  r.  fixed  reference  direction  in  the  plane  of 
rotation.  The  following  devices  are  frequently  used  to  obtain  information 
concerning  these  factors: 

(a)  Horizon  sensors 

These  devices  respond  to  the  difference  in  radiation  intensity  between 

the  background  of  the  sky  and  the  earth *s  surface  (or  atmosphere).  They  rotate 

with  the  satellite  and  indicate,  through  a  telemetry  system,  the  moment  at 

which  their  line  of  sight  points  to  the  boundary  line  between  earth  and  sky. 

This  direction  cen,  in  practice,  be  measured  to  an  accuracy  of  i2°;  however, 

10  7 

an  accuracy  of  ±0.1°  -  0.2°  is  obtainable  in  principle  *  . 

(b)  Solar-aspect  measuring  devices 

These  devices  work  on  a  photoelectric  principle,  and  give  the  angle 
between  the  direction  of  the  sun  and  a  reference  direction  of  the  satellite. 

Using  this  principle,  the  accuracy  of  angular  measurements  is,  for  Explorer  VTII, 
about  ±1°  (Ref. 13)  and  for  Telstar  I,  about  ±0.5°  (Ref .9). 

(9 


(c)  Television  cameras 


With  satellites,  such  as  weather  satellites,  which  carry  television  cameras 
onboard,  it  ic  possible  to  deduce  the  position  of  the  camera  or  satellite  axis 
from  the  positions  of  known  geographical  points  on  the  television  picture.  The 

li 

accuracy  of  the  estimate  of  this  direction  for  Tiros  I  is  a  few  degrees  . 

(d)  Antennae  with  directional  characteristics 

If  the  satellite  has  a  transmitting  antenna  with  ron-isotropic  character¬ 
istics,  field  strength  variations  occur  at  the  receiver  which  have  a  period 
corresponding  to  the  natural  rotation  of  the  satellite.  As  a  result  of  this 
phenomenon,  Naumann  was  able  to  estimate  the  rotation  axis  of  Explorer  IV 
(rotating  dipole)  co  within  i20°  (Ref.n)  and  that  of  Explorer  XI  to  within 
i3°  (Ref. 12). 

(e)  Reflecting  surfaces 

If  the  sun’s  light  is  reflected  by  mirrors  on  the  satellite  to  an  observ¬ 
ing  station,  information  can  be  gained  regarding  the  position  of  the  normals  to 
the  reflecting  surface  from  the  flash  of  light  reflected;  on  an  astronomical 
scale  of  brightness,  this  reflected  light  is  several  times  greater  than  the 
light  reaching  the  observer  in  other  ways.  The  present  paper  deals  with  a 
technique  based  on  this  principle. 

The  possibility  of  estimating  the  orientation  of  a  satellite  using  the 

reflection  law  has  already  been  investigated  (in  1957)  by  Davis  and  others^. 

!Rieir  technique,  however,  was  at  first  applied  only  to  satellites  of  cylindrical 

shape,  whose  rotation  is  not  discernible  and  which  have  a  reflective  coating, 

and  to  cases  involving  practically  simultaneous  observation  by  two  stations. 

The  problem  of  estimating  the  axis  of  rotation  of  such  a  satellite  by  rotation 

19 

about  an  axis  perpendicular  to  the  axis  of  symmetry  was  solved  by  Zessewich  , 
also  by  Notni  and  Oleak1^  for  the  booster  rocket  of  Sputnik  III  (accuracy  5°). 

A  system  with  reflectors,  whose  normals  are  inclined  at  68°  and  95°  to  the  axis 
of  symmetry,  was  recently  satisfactorily  developed  (accuracy  >0.5°)  for  the 
canmunication  satellite  Telstar  (Hill  19^5,  Ref  .8).  With  this  system,  the 
position  of  the  axis  of  rotation  can  be  calculated  from  the  times  of  two 
flashes  of  reflected  sunlight  as  sections  of  two  cones.  Working  independently, 
Giese,  at  the  Smithsonian  Astrophysical  Observatory  (Cambridge,  Mass.),  during 
investigations  into  the  possibility  of  estimating  the  spatial  orientation  from 
Baker-Nhnn  camera  photographs  or  from  visual  observations^  found  a  similar, 
but  simpler,  technique^.  If  one  furnishes  a  spin-stabilized  satellite  with  one 
or  more  mirrors,  in  such  a  way,  that  the  normal  to  the  mirrors  lies  perpendicu¬ 
lar  to  the  axis  of  rotation,  the  latter  can  be  determined  as  lines  of  intersection 
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of  two  planes,  providing  that  the  observed  position  of  two  (suitable)  reflec¬ 
tion  flashes  in  relation  to  the  stellar  background  is  known.  The  times  of  the 
observations  need  be  known  only  approximate  ly  and,  depending  on  the  precession 
velocity  of  the  axis  of  the  satellite,  may  be  separated  by  hours  or  days.  The 
technique  will  be  described  again  briefly  in  Section  2.  The  main  point  of  the 
investigation  here  however,  concerns  the  positions  of  the  sun,  satellite  orb5t 
and  satellite  axis  at  which,  in  general,  light  flashes  are  to  be  expected  when 
these  positions  are  observed  from  the  earth  (Section  3)  the  accurc.cy  with 
which  the  position  in  space  of  the  rotation  axis  can  be  determined  (Section  4) 

under  different  geometrical  conditions.  The  estimation  of  the  angle  of  rota- 

14  5 

tion  and  the  period  of  rotation  has  already  been  dealt  with  in  other  papers  ;  , 
and  will  not  be  discussed  further  in  this  investigation. 

1 .2  Co-ordinate  systems,  symbols,  conversion  formulae 

Since  it  is  necessary  in  this  work  to  use  co-ordinate  systems  when  carry¬ 
ing  out  calculations,  let  us  first  sunnarize  the  systems  used,  the  definitions 
ard  conversion  formulae . 

Co-ordinate  systems; 

The  equatorial  system  (Fig.l) 

The  equatorial  system  is  an  inertial  cartesian  system,  rotating  clockwise 

•4  H ¥ 

in  relation  uo  the  direction  of  its  axis,  with  the  unit  vectors  ci  >  C0,  C,  . 

The  equatorial  plane  is  parallel  to  the  datum  plane  bounded  by  Clf  C^;  C. 
points  towards  a  fixed,  particular  direction  in  this  plane  (\E-T.al  equinox)  and 
Cj  is  directed  towards  the  celestial  north  pole.  A  polar  co-ordinate  system, 
with  the  direction  angles  ®  (right  ascension,  here  in  degrees  or  radians)  and 
6  (declination)  is  associated  with  the  equatorial  system.  According  to  whether 
the  origin  of  the  system  is  at  the  centre  of  the  earth,  at  the  place  of  observa¬ 
tion  or  at  the  centre  of  mass  of  the  satellite,  the  system  is  termed  a  geocentric 
or  tqpocentric  one  or  one  referred  to  the  satellite. 

The  orbital  system  (Fig. 2) 

The  orbital  system  has  as  its  datum  plane  the  orbital  plane  of  the 
satellite.  The  position  of  thi3  plane  in  relation  to  the  equatorial  plane  is 
fixed  by  the  angle  fl  (right  ascension  of  the  ascending  nede)  and  i  (inclination). 

— f  |  “♦  f  | 

The  unit  vectors  of  the  corresponding  cartesian  system  are  C^,  Cg, 

(clockwise -rotating  system),  cj  then  points  to  the  ascending  node  and  to 
the  northern  celestial  hemisphere.  In  a  manner  analogous  to  that  described 
above,  angles  of  direction  a*  and  6*,  corresponding  to  CJ,  Cg,  are 
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associated  with  the  orbital  system.  The  angle  u*  (corresponding  to  a*)  is  used 
specifically  to  indicate  the  position  of  the  satellite  in  the  orbit.  The  origin 
of  the  orbital  system  is  the  centre  of  the  earth. 


The  spin  axis  system  (see  Fig  .6) 

The  spin  axis  system  has  as  its  datum  plane  the  plane  In  which  the  natural 
rotation  of  the  satellite  occurs.  Its  axis  unit  vectors  are  C.,  C_,  C, 

v  "  1  d  * 

(cartesian  system).  C.  points  towards  the  ascending  node  between  the  celestial 

-*  *» 

equator  and  the  plane  of  rotation  of  the  satellite;  points  in  the  direction 
of  the  spin  axis.  If,  during  any  period  of  time  considered  any  precession  and 

-♦ii  -»  n  -♦  H 

nutation  of  the  spin  axis  can  be  disregarded,  then  the  system  C^,  Cg,  can, 
within  this  time  Interval,  be  considered  fixed  in  space  relative  to  the  direc¬ 
tion  of  its  axes.  The  origin  of  the  spin  axis  system  is  the  centre  of  the 
satellite. 


The  auxiliary  plane  system 

This  cartesian  system  is  fixed  by  the  unit  vectors  E^,  Eg,  E^,  and  has  as 
its  datum  plane  a  plane  normal  to  the  position  vector  from  the  centre  of  the 
earth  to  the  satellite.  points  in  the  direction  of  this  vector.  E1  has  the 
same  direction  as  that  of  the  component  of  the  axis  vector  of  the  satellite 
which  lies  in  the  datum  plane.  Further  details  are  given  in  Section  4.4. 


Definitiors 


Vectors  will  always  be  indicated  by  arrows  written  above  their  symbols . 
Uhit  vectors  will  be  written  in  large  letters,  other  vectors  in  small  ones. 
Components  of  vectors  will  be  representeu  by  the  letters  X,  Y,  Z,  for  unit 
vectors,  by  x,  y,  z,  in  other  cases.  Above  or  below  is  indicated  the  character¬ 
istic  symbol  for  the  vector  or  the  co-ordinate  system.  No  apostrophe  signifies 
an  equatorial  system,  the  symbol  *  an  orbital  system  and  the  symbol  "  a  spin 
axis  system. 


The  following  vectors  are  used  particularly  frequently  (Fig.3)s 


— ♦ 
*B 

*R 


(unit  vector  X;  value  r) 

position  vector  from  the  centre  of  the  earth  to  the  satellite; 

(unit  vector  X^;  value  p) 

vector  from  the  observer  to  the  satellite; 

(unit  vector  X^;  value  R) 

vector  from  the  centre  of  the  earth  to  the  observer; 
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XA5  c3 


Unit  vector,  from  satellite  to  sun 

Unit  vector  in  the  direction  of  the  spin  axis,  independent  of  the 
direction  of  rotation  of  the  satellite  in  the  hemisphere  facing  the 
sun; 

Unit  vector  with  the  direction  of  the  vector  of  the  angular  velocity 
of  the  spin  motion.  Depending  on  the  direction  of  rotation,  known 
from  the  launching  conditions, 


A  -  X.  or  A  = 
A 


•> 

X. 


Conversion  formulae: 


For  the  conversion  of  the  equatorial  system  into  an  orbital  one,  the 
following  equation  applies :- 


cos 


fl 


sin  fi  cos  i 


sin  n 

cos  Cl  cos  i 


0  \  /  x 


\ 


sin  i 


sin  Cl  sin  i  -  cos  0  sin  i  cos  v  \  z 

— ♦ 

The  components  of  the  vector  x  are  in  the  orbital  system  :- 


0) 


^r  cos  u’\ 

y* 

|  a 

r  sin  u' 

\  Z  1  } 

\  0  / 

(2) 


in  the  equatorial  system:- 


P 

i 

a 

r  .  (cos  u’  cos  fi  -  sin  uT  cos  i  sin  fi)' 
r  .  (cos  u'  sin  n  +  sin  u*  cos  i  cos  fi) 


(5) 


A  unit  vector  in  the  direction  defined  by  a  and  8  has  the  components:- 


(4) 
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2 

2.1 


DETERMINING  THE  ORIENTATION  IN  SPACE 


Determining  the  normal  to  the  reflecting  surface 

In  Fig. 4  let  Xg  be  the  vector  from  the  observer  to  the  satellite  and  let 
XQ  be  the  unit  vector  from  the  satellite  to  the  sun  (see  1 .2) .  If  reflection 
of  the  sunlight  is  observed  in  a  reflecting  element  of  the  satellite  surface, 
then,  on  tne  basis  of  the  laws  of  reflection,  a  vector  n  at  right  angles  to  the 
reflecting  surface  can  be  determined.  This  is  given  by:- 


n  * 


+  X, 


(5) 


"b  '  “0 

-4  -4 

since  n  lies  in  the  plane  of  -Xg  and  XQ  and,  in  addition,  halves  the  angle 
included  by  these  two  unit  vectors .  The  unit  vector  in  the  direction  of  n  will, 

•4 

in  future,  be  denoted  by  N. 

-* 

The  vector  X^  can,  for  the  time  of  the  observation,  be  calculated  from  the 
vector  from  the  centre  of  the  earth  to  the  sun  and  from  the  geocentric  position 
vector  x  from  the  centre  of  the  earth  to  the  satellite.  For  earth  satellites, 
however,  a  very  close  approximation  can  be  obtained  by  substituting  solely  the 
geocentric  unit  vector  to  the  nun,  given  in  astronomical  ephemerides.  The 
resulting  deviation  in  direction  for  satellites  at  a  distance  of  1.1,  2.0,  or 
10  earth-radii  from  the  centre  of  the  earth  is  only  9.6,  1 8  or  88  seconds  of 
arc  respectively.  If  cnly  the  time  t,  of  a  reflection  flash  is  recorded,  the 
unit  vector  Xg  can  be  found  from  the  known  vector  from  the  centre  of  the 

earth  to  the  observer  and  from  x(t^)  with  the  help  of  orbital  determination 
programmes:- 


x(tj)  -  ^(tt) 

x(\)  -  XrO^)! 


(6) 


It  is  much  simpler  mathematically  to  calculate  Xg  directly  from  photo¬ 
graphs  or  visual  observations  of  the  position  of  the  reflection  flash  relative 
to  known  reference  stars  in  the  sky.  From  the  observed  an, jular  directions  a 

D 

and  of  the  reflection  flash  in  the  topocentric  equatorial  system,  one  obtains 
in  cartesian  equatorisd.  co-ordinates :- 


cos  5  cos  a  v 
B  B  \ 

X„  =  I  cos  6  sin  a 
B  I  B  b 


(7) 


sin  6 


B 
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N  can  then  be  calculated  u3ing  equation  (5)  and  it  is  assumed  from  now  on 
that  this  is  known. 

2 . 2  Determining  the  axis  of  rotat'  on 

For  the  practical  estimation  of  the  rotation  axis  of'  a  satellite,  the  cases 
represented  in  Figs. 5  end  6  may  be  of  particular  iirportance.  The  first  of  these 
figures  illustrates  the  geometrical  conditions  for  the  case  of  an  elongated 
cylinder  rotating  about  its  transverse  axis  (tumbling  motion)  with  reflecting 
plane  surfaces;  the  second  figure  represents  a  squat  cylinder  rotating  about 
the  axis  of  symmetry  (spinning  motion)  which  has  a  plane  mirror  fixed  to  its 
casing.  The  first  example  corresponds  in  principle  to  the  motion  of  a  satellite 
such  as  Explorer  XI,  the  second  example  to  a  type  such  as  Tiros .  If  A  is  to  be 
calculated  for  such  a  satellite  from  reflection  flashes,  one  observation  alone 
is  not  sufficient.  The  condition  NX  A  gives  only  one  plane,  that  in  which  A 
lies .  In  order  to  determine  the  axis  of  rotation  itself,  a  second  observation 
must  be  made  in  order  to  obtain  a  further  normal  vector  "n2 *  R2),  which 
fixes  another  plane  through  the  axis  of  retption,  and  hmee  A  a3  the  intersec¬ 
tion  line  of  the  two  planes.  This  gives 


t  [%,% 

"  |E»,,  n2]| 


(see  1 .2  for  choice  of  sign) .  Of  course,  when  deciding  on  the  times  at  which 
observations  are  to  be  made,  care  should  be  taken  to  ensure  that  the  time 
interval  At  =  t2  -  t^  between  the  two  flash  observations  taken  is  sufficiently 
small  for  any  precession  of  the  spin  axis  to  be  negligible .  If  At  is  of  the 
order  of  a  few  periods  of  rotation,  thi3  provision  becomes  a  rule,  (c.f.  1.1). 


From  the  angular  directions  aBy  8^  onda^,  6^  of  the  reflection  flasles 
directly  observed  in  the  topocentric  equatorial  system,  one  obtains 


r- > 

Liy  n2j  = 

XGi 

-♦ 

C1 

-  cosSB1cosuB1 

Y01 

-♦ 

C2 

-  cosftlsinaB1 

-¥ 

C5 

Z01 

C\J 

-  cosfiggCosa^  _ 

YQ2 

-  cos  ^sina^ 

Z02 

and 

_►  ng] 

±|[n1,  n2)| 


sir.6. 

sin6 


B1 


B2 

...  (9) 
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When  only  moderate  accuracy  is  required  (±30  minutes  of  arc),  it  is 
sufficient  to  take,  from  a  map  of  the  stars,  (e.g.  Bonner  catalogue),  the 
observed  positions  of  two  appropriate  reflection  flashes,  and  to  substitute 
for  vectors  X91  and  Xgg  the  geocentric  unit  vector  to  the  sun  for  time  t^  and  tg 
given  in  ephemerides.  Since  this  vector  varies  in  direction  only  by  about 
2.5  seconds  of  arc  per  minute  and  since  t  is  not  otherwise  necessary  for  evaluat¬ 
ing  equation  (9)>  only  on  approximate  (±  a  few  minutes)  knowledge  of  the 
observation  time  is  necessary.  For  the  determination  of  A  from  measurements  of 
the  direction  of  the  reflection  flashes,  the  instrumental  expenditure  is  limited 
to  a  telescope,  a  map  of  the  stars,  a  table  of  the  sun  from  an  astronomical 
almanac  and  a  pocket  watch.  Even  the  simple  mathematical  evaluation  of 
equation  (9)  can  be  carried  out  directly  by  a  smaller  station,  such  as  that  set 
tip  by  an  expedition.  On  the  other  hand,  to  determine  A  from  equation  (6),  not 
only  is  accurate  measurement  of  the  times  t1  and  tg  respectively  of  the  reflec¬ 
tion  flashes  necessary,  but  also  the  calculation  of  the  orbit  x(t. )  must  be 
carried  out  . 

3  VISIBILITY  OF  LIGHT  FLASHES 
3.1  Reflection  cone 

For  ihe  spin  axis  system  (see  1.2)  illustrated  in  Fig  .6,  the  reflection 
condition  of  equation  (5)  gives:* 
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Thus,  for  each  illuminated  position  a"  of  the  mirror  which  rotates  with  the 

Q 

satellite,  there  is  a  corresponding  reflection  beam  with  the  angular  directions 
u-g  and  8g.  During  the  rotation  of  the  satellite,  a/'  and  hence  vary  period¬ 
ically  with  time,  whereas  5R  =  and  remain  constant.  Since  the  dimensions 
of  the  satellite  can  be  disregarded  in  comparison  with  p,  the  totality  of  all 
reflected  rays  lies  in  very  close  approximation  to  the  part  01'  the  surface  of 
a  double  cone  which  is  turned  away  from  the  sun;  the  banerator  of  this  double 
cone  forms  an  angle  with  the  plane  of  rotation  and  its  veruex  lies  at  the 

centre  of  the  satellite.  The  semi -vertex  angle  (90°  -  6”)  of  the  cone  then 

-¥  .  w  -* 

corresponds  to  the  angle  between  XA  (see  Section  1.2)  and  XQ.  This  angle  will, 
in  future,  be  referred  to  as  0  and  the  corresponding  cone  formed  from  the 
reflected  rays  will  be  called  the  "reflection  cone".  If  not  only  a  surface 
element,  but  a  zone  forming  a  band  around  the  satellite  surface  has  a  reflective 
coating,  the  sunlight  (as  in  Ref  .3)  is  simultaneously  reflected  in  the  whole 
surface  of  the  reflection  cone  (Fig. 7),  since,  at  any  given  moment,  a  reflecting 
surface  element  exists  for  every  a ^  which  satisfies  equation  (13).  In  the 
following  discussion,  this  case  will  be  considered  first.  The  results  can 
without  loss  of  generality,  be  applied  to  the  rotation  velocity  and  orbital 
motion  of  a  satellite  by  suitable  adjustment  of  the  number  of  reflecting  areas 
used  in  practice. 

3  .2  A  sufficient  condition  for  the  avoidance  of  reflection  flashes 

If  a  station  observes,  from  the  earth,  a  flash  of  light,  it  must  at  the 
momenc  of  observation,  necessarily  lie  on  the  surface  of  the  reflection  cone 
and  thus  satisfy  the  condition  :- 

cos  f>  =  X  .  X  «  SL - A  .  in) 

6  A  \y 

-♦  *4 

Since  Xg  =  x  -  x^,  writing  in  terms  of  the  components,  the  following  relation¬ 
ship  holds :- 


(x-  Xp)  XA  +  (y-  yR)  YA  (z-  zR)  ZA 
J(x-xR)^+  (y-yR)2  +  (Z'ZR)2 


(X0  XA  +  Y©  YA  +  ZQ  V 


=  0 . 
05) 


Now  positions  of  the  satellite  and  of  its  axis  of  rotation  must  be  indica¬ 
ted  for  which  the  reflection  cone  by-passes  the  earth,  and  for  which,  therefore, 
no  flashes  of  light  can  be  observed  anywhere  on  the  earth’s  surface.  Such  a 
case  is  represented  schematically  in  Fig. 8.  In  the  illustration  chosen,  the 
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axis  of  rotation  and  the  centre  of  the  earth  lie  in  the  plane  of  the  paper.  The 
axis  of  the  reflection  cone  is  shown  by  a  dotted  line  and  the  lines  of  intersec¬ 
tion  of  the  plane  of  the  paper  and  the  reflection  cone  are  drawn  thickly.  The 
sun  is  to  be  Imagined  in  the  direction  of  a  surface  line  of  the  opposed  cone  - 
not  indicated  -  associated  with  the  reflection  cone.  The  circle  represents  the 
section  of  the  earth fs  sphere  in  the  plane  of  the  paper.  As  Fig. 8  shows,  the 
conditions  under  which  the  reflection  cone  by-passes  the  earth  are  dependent  on 

■4  *♦ 

XA,  p  and  on  the  position  of  the  vertex  of  the  cone  x.  They  can  be  derived  from 
the  figure  and  are  shown  in  detail  in  Table  1  where  eight  cases,  (a)  to  (h),  are 
considered. 

*♦  *4 

If  x  and  are  specifically  expressed  in  components  of  the  orbital  system 
with  the  help  of  the  angular  directions,  finally  the  required  equation  is 
obtained: - 

cos  n  s  cos  6^  .  cos  (uf  -  a^)  .  (21) 


Fig. 9  represents  a  case  where  the  earth  is  not  intersected  by  the  reflec¬ 
tion  cone  throughout  an  entire  revolution  of  the  satellite.  For  the  sake  of 
simplicity,  the  orbital  path  is  illustrated  as  a  circle  ABCDA.  During  a  revolu¬ 
tion,  the  apex  of  the  reflection  cone  moves  with  the  satellite  along  this  path 
whereas  the  shape  of  the  cone  (0)  and  its  orientation  in  space  remain  constant. 
Firstly,  the  ratellite  is  at  point  A,  where  case  (e)  applies.  The  corresponding 
condition 


sin  (y  -  3)  >  j 


(P.2) 


for  position  A  is:- 

sin  (|  6^|  -  3)  >  J 


(23) 


If  the  satellite  moves  from  A  to  B,  as  u1  increases,  according  to 
equation  (21),  y  also  increases.  At  point  B,  u*  -  a £  =  90°  and  hence  y  =  90° 
also.  Case  (f)  applies  here,  and  the  corresponding  condition 

sin  (90°  -  3)  >  |  (24) 

is  satisfied,  if  condition  (23)  was  also  satisfied.  If  the  satellite  moves 
through  B  towards  C,  "x  .  <  0  and  y  increases  from  90°  to  (l80°  -  |  6^|  ). 

Sin  (y  -  3)  varies  from  sin  (90°  -  3)  to  sin  (|6^|  +3).  Thus,  for  the  orbital 
arc  CB, 

1  sin  (y  -  3)  >  sin  (|fi^|  -  3)  >  |  (25) 
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when  the  conditions  for  case  (b)  or  (c)  are  obtained.  Similar  considerations 
apply  also  for  the  other  half  of  the  orbit.  Condition  (23)  is  thus  a 
sufficient  criterion  for  the  absence  of  reflection  flashes  at  all  positions 
of  the  satellite. 

3  *3  Conditions  for  the  appearance  of  light  flashes 

Whereas,  for  avoidance  of  reflection  flashes,  the  condition  that  the 
reflection  cone  shall  not  intersect  the  earth  is  a  sufficient  criterion, 
equation  (15)  alone,  on  the  other  hand,  does  not  guarantee  the  appearance  of 
a  reflection  flash.  For  this,  it  is  more  important  that  four  conditions 
should  be  fulfilled  simultaneously. 

(1)  The  satellite  must  be  illuminated  by  the  sun. 

(2)  The  satellite  must  le  above  the  horizon  of  the  observer. 

(3)  The  sun  must  be  below  the  horizon  of  the  observer  (dusk  or  nignt  at  the 
place  of  observation) . 

The  fourth  condition  is  identical  with  equation  (15).  It  can  be  seen  that 
this  is  necessary  but  not  sufficient,  since  it  describes  only  the  geometry 
of  the  intersection  of  a  cone  (reflection  cone)  with  the  earth's  sphere.  If 
for  example,  this  cone  penetrates  the  earth,  then  ev~~  the  penetration  curve 
on  the  side  of  the  earth  turned  away  from  the  satellite  satifies  equation  (15); 
this  fact,  of  course,  has  no  physical  significance.  The  first  condition  i3 
illustrated  by  Fig. 10.  The  earth's  shadow  is  represented  as  a  circular 
cylinder  of  radius  R.  In  reality,  the  radius  of  the  conical  umbra  of  the  earth 
at  a  distance  of  1 .1,  2  or  10  earth  radii  from  the  centre  of  the  earth  amounts 
to  O.995,  O.99  or  O.95  R  respectively.  If  the  satellite  is  in  the  earth's 
shadow,  frcm  Fig.  10  the  following  apply  simultaneously:- 

x  .  X0  <  0  (26) 

R2  >  r2  -  (x  .  X0)2  .  (27) 

The  second  condition  is  illustrated  by  Fig.  11.  If  the  satellite  is  above  the 
horizon  of  the  observer,  then:- 

VV  0  *  (2Q) 

A  more  precise  criterion,  which  makes  it  possible  to  exclude  cases  in  which 
the  satellite,  by  reason  of  insufficient  height  above  the  horizon,  is  observed 
only  poorly,  is  the  requirement  of  a  certain  minimum  height.  From  detailed 
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considerations,  and  in  accordance  with  Veis  (19^3)^,  a  minimum  height  of  15°  is 
required  here.  Taking  the  corresponding,  largest  possible  zenith  distance 
^max  following  relationship  replaces  equation  (28):- 

‘  b >  cos  5max  •  (29) 

The  third  condition  requires  that  the  sun  should  be  below  the  horizon  at 
the  place  of  observation,  and  so  the  relationship 

(30) 

holds  (c.f.  Fig.11).  For  the  cases  to  be  investigated  more  closely  in  this 

work,  this  requirement  is  more  sharply  defined  and  extended,  in  that  the  sun, 

at  the  place  of  observation,  is  at  least  12°  below  the  horizon  (nautical 

twilight).  With  the  required  minimum  zenith  distance  I  .  (102°),  the  follow- 

min 

ing  relationship  results:- 

W~  *-»  •  (51) 

3.4  Programme  for  the  calculation  of  flash  curves 

All  points  on  the  earth’s  surface,  at  which  a  reflection  flash  can  be 
observed  at  the  same  time,  together  give  a  curve,  which  in  future  will  be 
called  the  flash  curve.  It  is  that  part  of  the  penetration  curve  of  the  earth’s 
sphere  and  the  reflection  cone  which  satisfies  the  conditions  (l)  to  (3)  given 
above.  As  the  satellite  travels  round  the  orbit,  as  a  rule  the  shape,  length 
and  geographical  position  of  the  flash  curve  change,  or  the  flash  curve 
disappears . 

For  calculating  the  flash  curves  obtained  in  real  cases,  a  machine-code 
programme  has  been  written  for  the  Siemens  2002  digital  computer.  Its  struc¬ 
ture  is  shown  in  Fig. 12  in  a  simp] ified  form.  The  input  values  are  those  of  the 
radius  (a),  the  inclination  (i),  the  ascending  node  (0)  of  the  orbit,  the 
celestial  longitude  of  the  sun  (L)  and  the  right  ascension  (a^)  and  declination 

(6.)  of  the  satellite  axis.  The  calculation  is  limited  to  circular  orbits. 

A 

After  the  6tart,  the  programme  yielded  for  inspection  the  output  data,  the 
computed  height  of  the  satellite  above  the  surface  of  the  earth  and  also  its 
period  of  revolution.  The  machine  then  calculated  cos  P  from  the  positions  of 
the  sun  and  satellite  axis  previously  found  and  set  in  the  core  store  for  the 
later  calculations  a  table  of  values  of  trigonometrical  functions,  which  are 
frequently  needed.  The  main  part  of  the  programme  calculated  the  position  of 
the  satellite  (x)  for  equally  spaced  (AT)  times  (T),  checked  whether 
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condition  (l),  (nee  3.3)  vat  satisfied  and  then,  in  the  given  case,  sought  for 

points  on  the  earth  at  which  reflection  flashes  could  be  observed.  The  sign 

of  the  function  F(x^),  (see  equation  05))  was  thus  checked  at  equally  spaced 

(AS)  latitudes  (t^),  also  for  equally  spaced  (Act)  right  ascensions  (a^ ) .  When 

the  sign  changed,  the  zero  of  F(x^)  was  found  at  the  latitude  in  question,  and 

hence  the  right  ascension  or  the  geographical  longitude,  for  which  condition  (4), 

(see  3.3)  was  satisfied.  Test  calculations  and  the  results  showed  that  the 

values  of  latitude  used  for  finding  the  zeros  could  lie  several  degrees  apart, 

without  inpairing  the  feasibility  of  the  technique  through  the  zeros  lying  too 

close  together.  For  each  zero  found,  the  prograsme  checked  that  conditions  (2) 

and  (3)  were  satisfied  and  the  printing  of  tabulated  values  was  completed 

before  the  search  was  resumed.  From  these  values,  the  geographical  latitude, 

the  right  ascension  related  to  the  geocentric  equatorial  system,  and  the 

geographical  longitude  of  the  calculated  points  of  a  flash  curve,  related  to 

a  system  rotacing  with  the  earth,  were  given  for  the  corresponding  time  and 

position  (T,  u1)  of  the  satellite.  All  the  intervals  (Ar,AJ  ,  A  a),  the  limits 

cf  the  section  of  the  satellite  orbit  (uf,  a*  ),  and  the  region  of  the  earth's 

o-  max 

surface  considered  (Km1r|,  C  ^  amin>  away)  could  be  varied  at  will.  It  was 
therefore  possible  to  carry  out  investigations  either  involving  the  whole  earth 
or*  a  small  area  of  it,  (e.g.  Europe).  Similarly,  regions  could  be  excluded  in 
which  it  was  certain  that  no  light  flashes  would  occur,  in  this  case  by 
appropriate  choice  of  tne  limits  of  the  areas  investigated. 

3 . 5  Results 

In  erder  to  give  an  idea  of  the  flash  curves  produced  during  a  revolution, 
133  examples  were  worked  out  on  the  Siemens  2002  digital  computer.  The  input 
parameters  used,  L,  a/R,  i,  0,  a-A,  and  the  information  fed  in  by  means  of  a 
code  number  for  six  different  types  are  set  out  in  Table  2.  In  most  cases  it 
was  assumed  that  the  sun  was  at  the  equinoctial  point  (L  *  0  corresponding  to 
a0  =  =  B*nce  it  cannot  at  any  time  be  more  than  23.5°  above  or  below  the 

plane  of  the  equator  and  hence  L  -  0  represents  an  average  value .  Examples 
numbers  77  to  11 5,  however,  refer  to  other  positions  of  the  sun.  In  order  to 
study  the  effect  of  the  diameter  of  the  satellite  orbit,  two  larger  groups  of 
examples,  for  a  =  1.1R  and  a  =  2R,  were  considered.  Elliptical  orbits  were 
excluded  from  the  discussion,  since  they  would  contribute  no  basically  new 
results  to  the  earlier  fundamental  research.  For  the  inclination,  i  -  60°  or 
9C°  was  chosen,  in  order  that  only  those  cases  should  be  considered  in  which 
the  satellite  could  be  satisfactorily  oDserved  from  our  latitudes .  By  a 
systematic  variation  of  ft,  it  w«*i  ensured  that  widely  different  relative 
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positions  of  the  orbital  plane  and  the  sun  vere  covered.  Finally,  the  effect  of 
different  positions  of  A  was  investigated  by  systematically  varying  and  o^ 
for  otherwise  similar  satellite  orbits  and  solar  alignments. 

The  few  examples  in  which  there  is  no  flash  can  be  made  qualitatively 
plausible  by  considering  two  limiting  cases.  One  limiting  case  is  illustrated 

*4 

in  example  number  9.  Here  coincides  with  XQ  and  only  glancing  light  is 

reflected  in  the  direction  of  -X.  onto  the  sunlight  side  of  the  earth.  The 
other  limiting  case  arises  when  X0  and  lie  in  the  plane  of  the  orbit  and  the 
reflection  cone  degenerates  to  a  plane,  since  0  =  90°.  In  this  case,  the 
satellite  is  to  be  found  directly  in  the  earth’s  shadow  if  the  degenerate 
reflection  cone  meets  the  dark  side  of  the  earth.  Since  in  examples  numbers  3, 
97  and  1l8,  XQ  lies  exactly  in  the  orbital  plane,  P  =  90°,  and  XA  forms  a  small 
angle  (15°)  with  the  orbital  plane,  it  is  understandable  that,  under  the  more 
rigidly  defined  conditions  of  visibility  (equations  (29)  and  (31 ) )  no  flashes 
could  be  observed  here . 

Some  typical  flash  curves  are  reproduced  in  Figs.  1 3  to  19.  For  purposes 
of  illustration,  the  curves,  taking  into  account  the  rotation  of  the  earth,  are 
plotted  on  a  map  of  the  earth.  The  curve  parameter  is  the  time  In  minutes 
(related  to  an  arbitrary  initial  time).  As  a  result  of  the  rotation  of  the 
earth,  the  curves  are  repeated,  after  one  revolution  of  the  satellite,  over 
other  parts  of  the  earth  at  the  corresponding  geographical  latitudes.  The 
centres  of  these  areas  are  indicated  by  circles  containing  the  number  of  the 
revolution  concerned.  It  can  be  seen  that  for  a/R  =  2,  considerably  greater 
areas  of  the  earth’s  surface  are  covered  by  flash  curves  during  a  revolution 
than  when  a/R  -  1 .1,  where  the  distance  between  the  vertex  of  the  reflection 
cone  and  the  earth’s  surface  is  relatively  small.  On  the  other  hand,  in  the 
latter  case,  the  number  of  repetitions  per  day  is  greater,  due  to  the  corres¬ 
pondingly  small  period  of  revolution  of  the  satellite .  The  fact  that  the  flash 
curves  in  Figs.  1 6  and  18  travel  appreciably  more  slowly  over  the  earth ’s  surface 
can  be  attributed  to  the  small  deviation  (15°)  of  the  angle  between  "x^  and  the 
orbital  plane  from  a  right  angle.  In  the  extreme  case,  where  this  angle  is  90° 
and  in  addition,  P  *  90°,  all  flash  curves  in  the  intersection  circle  between 
the  orbital  plane  and  the  earth  would  coincide. 
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ACCURACY  OF  THE  TECHNHUE 


4.1  Errors  In  the  normal  to  the  reflecting  surface 

The  unit  vector  in  the  direction  from  the  observer  to  the  satellite  can, 

■4  *4 

in  practice,  be  estimated  to  within  an  error  of  AXg.  Instead  of  Xg,  a  false 

unit  vector  which  forms  an  angle  with  Xg,  is  obtained.  The  normal  vector 

ri‘,  calculated  using  X'  from  equation  (5),  then  also  deviates  from  the  true 
-+  B  -♦ 
normal  r  to  the  reflecting  area.  Let  the  angle  between  n  and  nT  be 

en(°  *  cn  <  180°) .  If  n  =•  X@  -  Xg  and  n*  =  XQ 

the  equation 

Iff,  n']| 


-  (Xg  +  AXg)  are  substituted  in 


tan  e 


n 


n .  n* 


then 


tan  c 


2  .  (1  -X9  .  Xg)  -  (x9-5^)  .ax,, 


(32) 


“f  -r 

This  general  formula  includes,  moreover,  the  familiar  case  where  X«,  X^  and 

Ax_  lie  in  one  plane  and  where,  independent  of  the  angle  6  between  X  end  Xg, 

a  \  is  always  true.  It  is  known,  however,  that  in  certain  circumstances 

the  error  e  of  the  calculated  normal  direction  can  be  considerably  greater 
n 

than  the  original  direction  error  c  of  the  observation.  If,  for  example, 

^  >4  -4 

AXg  is  perpendicular  to  the  plane  defined  by  Xfi  and  Xg,  then,  since 

|^|  =•  2  sin  and  j  X@  -  Xg|  =  2  sin  |  ,  and  since  the  second  term  in  the 
denominator  of  equation  (32)  disappears,  thiw  equation  becomes: 

e. 


tan  c 


sin  -~ 


sin 


T 


(33) 


For  a  given  value  of  e  ,  e  becomes  greater  as  the  angle  between  incident  and 

a  n 

reflected  rays  becomes  more  obtuse,  i.e.  as  0  becomes  smaller.  In  the  case  of 
glancing  incidence,  *  0),  the  technique  finally  becomes  completely  unusable. 

■4 

If  Cg  is  sufficiently  small,  AXg  can  be  replaced  by  the  vector 

neglecting  terms  of  order  c^;  this  vector  lies  at  right  angles  to  Xg  and  haa 

the  value  I.e  ,  Equation  (32)  then  becomes:- 
£ 


(*q-*b>31 

4  *in2  |-Xe.a1g 


w 
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Here  the  numerator  is  always  smaller  than  2  e  and  the  denominator  is 

2  0  "pfi 

always  larger  than  4  sin  2  "  eg>  as  long  as  e  <  4  sin  This  is  always  the 

case  in  practical  measurements  for  <  12%  since  as  explained  in  3.3  on  the 

grounds  of  technical  observation  requirements  a  minimum  value  of  9, 

i  .  -  t  =  27°  is  required.  Thus,  for  the  estimation  of  the  maximum  error, 

min  max 


tan  e  S 
n 


2* 


(35) 


In  cases  where,  in  addition  to  higher  powers  of 

e  in  the  denominator  of 

B  2  0 
comparison  with  4  sin  a  more  rigorous  formula  for 

error  in  the  direction  of  n  can  br  written 


equation  (34),  XQ  .  fiXg  «  cfi 


efi,  the  linear  term  in 
can  be  neglected  in 
estimating  the  maximum 


tan  c  i  — ~  .  (36) 

sin  g 

The  error  caused  by  neglecting  the  linear  term  in  the  denominator  is  less  than 
lit  if  0  >  27°  and  efi  <  7*. 

4.2  Errors  in  the  determination  of  the  axis 

Now  we  shall  investigate  the  effect  of  errors  in  the  direction  of  the 
two  normals  used  (index  1,2)  on  the  position  of  the  axis  vector  XA,  calculated 
from  equation  (8). 

Using  the  correct  normal  vectors,  first  a  vector  lying  along  the  satellite 
axis  is  obtained: - 


xA  *  sin  7  .  XA  (37) 

“4 

where  7  is  the  angle  between  and  IL.  On  the  other  hand,  taking  the  normal 

•4  “♦  *t 

vectors  found  by  observation,  which  are  subject  to  errors,  NJ  *  N.  .  AN.  and 

«  N2  .  AN2, 

xA  =  sin  7  .  XA  +  [*,,  ^2]  +  [ANV  +  [^,  AN23  .  (38) 

"4 

The  angle  e  between  x.  and  x.  represents  the  error  in  estimating  the  axial 
A  A  A 

direction  Since 
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tan  e.  = 
A 


it  follows  that : - 


|[xA  [»v  ^31  +  CxA  %33  +  [xA  US,,  ^3?23]| 

tan  =  ^  |  |  |  )  |  t  y  (59/ 

A  sin  7  +  (XA  N1  AN2)  +  (XA  AN1  l^)  +  (XA  4^  AN^) 


If  only  small  errors  in  the  normal  directions  fle^l  «  1,  |e  |  ^  ^ ) 
are  considered,  and  can  be  replaced  by  6^  and  Sit,,  where  has  the 
value  and  is  perpendicular  to  N^.  If  &N^is  divided  into  a  component  at 
right  angles  to  the  plane  defined  by  N-j  and  Ng  and  a  parallel  component,  then: 


e1  X 
Ni  XA 


(40) 


is  a  unit  vector  at  right  angles  to  in  the  1^  1^  plane.  If  this  relation 
ship  is  introduced  into  equation  (39)  and  the  quadratic  terms  in  are 
negligible  in  comparison  with  the  linear  ones,  then:- 


tan  £, 


K  e1  “ 
"l  *!G 


-vii 


sin  7  ♦  cos  7  •  (e 


l  .T) 

re  ni' 


(41) 


Since,  in  general,  nothing  is  known  about  the  sign  and  relationship  of  the 

errora  and  ej^,  the  most  unfavourable  case  will  be  assumed  for  the  present, 

and  the  numerator  in  equation  (41 )  be  made  as  large  as  possible,  or  the 

denominator  as  small  as  possible.  In  addition,  the  c Deponents  •d.  ,  ejj.  may  be 

m  Ni 

replaced  by  the  largest  possible  error,  in  the  normal  direction.  For  the 
upper  limit  of  the  error  of  one  obtains :  - 


tan  ®A  < 


lcwl  «  feP  +  I  cos  7\) 
■in  7  -  2  .  |eN  cos  j\ 


(42) 


this  is  valid  as  long  as  sin  7  >  2je  cos  7|  .  Thus,  in  agreement  with 
intuition,  the  error  in  the  calculated  axial  direction  for  a  given  error  of 

a 

the  normal  becomes  smaller  as  7  approaches  a  right  angle,  and  greeter  as  the 
•*  «♦ 

angle  contained  by  N^  and  1^  becomes  more  acute  or  obtuse.  In  the  limiting 
case  where  7  ■  0  or  7  ■  180*,  lies  in  the  direction  or  counter-direction  of 
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and  the  technique  de scribed  for  detennlning  XA  fails  completely.  An  idea  of 
the  dependence  of  the  maximum  error  e A  of  the  axial  direction  on  7  and  can  be 
obtained  from  the  values  shown  in  Table  4,  which  were  calculated  using 
equation  (42). 

Since  the  normal  angle  is  of  great  significance  where  the  accuracy  of  the 
method  described  here  is  concerned,  in  the  following  section  formulae  for  the 
calculation  and  evaluation  of  7  will  be  derived,,  and  some  numerical  examples 
given. 

.  4.3  The  normal  angle 

The  angle  (0  *  7  *  180*)  between  the  two  normals  obtained  from  two 
observations  is  given  directly  by  the  scalar  product :- 


cos  7 


(43) 


If  one  substitutes  for  n^  and  using  equation  (5),  then,  after  multiplying  out 
the  denominator  it  follows  that:- 


cos  7 


_ ~  •  (XQ2  ~  _ 

2  .Jo  .5^)0  -4.  •  V 


m 


Thus,  for  each  pair  of  observations,  7  can  be  obtained  without  further 
difficulty. 

If  Xq,  Xg.j  and  Xjg  &re  each  divided  into  a  component  parallel  to  XA  and 
one  perpendicular  to  the  satellite  axis  (unit  vectors  N^,  N^,  K  ,  -  see  Fig. 20) 
it  follows  that:-  *  X^g  =  XQ) 


* 

sin  P  +  cos  P  1(A 

(45) 

*B1  - 

-  sin  P  cos  P  XA 

(46) 

4  " 

-  sin  P  cos  P  XA 

(47) 

From  this,  by  substituting  in  equation  (3)  one  obtains: 

nl  •  Iiil  P  (*9*  V 


(48) 
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Bln  3  (SQ  *  Kgg) 


(49) 


The  vectors  (N^  +  1^^)  and  (N^  +  Njg)  have  the  same  direction  as  the  normal 
unit  vectors  and  Ng  respectively  (Fig.21).  By  substituting  equations  (46) 
and  (49)  in  equation  (43),  an  equation  analogous  to  (44)  is  obtained:- 


cos  7  = 


(V 


V 


("a  *  V 


.  J(i  +  it0  .  ^,)(i  +  S8  .  Njj) 

If  7  is  a  right  angle  (cos  7  *  0),  it  follows*  that:- 


(50) 


*BS  -  -%1  ■  (5,) 

Equation  (51 )  means  that  observations  must  be  made  at  points  on  opposing 
generators  of  the  reflection  cone  if  7  =  90°.  If  equation  (51)  is  substituted 
In  eqviation  (461,  then,  by  adding  equations  (47)  and  (46)  a  further  form  for 
the  condition  7  =  90°  results 


Xjg  =  -  +  2  cos  3  XA  .  (52) 

As  a  rule,  one  is  not  able  to  obtain  pairs  of  observations  for  which  the 
observing  stations  lie  exactly  on  opposing  generators  of  the  reflection  cone. 
Often  the  generators  lying  opposite  each  other  at  a  particular  time  do  not  meet 
the  earth  at  all,  so  that  it  is  possible  to  combine  only  those  observations  for 
which  |  90°  -  7|  >  0.  It  is  clear  that  7  becomes  smaller  as  the  generators  of 
the  reflection  cone  corresponding  to  the  combined  observations  come  closer 
together.  In  the  limiting  case  (which  has  no  practical  significance),  where 
the  two  observations  are  made  cn  the  same  generator  of  the  reflection  cone, 

7  =  0. 

If  the  angle  between  the  generators  is  denoted  by  then 

008  r  *  *B,  -*K>  •  (55) 


*  The  cases  =  -Nq  and  =  -lfQ  are  without  practical  significance,  since 

here  the  observer  would  find  that  the  bearing  of  the  satellite  coincided  with 
that  of  the  sun. 
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In  order  to  derive  a  quantitative  relationship  between  r  and  7,  the 
following  equation  is  obtained  directly  from  Fig  .21,  which  represents  a  plane 

•4  *4  *4 

perpendicular  to  the  satellite  axis  with  the  unit  vectors  N^,  and 


7  =  \  (*2  -  V  (54) 

“4  ■+ 

since  the  normals  n^  and  n2  bisect  the  angles  ^  and  $2  respectively.  On  the 
other  hand,  ($2  -  3^)  is  the  angle  between  Ng  end  .  It  thus  follows  that:- 


cos  2  7  =  ^  .  Ngg  .  (55) 

and  3^2 

respectively,  and  X^,  using  equations  (46)  and  (47),  then  introduces  the 
relationship  given  in  equation  (53 );  and  notes,  when  evaluating  the  scalar 
product,  that  Xg  .  XA  =  cos  p,  then  the  desired  equation  is  obtained:- 

„  cos  T  -  cob  P  lrc\ 

cos  2  7  =  - = - -  .  (56) 

sin  P 

This  equation  shows  that  7,  starting  with  7*0  when  r  =  0,  increases 
with  increasing  r  and  finally,  for  T  =  2P,  reaches  the  value  7  =  90°.  The 
last  resalt  agrees  with  that  obtained  from  equation  (53);  according  to  which 
for  7  *  90°,  the  observing  stations  lie  on  opposite  generators  of  the  reflection 

■4 

cone,  i.e.  in  uhis  case,  the  angle  between  Xfil  and  X^  is  2P.  Further, 
equation  (56)  shows  that  for  a  given  value  of  T  (2P>  T),  the  angle  between  the 
normals  increases  as  the  vertex  angle  of  the  reflection  cone  decreases. 

Numerical  data,  calculated  from  equation  (56),  concerning  the  relationship 

between  r,  3,  and  7,  are  given  in  Table  5.  The  alternative  solution 

(l80°  -  7  ),  also  obtained  from  equation  (56),  gives  the  same  error  e  as 

max  a 

7 max  when  8ut>s4i4uted  in  equation  (42). 

Tne  nwirtwnan  angle  which  can  be  achieved  with  two  simultaneous 

observations,  is,  (mart  from  the  limitation  23  *  T  ) 

max 

(0  *  <,eo-> .  (57) 

Tt  is  equal  to  the  largest  anglt  at  which  two  stations  could  be  observed 
simultaneously  from  the  satellite.  This  angle  uas  a  value  of  131  60°  or 

11.5*  for  orbital  radii  of  1  1R,  2R  or  TOR  respectively.  The  corresponding 
maximum  angles  between  the  normals,  I*  y  or  180  -  7nmv,  for  P  =  45°,  7mpY  =  90°, 


If  one  resolves  %]  and  Ngg  into  components  in  the  directions  X^ 
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45*  or  8®,  and  for  P  «  90%  only  65.5%  30®  or  6®.  When  simultaneous*  observa¬ 
tions  are  made,  from  Table  4,  for  satellites  with  large  orbital  radii  (a/R  >2), 
it  can  be  expected  that  the  estimation  of  X^  will  be  insufficiently  accurate . 

The  accuracy  of  is  however  satisfactory,  even  with  distant  satellites, 
if  the  earth,  during  each  revolution,  is  intersected  twice  by  the  reflection 
ccne,  and  one  is  not  then  dependent  on  the  combination  of  "simultaneous"  or 
"almost  simultaneous"  observations.  This  case  is  illustrated  in  Table  2  and 
denoted  by  the  code  number  4  or  5.  If  the  generators  of  the  reflection  cone 
produced  by  the  first  group  of  flash  curves  form  a  sufficiently  large  angle 
with  the  generators  resulting  from  the  second  series  of  flash  curves,  observa¬ 
tions  from  two  different  series  of  flash  curves  at  any  time  can  be  combined. 
Favourable  combinations  -  leaving  aside  limitations  imposed  by  conditions  of 
visibility  are  possible  when  the  angle  the  satellite  axis  forms  with  the 
orbital  plane  is  considerably  smaller  than  P .  If  X^  lies  in  the  orbital  plane, 
a  sufficient  number  of  opposing  generators  of  the  reflection  cone  meet  the 
earth  during  one  revolution  even  when  a»  R. 

In  order  to  obtain  an  idea  of  the  values  which  7  can  attain  in  practical 
cases  by  considering  all  conditions  of  visibility,  positions  of  observing 
stations  are  indicated  on  Figs. 13  to  19*  and  7  is  calculated  for  the  combination 
of  observations  made  by  all  pairs  of  such  stations.  The  results  are  shown  in 
Table  6.  They  show,  in  agreement  with  equation  (57),  that  7  increases  as  the 
distance  between  the  stations  on  the  flash  curves  Increases,  i.e.  as  T  becomes 
greater.  Bearing  this  fact  in  mind,  the  angle  7  was  calculated  for  each  pair 
of  observations  for  00  distinct  casen  from  examples  numbers  1-76  and  116-126. 

In  57  cases,  for  a/R  =  1.1,  this  angle  was  greater  than  45®,  and  in  21  cases 
even  greater  than  75°.  On  the  other  hand,  for  a/R  =»  2,  with  combinations  from 
only  one  series  of  flash  curves,  7  remained  at  30®,  and  with  combinations  of 
observations  from  two  different  series  of  flash  curves  7  was  70®,  The  conclu¬ 
sion  that,  for  a/R  ■  2,  values  of  7  >  60®  can  only  be  attained  using  combina¬ 
tions  of  observations  from  different  sets  of  flash  curves,  agrees  with  the 
estimate  made  in  equation  (57). 

4.4  Effect  of  the  finite  diameter  of  the  sun *8  disc 

In  the  discussion  so  far,  the  sun  has  been  regarded,  by  duplication,  as  a 
point-source  of  light.  Since,  however,  the  apparent  diameter  of  the  sun’s  disc 
i3  about  0.5°,  the  flash  curves  broaden  into  bands  bounded  by  two  reflection 
cones  with  semi -vertex  angles  P  and  (P  -  0.5®).  These  cones  will,  in  future, 

*  or  an  exact  number  of  revolutions  apart. 
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be  referred  to  as  limiting  cones.  If  a  satellite  at  a  height  of  H  =  R(a/R  -  1) 
is  at  the  zenith  of  the  observer,  who  is  in  the  range  of  one  limiting  cone,  the 
distance  of  the  other  limiting  cone  from  the  observer  on  the  earth’s  surface  is 
found  to  be  b  =  H  tan  (0.5°).  For  a/R  =  1 .1  or  a/R  =  2,  the  corresponding  width 
of  the  visibility  band  on  the  earth^s  surface  is  5.6  or  55.5  km  respectively. 

If  the  satellite  moves  in  a  direction  perpendicular  to  the  tangent  to  the 

visibility  band  over  the  place  of  observation,  the  band  sweeps  over  the  observ¬ 
ing  station  in  the  shortest  possible  time.  If  the  satellite’s  velocity  is  7.5 
or  5-6  km/sec  in  the  example  taken  above,  the  approximate  duration  of  the  flash 

is  0.7  or  10  seconds.  In  both  cases,  the  length  of  the  corresponding  flash 

trail  seen  by  the  observer  is  0.5°. 

As  a  rule,  however,  the  satellite  is  not  observed  at  zenith,  nor  does  the 
visibility  band  sweep  over  the  station  in  the  way  described  above.  Therefore, 
using  the  modified  programme  described  in  Section  3.4,  two  examples  (Table  2, 
numbers  2  and  7)  were  confuted  between  K  =  47“  and  £  =  54.5°;  the  development 
of  the  visibility  bands  was  evaluated,  taking  into  account  the  earth’s  rotation 
and  plotted  on  a  map  of  Germany  (Figs. 22  and  23).  In  the  case  of  Fig. 22,  the 
width  of  the  band  is  about  12  to  13  km.  Since  the  visibility  band  for  K  “  50° 
sweeps  over  the  surface  of  the  earth  with  a  velocity  of  about  6.5  kn^sec,  a 
flash  time  of  about  2  seconds  is  obtained.  The  upper  limit  for  the  length  of 
the  flash  trail  seen  by  the  observer  can  be  found  from  this.  It  i3  equal,  to  the 
angle  which  the  satellite  would  cover  (hiring  the  period  of  a  flash  if  it  w^re  at 
the  zenith.  In  the  case  considered,  this  angle  is  about  1.4°. 

In  the  example  given  in  Fig  .23,  the  width  of  the  visibility  band  is 
approximately  the  same  as  that  for  Fig. 22.  However,  it  now  sweeps  over  the 
station  with  a  considerably  smaller  velocity  of  about  0.73  km/sec.  This  leads 
to  a  flash  duration  of  about  17  seconds.  The  upper  limit  for  the  length  of  the 
observed  flash  trail  is,  in  this  case,  circa  12°.  If  it  is  assumed  that  the 
mid-point  of  a  flash  trail  can  be  estimated  to  an  accuracy  of  ±10$  of  its  length, 
the  direction  of  can  be  estimated  to  within  ±0.14°  in  the  case  illustrated  by 
Fig. 22,  and  in  the  very  unfavourable  example  represented  by  Fig. 23,  still  to 
within  ±1.2°. 

In  order  to  study  more  generally  the  behaviour  of  the  reflection  cone  as 
it  sweeps  over  an  observing  station,  an  auxiliary  plane  will  now  be  considered. 
This  plane,  which,  as  a  rule,  intersects  the  earth’s  sphere,  contains  the  place 
of  observation  B  and,  in  addition,  has  x  as  its  normal  (Fig. 24).  Let  the  height 
of  the  satellite  above  the  auxiliary  plane  be  h.  In  the  following  discussion,  the 
co-ordinate  system  used  will  be  the  auxiliary  plane  system  defined  in  Section  1 .2. 
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Further,  let  the  angle  between  X  and  the  plane  be  denoted  by  tj.  In  the 

A 

(E^  Eg,  Ej)  system :- 


I  cos  n 
0 

\  sin  t| 


(50) 


and  If  the  corresponding  components  of  the  vector  from  the  origin  of  the 

*  »♦ 

(E.j,  Eg,  EL)  system  to  the  observing  station  are  indicated  by  x,  y  and  o, 
then:- 


(59) 


All  points  in  the  plane  which  lie  on  the  reflection  cone  must  satisfy 

r-fr 

equation  (14).  In  the  (fi^.  Eg,  F^)  system,  this  becomes:- 


cosH  .  -  X  oo.  TI  -  h  SlDJ  _  (6o) 

177/  ♦  h2 

After  squaring  and  combining  the  corresponding  equations  of  the  intersection 
cvrves  in  the  auxiliary  plane,  we  obtain:-  for  |  n|  =  3  (parabola) 


( Z)2  _  sin  2r\  (x  cos2  0  -  sin2  (A 
W  "  "  cos2  p  VH  sin  2ti  ) 


For  |t)|  >  P  (ellipse) 


(h 


\  sin  2  t) 


2  2 
cos  P  -  cos 


-J 


( 


h  Bln  23 


2  2 
cos  0  -  cos 


-)  ( 


sin  P 


I  2  2 

i cos  0  -  cos 


pJ 


(60 


(62) 
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For  |n|  <  P  (hyperbola) 


(h- 


\  sin  2n 


2  2 
cos^  n  -  cos  P 


1® 


^  h  sin  2P 


2  2 
cos  h  -  cos  P 


)2  (■ 


sin  P 


2  2 
cos  n  -  cos  P 


7 


.  (63) 


If  x  and  y  are  measured  In  units  of  h,  the  conics,  as  derived  from 
equations  (62)  and  (63),  depend  only  on  P  and  on  the  angle  n  at  which  the 
satellite  axis  is  inclined  to  the  auxiliary  plane.  In  Fig. 26  are  shown,  for 
different  values  of  n,  the  curves  resulting  from  the  intersection  of  the 
auxiliary  plane  with,  in  each  case,  two  limiting  cones  of  semi-vertex  angles 
P  n  6o°  and  P  =  59*5°.  The  geometrical  position  of  all  points  which  are  cannon 
to  the  auxiliary  plane  and  the  earth's  surface,  and  at  which  the  satellite 
appears  at  a  zenith  distance  5,  must  now  be  located  in  the  h-normalised 
auxiliary  plane  representation.  From  Fig. 24:- 


R  _  sin  <  BSP  _  b/Jb2  +  h2 
a  sin  <  SBO  “  sin  i 


(64) 


From  this  it  follows  that:- 


s 

h 


sin  5 


(65) 


The  position  desired  thus  lies  on  a  circle  about  the  origin  of  the 

. 

(E^,  E*,,  EL)  system,  which,  in  the  illustration  choren,  has  a  radius  of  s/h. 

For  i  =  0°,  s/h  is  always  0.  If  the  satellite  is  at  an  elevation  of  12°  above 
the  horizon,  for  a/R  =  1.1  or  2.0,  the  radius  s/h  =  1.9  or  0. 56.  If  the 
satellite  is  30°  above  the  horizon,  the  radius  of  the  circle  (shown  in  Fig. 26) 
is  s/h  =  1.28  or  0.48.  Cf  all  the  stations  on  such  a  circle,  only  those  which, 
in  addition,  lie  between  the  two  limiting  cones,  will  see  a  flash.  These  points 
are  to  be  found  in  the  auxiliary  plane  between  the  two  adjacent  limiting  curves 
which  correspond  to  the  parameter  in  question.  Let  the  co-ordinates  of  one  such 
station  be,  in  the  auxiliary  plane  system  (Fig. 25),  xB  =  s  cos  and 
y_  =  s  sin  Then: 

D  O 


(-  B  COS 

-  s  sin 

b 
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(66) 


In  orle  r  to  evaluate  the  length  of  the  flash  trail  seen  by  the  observer 
at  B,  the  auxiliary  plane  system  is  now  regarded  as  being  fixed  in  space,  and 
the  satellite  moves  from  the  point  (0,  0,  h)  along  its  path  in  the  direction 
indicated  by  (Fig. 25)  by  a  small  amount  d,  such  that  the  component  of  its 
movement  perpendicular  to  the  auxiliary  plane  can  be  neglected.  The  change  in 
position  of  the  satellite,  and  therefore  of  the  limiting  curves  in  the 
auxiliary  plane,  is 


d  cos 

( 

d  sin  i|r 


(67) 


The  angle  between  the  old 
satellite  is  given  by:- 


end  new  directions  from  the  observer  to  the 


^  f^jl 

•&» 


(68) 


Finally,  after  substituting  for  the  vectors  in  the  form  given  above 


tan  e . 


-  [l  +(*)2 


Bin2  (\|r. 


B 


1  + 


(hJ  ’  E  *  h  *  C0S  "  V 


(69) 


Using  this  formula,  calculations  were  made  to  find  the  interval  d  by 
which  the  satellite  can  move  in  various  directions  relative  to  the  auxiliary 
plane  before  the  length  of  the  flash  trail  observed  in  the  sky  reaches  a  value 
of  5°.  The  smallest  tolerance  values  found,  (d/h)  min,  are  indicated  by  circles 
round  B^  and  Bg  in  Fig. 26.  They  represent  the  case  where  the  satellite  moves 
perpendicular  to  the  direction  from  which  the  flash  is  observed,  so  that  no 
fore -shortening  of  the  flash  trail  occurs. 

Since,  in  the  representation  of  the  auxiliary  plane,  when  the  satellite 
moves  by  a  further  interval  of  d,  the  limiting  curves  are  also  displaced  by  a 
distance  d  in  the  auxiliary  plane,  the  range  of  the  angle  ^  for  which 
remains  below  5°  can  be  read  directly  from  Fig. 26,  If,  from  B^  or  Bg,  one 
marks  off  in  different  directions  the  circles  of  radius  (d/h)  min  shown,  the 
condition  ed<  5°  is  satisfied  for  all  radii  which  intersect  one  limiting  curve. 
If  it  is  further  assumed  that  the  mid-point  of  a  flash  trail  can  be  estimated 
to  an  accuracy  of  ±1 0$,  the  chosen  example  shows  that  the  technique  makes  it 
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-4 

possible  in  most  case’s  to  estimate  Xg  to  within  ±0.5°.  Only  a  range  of  8°  or  7° 
(about  the  tangent  to  the  limiting  curve)  is  unusable  in  the  case  of  observing 
positions  B1  or  Bg. 

These  results  are  not  significantly  impaired  if,  instead  of  a  complete 
reflecting  ring  on  the  surface  of  the  satellite  (c.f.  3*1)  several  individual 
plane  mirrors  are  used,  when  the  reflection  flash  is  cut  into  several  short 
single  flashes.  According  to  data  in  the  bibliography^1  ^ ^ the 
rotational  velocity  of  a  spin-stabilized  satellite  lies  between  6  and  0.2 
rotations  per  second.  In  order  to  obtain  about  five  single  flashes  per  second 
for  a  total  flash  duration  of  about  one  second,  between  one  and  twenty-five 
equally  spaced  mirrors  must  be  fixed  to  the  casing  of  the  satellite.  The  number 
of  mirrors  required  for  a  longer  duration  of  the  total  flash  cycle  is  correspond¬ 
ingly  smaller. 

5  OORCLUSION 

If  mirrors  are  fixed  to  a  rotating  satellite,  with  their  normals  perpendicu- 
lnr  to  the  axis  of  rotation,  the  satellite's  orientation  in  space  can  be  deter¬ 
mined  by  means  of  observed  flashes  of  reflected  sunlight.  This  technique  is 
distinguished  by  the  relatively  small  expenditure  on  computation  and  instruments . 
As  an  evaluation  of  the  possibilities  of  practical  application,  on  the  one  hand, 
conditions  are  given  under  whicn  the  method  fails.  On  the  other  hand,  it  was 
chown  that  the  technique  is  at  least  as  accurate  (+0.5°)  as  familiar  methods 
described  in  the  literature,  particularly  when  the  angle  between  the  satellite 
axis  and  the  direction  of  the  sun  is  between  45°  and  90°,  and  for  satellite  orbits 
having  semi -major  axes  of  less  than  2  earth  radii.  135  numerical  examples  were 
worked  out  and  quantitative  data  is  given  concerning  the  errors  to  be  expected 
and  their  evaluation. 
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Table  1 


SUFFICIENT  CONDITIONS  FOR  THE  ABSENCE  OF  REFLECTION  FLASHES 
(REFLECTION  CONE  DOES  NOT  INTERSECT  THE  EARTH'S  SURFACE) 

P  semi -vertex  angle  of  reflection  cone 

M  angle  between  radius  vector  x  and  axis  of  cone  X. 

-*  -*  A 
x  .  XA  projection  of  the  radius  vector  on  the  axis  of  the  cone 

R  radius  of  earth  (c.f.  Fig. 8). 


Case 

Vertex  of  cone 

Angular  range 

Conditions 

a 

x  .  X  <  -  R 
A 

0  <  0  S  90° 

no  further  conditions 

b 

-  R  *  x  *  X A  <  0 

A 

0  >  p  -  90° 

sin  -  0)  >  | 

— ♦ 

c 

-  R  *  x  .  X.  <  0 

A 

0  <  p  -  90° 

1  >  S 

p 

d 

— *  "♦  . 

x  .  X.  >  0 

A 

0  >  n 

sin  (0  -  |i)  >  § 

e 

x  .  XA  >  0 

A 

0  <  n 

sin  (n  -  0)  >  £ 

f 

X 

• 

>* 

n 

o 

M  =  90° 

cos  0  >  s 

r 

-»  -» 

g 

x  .  X  =  r 

A 

o 

u 

n. 

sin  0  >  S 

h 

-♦  -♦ 

x  .  X.  **  -  r 
A 

0  <  3  =  9C°;  p  =  T8o8 

,  >i 

r 

Note:  I  wish  to  thank  Professor  Lenz  for  the  information  that  the  areas  in 
the  plane  of  drawing,  in  which  the  vertex  of  the  reflection  cone  must  lie  if 
no  light  flash  is  to  occur,  are  defined  by  the  three  conditions : 


0) 

H  >  0  +  90° 

06) 

(2) 

0  <  pi  *  p  +  90° 

07) 

and 

sin  0*  -  P)>  S 

03) 

(3) 

0  <  n  <  0 

(•9) 

and 

sin  (0  -  n)  >  S 

(20) 

and  are  illustrated  particularly  clearly  in  Fig. 8(a). 
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Table  2 


WORKED  EXAMPLES 


a  orbital  radius,  i  inclination, 

right  ascension  of  the  spin  axis  vector, 

6^  declination  of  the  spin  axis  vector, 

L  celestial  longitude  of  the  sun, 

EL  possibilities  of  observing  reflection  flashes: 

0  =  none, 

1  =  once  per  revolution  in  the  northern  hemisphere, 

2  =  once  per  revolution  in  the  southern  hemisphere, 

3  -  once  per  revolution  in  the  vicinity  of  the  equator, 

4  =  once  in  the  northern  and  once  in  the  southern  hemisphere  per 

revolution 

5  =  twice  per  revolution  in  the  vicinity  of  the  equator 
()  =  test  carried  out  only  with  the  visibility  restriction 


max 


i 


min 


90' 


No. 

L 

a 

R 

i  i 

1 

n 

aA 

5a 

B1 

1 

0 

1.1 

60 

0 

0 

45 

1 

2 

0 

1.1 

60 

0 

45 

45 

1 

3 

0 

1.1 

60 

0 

90 

45 

0 

4 

0 

1.1 

60 

0 

135 

45 

2 

5 

0 

1.1 

60 

0 

180 

45 

2 

6 

0 

1.1 

60 

0 

225 

45 

4 

7 

0 

1.1 

60 

0 

270 

45 

4 

8 

0 

1.1 

60 

0 

315 

45 

4 

9 

0 

1.1 

60 

0 

0 

0 

0 

10 

0 

1.1 

60 

0 

45 

0 

1 

n 

0 

1.1 

60 

0 

90 

0 

4 

12 

0 

1.1 

60 

0 

135 

0 

2 

13 

0 

1.1 

60 

0 

0 

90 

0 

14 

0 

1.1 

60 

60 

0 

45 

1 

15 

0 

1 .1 

60 

60 

45 

45 

1 

l6 

0 

1.1 

60 

60 

90 

45 

2 

17 

0 

1.1 

60 

60 

135 

45 

2 

l8 

0 

1.1 

60 

60 

180 

45 

2 

19 

0 

1.1 

60 

60 

225 

45 

4 

20 

0 

1.1 

60 

60 

270 

45 

4 

21 

0 

1.1 

60 

60 

315 

45 

4 

22 

0 

1.1 

60 

60 

45 

0 

1 

23 

0 

1.1 

60 

60 

90 

0 

1  I 

24 

0 

1.1 

60 

60 

135 

0 

2 

25 

0 

1.1 

60 

120 

0 

45 

1  1 

26  1 

0  ' 

1.1 

60 

120 

45 

45 

4  ; 
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Table  2  (Contd) 


No. 

L 

a 

R 

i 

n 

aA 

6a 

— 

B1 

150 

0 

2.0 

90 

0 

135 

45 

[2) 

131 

0 

2.0 

90 

0 

225 

45 

J21 

132 

0 

2.0 

90 

0 

270 

45 

2) 

133 

0 

2.0 

90 

0 

315 

45 

!d 

Table  5 


Dependence  of  the  value  of  the  directional  error  (minutes  of  arc)  in 

the  calculated  normal  to  the  reflecting  surface  on  the  error  e  (minutes  of  arc) 

B 

in  the  observed  direction  of  the  flash  of  light  and  on  the  angle  0  (degrees) 
between  the  direction  of  the  sun  Xq  and  the  direction  from  the  observer  to  the 
satellite  (X^) .  (Using  equation  (55).) 


e° 

\ 

30° 

vn 

C 

0 

H 

110° 

130° 

O 

O 

IT\ 

o 

O 

1 

7.472 

2.800 

1.210 

1 .000 

0.752 

0.609 

0.534 

0.504 

5 

37.52 

14.03 

7.607 

5.004 

3.727 

3.045 

2.681 

2.520 

10 

75.45 

28.11 

15.23 

1 .001 

7.46c 

6.093 

5.363 

5.042 

15 

113.7 

42.24 

22.87 

15.03 

11 .20 

9.143 

8.C48 

7.566 

30 

231 .1 

85.OO 

45.90 

30.13 

22.43 

18.31 

16. 11 

15.15 

45 

351.9 

128.3 

69.17 

45.29 

33.69 

27.50 

24.20 

22.75 

60 

476.0 

172.0 

92.39 

60.52 

45.00 

36.72 

32.30 

30.36 

75 

603.2 

216.5 

115-9 

75.81 

56.34 

46.00 

40.42 

33.00 

90 

733.2 

261.0 

139.5 

91.17 

67.71 

55.22 

48.57 

45.64 

105 

865.8 

306.3 

163.2 

106.6 

79.13 

64.51 

56.73 

53-31 

120 

1000 

352.0 

187.2 

122.1 

90.57 

73.82 

64.90 

61 .00 

180 

1554 

5390 

184.5 

184.7  ■ 

136.7 

ill  .3 

97.81 

91.88 

240 

2110 

733.3 

383.6 

248.2  ' 

183.4  : 

149  2 

131.0 

123.3 

300 

933. 

485.0  ' 

312.8 

230.7  1 

187.4 

164.5 

15M 

99 


Table  4 


dependence  of  the  angular  error  (minutes  of  arc)  of  the  satellite  axis 
on  the  directional  error  (minutes  of  arc)  of  the  normals  to  te  used  in  the 
calculation  and  on  the  angle  7  (or  180°  -  7)  between  these  normals. 


\  7 

cn\ 

10® 

20° 

30° 

40® 

50° 

H 

B 

90° 

1 

11.5 

5.77 

3.87 

2.93 

2.37 

2.00 

1.74 

1.56 

l  .41 

5 

58.3 

29.0 

19.4 

14.7 

n.9 

1  .00 

1.73 

7.78 

l  .07 

10 

118 

58.5 

39.0 

29.4 

23.8 

20.1 

17.5 

15.6 

14.1 

15 

181 

88.5 

58.8 

44.3 

35.8 

30.2 

26.2 

23.4 

21 .2 

50 

381 

181 

120 

89.6 

7?.0 

60.6 

52.6 

46.8 

42.4 

45 

600 

279 

182 

136 

109 

91.4 

79.2 

70.3 

63.6 

841 

381 

246 

183 

146 

122 

106 

93.9 

84.8 

75 

1104 

488 

313 

231 

184 

154 

133 

118 

106 

90 

1388 

599 

381 

280 

222 

185 

l6o 

1 41 

127 

10? 

1692 

716 

451 

330 

261 

217 

187 

165 

146 

120 

2013 

838 

523 

381 

301 

250 

214 

189 

170 

180 

3557 

1377 

833 

595 

464 

382 

325 

285 

254 

240 

4526 

1987 

1175 

826 

635 

51  d 

438 

381 

338 

300 

5365 

2640 

1548 

1071 

8l6 

659 

553 

478 

422 

Table  5 

Angle  T  (degree! )  between  two  generators  of  the  reflection  cone  as  a 
function  of  P  (degrees)  and  7  (degrees).  Prom  equation  (56). 


X 

10® 

20® 

a 

a 

CD 

O 

0 

VO 

0 

0 

5 

1 .74 

3.42 

5.00 

6.42 

7.66 

8.66 

9.40 

9.85 

10.0 

10 

3.46 

6.81 

9.96 

12.8 

15.3 

17.3 

18.8 

19.7 

20.0 

15 

5.15 

10.2 

14, 9 

19.2 

22.9 

25.9 

28.2 

29.5 

30.0 

20 

6.8i 

13.4 

19.7 

25.4 

30. 4, 

34.5 

37.5 

39.4 

40.0 

25 

8.42 

l6.6 

24.4 

31.5 

37.8 

42.9 

46.8 

49.2 

50.0 

30 

9.96 

19.7 

29.0 

37.5 

45.0 

51.3 

56.1 

59.0 

60.0 

35 

11.4 

22.6 

33.3 

43.3 

52.1 

59.6 

65 .2 

68.8 

70.0 

40 

12.8 

25.4 

37.5 

48.8 

59.0 

67.7 

74.3 

78.6 

80.0 

45 

14.1 

28.0 

41 .4 

65.6 

75.5 

83.3 

88.3 

90.0 

50 

15.3 

30.4 

45.0 

11™ 

71.9 

83.1 

92.1 

98.0 

100 

55 

16.4 

32.5 

48,4 

63.5 

77.7 

90.4 

101 

108 

no 

60 

17.3 

34.5 

51.3 

67.7 

83.1 

97.2 

109 

117 

120 

65 

18.1 

36.1 

53.9 

71.3 

87.9 

103 

117 

126 

130 

70 

18.8 

37.5 

56.1 

74.3 

92.1 

109 

124 

136 

140 

75 

19.3 

38.6 

57.8 

76.8 

95.5 

BBS 

130 

144 

150 

80 

19.7 

39-4 

59.0 

78.6 

98.0 

mi 

136 

152 

160 

85 

39.8 

59.8 

79.6 

99. 5 

139 

158 

170 

90 

40.0 

_ 

60.0 

80.0 

100 

HI 

140 

l6o 

18c 

Table  6 


Angle  between  the  normals  (?)  from  combinations  of  observations  from 
observing  stations  indicated  in  Figs. 13  to  19. 


Figure  No. 

Example  No. 

Station  1 

Station  2 

1 

g 

13 

H7 

1 

a 

b 

33 

14 

2 

a 

b 

11 

b 

c 

37 

c 

d 

30 

a 

d 

78 

e 

g 

76 

h 

k 

74 

f 

h 

22 

i 

d 

59 

a 

k 

77 

15 

125 

a 

b 

22 

c 

d 

25 

a 

d 

-- 

55 

1 6 

122 

a 

b 

27 

c 

d 

28 

b 

« 

6i 

18 

7 

a 

b 

24 

b 

c 

33 

a 

c 

57 

19 

112 

a 

b 

39 

b  i 

c 

37 

i 

a 

c 

76 

t 

1 

! 
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ABSTRACT 


OPTICAL  STUDIES  OF  SPACECRAFT 
FROM  THE  C  LOU  DC  R  OFT  FACILITY 

E.  T.  TYSON 

Date,  consisting  primarily  of  resolved  and  quasi-resolved 
images,  are  collected  with  a  48-inch  telescope  whose  effective 
tocal  length  can  be  varied  between  312  and  6000  inches.  Although 
image  recording  has  been  tried  using  vicon,  image  orthicon,  and 
direct  photography,  none  of  these  methods  has  been  found  to  match 
the  eye  of  a  trained  observer  for  signal  separation  from  noise  and 
adaptability  to  rapid  and  extreme  brightness  fluctuations. 

A  16-mm  movie  made  from  a  video  tape  recording  of  a  Saturn 
S-IVB  stage  was  shown,  which  provided  an  example  of  resolution 
quality,  and  the  effect  of  large  fluctuations  in  brightness. 
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University  of  London 
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ABSTRACT 


THE  DESIGN  AND  OPERATION  OF  SOME  DEVICES 
FOR  PHOTO-ELECTRIC  AND  PHOTOGRAPHIC 
REGISTRATION  OF  ARTIFICIAL  SATELLITES 

Kenneth  Fea 

Photoelectric  techniques  are  applied  to  satellite  position 
measurement,  where  precision  is  comparable  to  that  of  the 
highest  quality  Baker-Ninn  photographic  techniques.  The 
method  described  employs  a  6 -inch  f/  8  objective  with  a 
10-inch  tri-alkali  photocathode  behind  an  M-shaped  linear 
slit  system.  Satellite  images  trail  across  the  slits  producing 
pulses  in  the  photoelectric  output.  Intervals  between  pulses 
are  related  to  image  position  through  the  M  pattern. 
Simultaneous  photography  and  photoelectric  sensing  of  stars 
is  used  for  equipment  calibration.  Resulting  satellite  positions 
are  measurable  to  one  arc  second. 
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RESULTS  OF  PHOTOMETRIC  AND  POLARIMETRIC 
MEASUREMENTS  OF  VARIOUS  SATELLITES 
USING  NASA'S  MOBILE  PHOTOMETRIC  OBSERVATORY 


Richard  H.  Emmons  and  Raymon  J.  Preski 
Goodyear  Aerospace  Corporation 
Akron,  Ohio 
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RESULTS  OF  PHOTOMETRIC  AND  PCIARIMETRIC  MEASUREMENTS  OF 


VARIOUS  SATELLITES  USING  NASA'S  MOBILE  PHOTOMETRIC  OBSERVATORY* 

Richard  H.  Emmons+  and  Raymond  J.  Preski+ 

Goodyear  Aerospace  Coiporation 
Akron,  Ohio 


1.  INTRODUCTION 

Under  contracts  NASl-6189  and  NAS1-6436  with  the  National  Aeronautics  and 
Space  Administration  (NASA)  Langley  Research  Center,  Goodyear  Aerospace 
Corporation  has  obtained  and  analyzed  multicolor  photometric  measurements  of 
the  Echo  I,  Echo  II  and  PAGEOS  I  balloon-type  satellites,  to  establish  a 
time-history  of  their  approximate  size,  figure  and  their  surface  reflectance, 
specularity  and  polarization  as  functions  of  wavelength.  These  observations 
were  made  using  NASA's  Mobile  Satellite  Photometric  Observatory  which  was 
developed  for  the  Langley  Research  Center  by  Goodyear  Aerospace  Corporation 
under  Contract  NAS  I-558O. 


2.  INSTRUMENTATION 

The  mobile  observatory's  new,  cryogenically -cooled  UBVRI  photometer  head, 
shown  in  Figure  1,  provides  up  to  four  simultaneous  intensity  signal  channels 
with  two  pairs  of  matched  S-20  and  S-l  photomultipliers ,  arranged  to  beam-share 
and  simultaneously  photometer  both  the  target  field  and  the  adjacent  sky. 

For  the  photometric  mode,  matched  pairs  of  field  stops  from  30  arc  seconds 
to  4  arc  minutes  are  provided.  UBV  color  filters  are  time-shared  on  the  S-20 
detectors,  while  VR1  filters  are  time-shared  on  the  S-l  detectors.  Polarimetric 
measurements  may  alternatively  be  made  in  each  of  these  colors  by  re-positioning 
both  the  single  2- minute  arc  central  field  stop  and  the  stowed  Wollaston  prism 
into  the  light  path,  and  then  manually  rotating  the  entire  photometer  head  on  a 


This  work  is  supported  by  the  Langley  Research  Center  (NASA),  Hampton,  Va.  , 
under  Contract  NAS1-6436. 

+Space  Systems  and  Analytics  Division,  Goodyear  A<  rospace  Corporation, 

Akron,  Ohio. 


precision  bearing  through  seven  keyed  36-degree  steps.  UBV  polarimetric 
measurements  can  also  be  made  automatically  by  repositioning  the  stowed 
achromatic  half-wave  plate  into  the  optical  path  and  initiating  its  automatic 
36-degree  steppec  rotations.  Signals  indicating  the  intensities,  amplifier 
gains,  color  filters,  and  half-wave  plate  orientations  are  fed,  together 
with  radio  time  signals,  to  a  recorder  in  the  instrumentation  room. 

Figure  2  presents  the  normalized  response  curves  and  effective  wave¬ 
lengths  for  the  standard  UBVRI  color  bands.  These  five  bands  span  about 
70  percent  of  the  solar  energy.  The  availability  of  hundreds  of  photometric 
reference  stars  which  have  been  repeatedly  measured  in  these  colors  permits 
highly  accurate  calibration  of  the  UBVRI  satellite  'observations . 


The  principle  by  which  sunlit  spherical  satellite's  specularity  may  be 
determined  from  the  photometric  observations  is  presented  in  the  following 
equations . 

Illuminance  ratios  of  spherical  satellites: 


F  E 

Specular:  sp  =  Antilog  (-O.^m)  =  l/b  _s 

E  E  ’ 

o  o 

Diffuse:  ji  =  Antilog  ( -0.4m)  =  2/3  fs 


R 

IT 


M  f  ( 

D 


f  (  0  )  =  i  £sin  0  +  (jt  •  0  )  cos  0  j 


Regression  equation  for  determining  specularity: 

Antilog  (-0.4m)  =  1/4a  +  2/?B  f(  0  ) 

A 


Specularity 


A  +  B 


The  first  two  equations  (from  Tousey,  Reference  l)  predict  the  extra- 
atmospheric  stellar  magnitude,  m  (or  illuminance  ratios,  —  )  of  100  percent 
reflecting  specular  and  diffuse  spheres  of  radius  R  and  distance  D.  The 
phase  function,  f  (  0  )  (from  Russell,  Reference  2)  gives  the  dependence  of 
the  illuminance  upon  phase  angle  for  a  perfectly  diffuse  sphere,  whose  surface 
reflects  in  accordance  with  Lambert's  cosine  law.  By  eliminating  the  common 
factors  in  these  equations,  the  regression  equation  is  found  which  weighs  and 
sums  the  specular  and  diffuse  components  in  relation  to  the  processed  satellite 
observations  (Reference  3).  Prior  to  performing  this  regression  analysis,  the 
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satellite's  reduced  magnitudes  are  corrected  for  Earth  albedo  and  normalized 

to  a  uniform  distance  and  the  instantaneous  phase  angles  are  computed.  After 

the  best-fit  solution  is  obtained  for  the  two  weighting  coefficients  (specular 

(A)  and  diffise  (B))  the  "specularity",  as  a  fraction  of  the  total  reflected 

A 

light,  is  found  simply  by  eva.uating  - - —  . 

For  spherical  satellites  that  sre  found  to  be  highly  specular,  the  small 
diffuse  component  of  the  reflected  sunlight  can  be  computed  and  subtracted 
from  the  illuminance.  This  net  "specular  magnitude"  can  then  be  employed  in 
the  theory  to  determine  either  the  satellite's  reflectance  or  radius  of 
'curvaTure*  of 'parametrically,  both. 

3.  OBSERVATIONS  AND  RESULTS 

The  first  observations  with  NASA's  Mobile  Satellite  Photometric  Observa¬ 
tory  were  made  during  the  summer  of  19 66  from  Palomar  Mountain,  California, 
from  a  peak  opposite  that  occupied  by  the  famous  200-mch  Hale  telescope. 

The  mobile  observatory's  original  single  channel  UBV  photometer  was  used 
there  to  observe  both  the  Echo  I  and  PAGEOS  satellites,  the  latter  during 
the  first  60  days  of  its  orbital  life  (Reference  4).  Upon  completion  of  the 
new  UBVRI  photometer  and  polarimeter  these  satellites  were  re-observed,  as 
well  as  Echo  II,  from  the  desert  floor  near  Yuma,  Arizona,  in  the  winter  of 
1967.  Table  I  presents  the  results  of  all  these  photometric  observations, 
together  with  seme  early  (1964)  visual-band  results  for  Echo  I. 

The  principal  findings  from  the  visual-band  photometric  data  for  the 
three  satellites  are  summarized  in  Table  TI. 

The  measured  color  indices,  in  stellar  magnitudes,  are  presented  in 
Table  III  together  with  the  reference  solar  values.  The  ultraviolet  deficiency 
in  the  light  reflected  from  Echo  I  and  PAGEOS  is  as  expected,  and  appears  to 
increase  with  time. 

Figure  3  presents  PAGEOS'  indicated  mean  and  extreme  radii  of  curvature 
versus  1 966  date,  from  the  Palomar  observations.  The  effect  of  the  first 
contact  with  the  Earth's  shadow  on  July  l4  is  clearly  shown.  The  relative 
stability  regained  in  early  August,  after  six  weeks  in  orbit,  may  be  associated 
with  the  loss  of  most  or  all  of  the  more  volatile  inflation  gases. 


no 


Figure  U  presents  the  visual  band  polarimetric  determinations  versus 
phase  angle.  The  solid  curve  is  the  theoretical  curve  for  bare  aluminum. 
There  seems  to  be  an  indication  here  that  Echo  II  may  have  lost  some  of  its 
amorphous  phosphate  thermal  balance  coating. 
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Figure  2  -  Normalized  Response  Curves  and  Effective 

Wavelengths  for  the  standard  UBVRI  Color  Bands 
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SOLAR  ENERGY,  PERCENT 


Figure  3  - 


LGECS 1  INDICATED  MEAN  AND  EXTREME  RADII  OF 
TRVATURE  VERSUS  1966  DATE,  FROM  PALOMAR 


OBSERVATIONS 


134 


V  -  BANO  POAMMETRC  OBSERVATIONS 


PHASE  ANCLE  (DECREES) 


Figure  4  -  Visual  Band  Polarimetric  Determinations 
versus  Phase  Angle 
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TABLE  I.  ANALYSIS  OF  PHOTOMETRIC  OBSERVATIONS 
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Based  on  laboratory -measured  reflectances  of  surface  materials. 
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AN  APPROXIMATE 


’.lETHOD  OF  ESTIMATING-  THE  REFLECTION 


CHARACTERISTICS  OP  A  SPHERICAL  SATELLITE 

by 

D.  E.  Smith 


Summary 

The  sunlight  reflected  from  a  spherical  satellite  is  assumed  to  be  composed 
of  an  amount  that  is  specularly  reflected  and  an  amount  that  is  diffusely  reflected. 
It  is  shown  that  the  amount  of  reflected  sunlight  falling  on  a  unit  area  a  great 
distance  from  the  satellite  is  directly  related  to  the  amount  that  is  specularly 
reflected  and  by  inference  to  the  amount  that  is  diffusely  reflected.  The 
apparent  stellar  magnitude  of  the  satellite  is  used  to  estimate  the  fraction  of 
solar  radiation  scattered  by  the  satellite  that  is  specularly  reflected.  The 
method  has  been  applied  to  the  Echo  2  satellite. 
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Introduction 


The  orbits  of  large  balloon  satellites,  such  as  Echo  1,  Echo  2  and  Pageos, 
are  subject  to  perturbations  hy  the  forces  of  3olar  radiation  pressure  acting 

o 

on  the  satellite.  These  perturbations  are  proportional  to  the  area-to-mass 
ratio  of  the  satellite  and  to  a  dimensionless  quantity,  often  termed  the  "reflection 
coefficient".  This  coefficient  is  a  measure  of  the  energy  transferred  from  the 
photons  to  the  satellite  and  depends  on  whether  the  incident  radiation  is  absorbed 
or  reflected;  and  if  refl'  cted,  whether  the  reflection  is  specular  or  diffuse^ 

For  example,  if  the  satellite  i3  a  perfect  sphere  and  specularly  reflecting,  so 
that  the  incident  radiation  falling  on  the  satellite  is  reflected  with  equal 
intensity  in  all  directions,  this  quantity  is  unity.  If,  however,  all  the 
incident  radiation  is  uniformly  diffusely  reflected  from  the  satellite  the 
"reflection  coefficient"  is  about  1.44  • 

For  the  accurate  computation  of  solar  radiation  pressure  perturbations  of 
a  satellite  orbit  it  is  obvious  that  a  knowledge  of  the  way  sunlight  is  reflected 
from  the  satellite  is  essential.  In  addition,  a  knowledge  of  the  way  light  ia 
reflected  can  be  an  indication  of  the  condition  of  the  surface  of  the  satellite. 

In  the  present  3tudy  the  amount  of  sunlight  reflected  from  a  spherical 
satellite  of  known  radius  and  falling  on  a  unit  area  at  a  large  but  known  distance 
from  the  satellite  is  investigated  as  a  function  of  the  fraction  of  the  incident 
sunlight  that  is  reflected  specularly.  It  is  assumed  that  a  fraction  a  of  the 
incident  radiation  is  reflected  and  that  a  is  composed  of  a.  the  specular  albedo 
and  a ^  the  diffuse  albedo.  It  is  assumed  that  the  fraction  1  -  a  ia  absorbed. 

2.  Specular  Rejection 

Let  I  be  the  intensity  of  solar  illumination  of  a  unit  area  normal  to  the 
direction  of  the  sun  and  F  the  total  radiation  falling  on  a  spherical  satellite 
of  radius  a,  then 

F  =  rra2  I  (1) 

If  ine  satellite  is  a  perfect  specular  reflector,  all  the  incident  radiation  is 

reflected  uniformly  over  the  entire  unit  sphere.  The  flux  of  radiation  reflected 

cpcculnr_y  in  a  unit  solid  angle  is  therefore 

_F 

Uir 
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and  the  intensity  of  illumination,  E  ,  of  a  unit  area  distance  r  from  the  satellite 
(where  r  is  considerably  greater  than  a,  the  radius  of  the  satellite)  is  given  by 


If  only  a  fraction  a of  the  incident  radiation  is  specularly  reflected. 


(2) 

then 


3. 


asF 

E  s  — 2— , 

s  Zftfr 


^  (&  Y 


(3) 


Diffuse  keflection 


A  spherical  coordinate  system  centred  on  the  satellite  can  be  constructed  such 
that  the  fundamental  reference  plane  contains  the  sun  and  the  eye  of  the  observer, 
and  the  normal  to  the  reference  plane  is  the  polar  axis  of  the  spherical  system. 

(see  Pig.  l). 

Let  the  angle  between  the  sun  and  the  observer,  distant  r  from  the  centre  of  the 
satellite,  as  measured  in  the  reference  plane  be  0  and  any  point  on  the  satellite 
be  designated  by  A,  the  longitude  measured  in  the  reference  plane  from  the  satellite- 
sun  line,  and  9,  the  latitude  of  the  point  with  respect  to  the  reference  plane.  An 
elementary  area,  ds,  is  therefore  given  by 

ds  =  a3  cos  <p  d  <p  dA  (4) 

The  solar  radiation  (d?)  falling  on  ds  is 

dP  =  la3  cos  9  cos  /?  d  9  dA  (3) 

where  /3  is  the  angle  between  the  normal  to  ds  and  the  sun  and  from  spherical 
trigonometrical  considerations  in  Pig.  1,  we  have 

cos  /3  =  cos  9  cos  A  (6) 


and  therefore 

dF  =  la3  cos3  9  cos  A  d9  dA  (7) 

If  the  radiation  falling  on  the  elementary  area  ds  is  diffusely  reflected 

according  to  a  general  function  f(y)  where  y  1*  the  angle  between  the  normal  to 

ds  and  the  observer  (see  Pig.  l),  then  the  intensity  of  illumination,  dE^,  of  a 

unit  area  at  the  observer,  distance  r  from  the  satellite  (r  considerably  theater 
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tr.an  a )  ,  13 


(8) 


=  T7  003  Y  f(  y)  dJ 


where  is  the  diffuse  albedo.  The  total  diffusely  reflected  solar  radiation, 
~D*  fu-l--ne  on  a  unit  area  distance  r  from  the  satellite  is,  therefore,  given  by 
the  integral  of  dE^  over  that  part  of  the  satellite  that  is  illuminated  and  can  be 
seen  from  the  unit  area*  Thus 


ir 

r*  =  2 


*  -j 


%  -  — 

-2 

!  d^ 


.  1  { 
2  A 


(9) 


a  n 
=  9  "  2 


=  °D* 


TT 

9  =  2 


X  -  - 

,  2 


J  IT  1 

9  *  -  o  A, 


cosa  9  cos  X  cos  y  f(r)  d<P  dX 

-•-I 


(10) 


From  Fig.  1  we  can  obtain 

cos  y  =  cos  <p  cos  (9  -  X) 

and  hence 


(11) 


,X  “  2 


COS 
as  0 


9  cos  X  cos  (9  -X)f(Y)  &9  dX 

ir 

'  2 


(a)  perfect  diffuser 

If  the  surface  of  the  satellite  is  a  uniformly  diffuse  reflector,  the 
luminance  is  the  same  in  all  directions  and  f(Y)  takes  its  simplest  form,  that 
is,  a  constant*  We  may  then  write 


f(r)  =  *, 

where  is  unknown  and  needs  to  be  evaluated* 

The  total  illumination  on  the  surface  of  a  hemisphere  of 
radius  r  with  ds  at  the  centre  is  equel  to  the  total  radiation  diffusely 
reflected  (ot^dF)  from  ds.  We  therefore  have 


(12) 


(13) 


ot^dF  =  dEp  r3  sin  y  dY  d t/i 


(14) 


hemisphere 


143 


where  ip  As  the  longitude  of  a  point  on  the  surface  of  the  hemisphere  as  shown  in 
Fig.  2.  p'rom  equations  8  and  12*.  v;e  obtain 


ip  =  2n 

“D^  =  J 

ip  -  0 


Y  =  2 

an  sin  y  cos  Y  ^(y)  dF  dy  dip 

jD 

Y  =  0 


and  hence 


2ir 


r  =  2 

j  sin  y  Y  *(y)  dy  »  1 
Y  =  0 


(15) 


when  f ( y)  =  ,  we  ob  tain  from  equation  1 5 


k 


1 


1_ 


(16) 


Performing  the  integration  in  equation  12  we  obtain 


*D 


21a3 

Jwr* 


°D 


(w  -  d)  cos  9  ♦  sin  9 


(17) 


The  above  equation  gives  the  total  diffusely  reflected  radiation  falling  on  a 
unit  area  distance  r  from  the  satellite  and  angle  9  from  the  sun  when  the  surface 
of  the  satellite  is  a  uniform  diffuser* 

(b)  imperfect  diffuser 


144 


■'r:  general,  f(y)  is  not  constant  and  the  luminance  is  therefore  not  the 
same  in  all  directions.  Let  us  consider  the  case  of 


f(y)  =  k2  cos  y 

From  eru.ition  15*  we  obtain 


and  fro..:  equations  11  and  12 


(18) 


(19) 


,  3Iaa 
“  2u?r 


*  2 

ap  cos4  <|<d<f> 
J  ir 

<p  =  -  *2 


.  n 

Ks2 


cos  A  cosa  (0  -  A)  dA 


(20) 


A  =  d  -  ^ 


which  leads  to 


ED  =  itF"  aD  ^  *  coa  (21 ) 

This  ecuation  gives  the  total  diffusely  reflected  radiation  falling  on  a  unit 
area  when  the  reflection  follows  a  cos  flaw. 


4»  Apparent  Brightness  of  .Satellite 

In  jenerul,  the  radiation  reflected  from  a  satellite  will  be  composed  of 
radiation  that  is  specularly  reflected  and  radiation  that  is  diffusely  reflected. 
Thus  if  is  the  total  reflected  radiation  falling  on  a  unit  area  at  a  distance 
r  from  the  satellite  we  have 


E 


o 


E 

3 


♦  £4 


(22) 


Further,  if  we  assume  that  a  fraction  a  of  the  incident  radiation  is  reflected, 
specularly  .and  diffusely,  and  8  is  the  fraction  of  a  that  is  specularly  reflected, 

then 


a  s  aS 
s 

(23) 

a-  =  u  (1  -8) 

U 

(24) 
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If  the  type  of  diffuse  reflection  is  known  and  values  of  Eq  are  obtained 

fron  say,  a  photo-electric  photometer  measurement,  we  can  derive  a  value  for  8. 

For  the  Echo  2  balloon  satellite  (1964-O^A),  which  we  shall  assume  to  be  perfectly 

(2) 


spherical,  we  have  a  =  20.5  metres 
from  equation  3 


When  r  =  1000  km,  we  therefore  have 


•sr  =  1.051  x  10"10  x  a 

1  9 


anti  from  equation  1  7 


(25) 


=  0.692  x  10"10  x  Ojj  (w  -  d)  cos  9  +  sin 

when  the  satellite  is  a  perfect  diffuser,  that  is  f(y)  =  ^  * 
and  from  equation  21 


(26) 


0.768  x  10-10  x 


Ojj  1  ♦  COS 


when  f(y)  =  ^  cos  Thus,  from  equations  Z  to  27,  we  have 


(27) 


x  1010  =  1.051  aB  +  0.892  a  (1  -  B)  (jt  -  9 )  cos  9  ♦  sin 

when  f(y)  =  “  ,  and 

x  1010  =  1.051  aB  ♦  0.788  a  (1  -  B)  £  1  ♦  cos  9  J  * 
when  f( y)  cos  y» 

The  variation  of  fo  with  9  for  8=0,  0.5  and  1  is  shown  in  Fig.  3.  For 

“T  1 

values  cf  0,  phase  angle,  greater  than  2  there  is  only  a  small  difference 

between  the  two  sets  of  curves.  For  small  values  of  9  the  difference  is  larger 

and  for  constant  9  increases  as  B  decreases.  Tt  is  interesting  to  note  that 

E 

the  graphs  show  that  for  0=0,'—  is  about  3  times  larger  for  completely 

diffuse  reflection  (B  =  0)  than  for  a  completely  specular  reflection  (B  =  l). 

l4> 


(28) 


(29) 


tn  *  ay  30th,  1966  the  Echo  2  satellite  was  tracked  by  the  U.S.A.F.  at  the 
fi'right-Patterson  Air  Force  Base,  Ohio  using  a  photo-electric  photometer^ 

The  measurements  obtained  at  Ohio  on  this  pass  have  been  compared  with  the  theory 
outlined  in  this  and  previous  sections  of  this  paper  and  used  to  estimate  a  value 
of  the  fraction  of  scattered  radiation  that  was  specularly  reflected. 

The  transit  of  Echo  2  at  about  0420  U.T.  on  Hay  30th,  1966  passed  almost 
overhead  of  Ohio,  rising  in  the  north  and  entering  eclipse  in  the  south.  Measure¬ 
ments  of  the  apparent  stellar  magnitude  were  made  continuously  for  over  five  minutes. 
The  magnitude  changed  steadily  throughout  this  time  and  showed  no  violent 
oscillation  or  jumps  and  the  transit  appeared  to  be  a  very  suitable  one  for 
comparing  with  the  analysis  already  described. 

Six  values  for  the  apparent  stellar  magnitude,  at  minute  intervals,  were 
taken  from  the  photometer  measurements  and  corrected  for  almospheric  extinction 

and  range.  The  extinction  correction.  Am  ,  was  calculated  on  the  basis  of  0.5 

© 

stellar  magnitude  par  air  mass  (measured  at  Ohio  at  the  time  of  the  observations) 
from  the  equation 

Extinction  correction.  Am  »  0.5/  sin  E  magnitudes  (30) 

where  E  is  the  elevation  of  the  satellite.  The  range  correction  was  calculated 
for  a  standard  range  of  1000  las  on  the  basis  that  the  brightness  of  the  satellite 
varied  inversely  as  the  square  of  the  range.  The  apparent  magnitude  (m)  of  an 
object  is  defined  in  terms  of  the  brightness  (B),  as  follows^^ 

m  =  -  2.5  1o&j0(B)  (31) 

Thus,  two  brightness  values  (suffixes  0,  l)  are  related  in  apparent  magnitude  by 

£ 

ni  -  no e  2-5  lo«io  (b;)  (3z) 

and  since 

!»  -  ft i\‘ 

'  VJ 
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when.  p^  y  p^  arc  ‘he  ranges  ol’  the  objoct,  we  have 


°1  -  »0  -  5.0  lo£10 

If  in1  i-5  the  stellar  magnitude  at  a  range  of  1000  ha  (*  )  and  m0  is  ths  stellar 

lanpr.i ‘ui<-  at  the  observed  range  (pQ)»  then 

range  correction,  Amr  =  5.0  ■*,o®10 

and  ‘he  corrected  apparent  stellar  magnitude  (m  )  is  therefore  given  by 

c 

m  =  m  —  Am  —  Am 
co  e  r 

Detailc  of  the  satellite  pass,  the  observod  magnitudes  of  the  satellite,  the 
corrections  applied  and  the  values  of  mc  are  shown  in  Table  1, 

TA3L1;  1«  Details  of  hcho  2  observed  from  Ohio  on  May  30th.  1966 


1 - 

Time 

m 

E 

>o 

Phase 

Ain 

Am 

m 

!  (U.T.) 

1 

0 

(degs) 

(km) 

9  (degs) 

e 

.  _ 

r 

c 

0418 

♦ 

0.82 

31.6 

1825 

112 

i 

♦  0.95 

♦  1.31 

-  1.44 

04 19 

♦ 

o 

• 

o 

42.5 

1533 

104 

♦  0.74 

♦  0.93 

-  1.60 

0420 

- 

0.63 

57.9 

1304 

94 

♦  0.59 

+  0.58 

o 

CO 

. 

i 

0421 

- 

1.12 

77.0  , 

1162 

80 

+  0.51 

+  0.33 

-  1.96 

04 22 

- 

1.24 

81.1 

1136 

64 

+  0.51 

+  0.28 

-  2.03 

0425 

- 

1.10 

_ 

61 .6 

1 

_ 

1277  j 

1 

53 

+  0.58 

I 

♦  0.53 

. 

-  2.21 

The  sources  of  error  in  the  values  of  m  arise  from  errors  in  the  observed 

c 

majii;tu*l«  (m  ),  the  extinction  correction  (Am  )  and  the  range  correction  (Am  ). 

o  ox* 

The  noise  in  the  photo-electric  photometer  recording  had  an  overall  width 
of  about  0.2  magnitudes  and  was  approximately  constant  throughout  the  record. 
The  central  value  taken  from  the  record  had  a  standard  deviation  of  about  0.05 


The  error  in  the  extinction  correction  is  primarily  due  to  any  error  in  the 
measured  value  of  0.5  magnitudes  per  air  mass.  The  accuracy  of  this  figure  is 
probably  of  the  same  order  as  the  accuracy  of  the  satellite  magnitude,  that  is, 

O.05  Magnitudes  but  for  the  extinction  correction  the  error  is  systematic.  Assuming 
a  figure  of  0.05  magnitudes  the  systematic  error  of  the  extinction  correction  lies 
between  C.05  and  0.1  magnitudes,  depending  on  the  elevation.  Additional  errors 
in  the  extinction  correction  arising  from  errors  in  the  computed  elevation  of  the 
satellite,  whichAelieved  to  be  accurate  to  about  0.2  degrees,  can  be  shown 
to  be  negligible. 

The  accuracy  of  the  computed  range  of  the  satellite  is  thought  to  be  better 
than  10  km  and  is  a  systematic  error  rather  than  a  random  error  because  the 
relative  accuracy  of  one  predicted  range  with  another  is  probably  of  the  order  of 
1  km.  A  systematic  error  of  10  km  introduces  a  systematic  error  in  the  range 
correction  of  0.03  to  0.04  magnitudes,  depending  on  the  range.  The  error  in  mc, 
the  corrected  stellar  magnitude,  is  therefore  composed  of  a  standard  deviation 
of  about  0.05  magnitudes  and  a  possible  systematic  error  between  0.09  magnitudes 
and  0.12  magnitudes,  depending  on  the  elevation  and  range  of  the  satellite. 

The  effect  of  a  systematic  error  on  the  corrected  stellar  magnitude  is  a 
reduction  or  an  increase  of  all  the  values  by  approximately  the  same  amount. 

Thus,  a  systematic  error  will  change  the  individual  values  of  mc  but  the  shape 
of  the  magnitude  versus  phase  curve  will  be  largely  unaffected. 

The  theoretical  apparent  stellar  magnitude  (m^)  can  be  obtained  from  the 
theoretical  brightness  of  the  satellite  in  the  following  way.  Prom  equation  32, 
we  obtain 

mt  =  M  +  2.5  log1£) 

where  K  is  the  apparent  stellar  magnitude  of  the  sun  (equal  to-26.78^^ ) ,  EQ  is 
the  intensity  of  illumination  of  a  unit  area,  distant  r,  by  sunlight  reflected  from 

f’e  sa'-  Llite  and  I  is  the  intensity  of  illumination  of  direct  sunlight. 
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The  values  of  m  are  shown  in  Pigures  4  And  5  superimposed  on  curves  of  m. 

c  t 

for  8*0,  0.5,  0.6  and  1.0  when  a  ■  1*  The  standard  deviation  of  each  value 
of  o>c  is  indicated  by  the  length  of  the  vertical  lines  above  and  below  the  point* 

■Hie  systematic  error  is  not  shown.  Curves  >f  m^  when  a  F  1  are  identical  in  shape 
to  those  for  a  *  1  but  are  displaced  so  as  to  decrease  the  theoretical  magnitude. 

A  comparison  of  the  slope  of  the  observed  points  with  the  theoretical  curves  can 
provide  an  estimate  of  8,  the  fraction  of  scattered  sunlight  that  is  specularly 
reflected. 

An  inspeotion  of  Figure  4,  in  which  the  theoretical  curve  is  based  on  f(y)  ■ 
suggests  that  •  lies  in  the  range  0.5  to  0.6  and  Figure  5,  based  on  f(y)  ■  3oos  y/2t, 
suggests  that  fl  probably  lies  in  the  range  0.55  to  0.65*  In  addition,  there  is 
some  indication  that  there  is  better  agreement  between  the  observed  and  theoretical 
values  in  Figure  4  th.vn  in  Figure  5,  suggesting  that  the  fraction  of  sunlight  that 
is  diffusely  reflected  is  being  reflected  nearly  uniformly  in  all  directions. 

It  is  important  to  remember  that  the  assumption  that  the  satellite  is  perfectly 
spherical  may  not  be  entirely  true  and,  if  not,  could  introduce  systematic  errors. 
Errors  of  several  percent  from  this  source  cannot  be  ruled  out. 

Neglecting  these  possible  errors,  the  present  study  suggests  that  tf,  the 
fraction  of  scattered  solar  radiation  that  is  refleoted  specularly  lies  in  the 
range  0.50  te  0.65,  say 

8  *  0.57  ±  0.07  (38) 

6.  Sumnary  and  Conclusions 

The  apparent  brightness  of  a  spherical  satellite  is  shown  to  be  s  measure 
of  the  fractional  part  of  the  scattered  solar  radiation  that  is  speoularly  reflected. 
Expressions  have  been  derived  for  the  amount  of  reflected  radiation  falling  on  a 

unit  area  for  specular  reflection  and  for  two  models  of  diffuse  reflection. 
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It  h.is  been  assumed  that,  in  general,  a  fraction  a  of  the  solar  radiation  incident 

on  the  satellite  is  reflected  and  is  composed  of  a  fraction  a  that  is  specularly 

s 

reflecting  and  a  fraction  that  is  diffusely  reflected. 

The  theoretioal  stellar  magnitude  of  the  Ecno  2  satellite  at  a  range  of  1000  km 
has  been  compared  with  the  observed  magnitude  (corrected  for  this  range)  for 
different  phase  angles  and  has  indicated  that  about  57  percent  of  the  solar  radiation 
scattered  by  Echo  2  is  specularly  reflected,  the  remaining  43  percent  being 
diffusely  reflected. 

Any  major  deterioration  in  the  surface  of  the  satellite  or  extensive  small 
scale  changes  in  the  shape,  such  as  wrinkling,  will  almost  certainly  affect  the 
reflection  characteristics  of  the  satellite  and  therefore  its  apparent  brightness  as 
a  function  of  range  and  phase.  The  present  paper  has  described  a  simple  method 
of  estimating  these  reflection  characteristics  from  ground  based  measurements  of 
the  apparent  stellar  magnitude  which  might  be  used  as  a  measure  of  the  deterioration 
(if  any)  of  the  surface  of  the  satellite. 
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COORDINATE  SYSTEM  FOR  DIFFUSE  REFLECTION 
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FIGURE  2 
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FIGURE  3 
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PHASE  0 

THEORETICAL  VARIATION  OF  REFLECTED 
RADIATION  FROM  ECHO  2  WITH  PHASE 
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FIGURE  5 
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DESIGN  OF  TELESCOPE  OPTICAL  SYSTEMS 
FOR 

VISUAL  OBSERVATIONS  OF  ORBITING  SATELLITES 


By 

Arthur  S.  Leonard* 

The  first  telescopes  used  for  visual  satellite  tracking  were  designed 
as  Richest-Field^ ^  instruments.  The  purpose  of  a  "Richest-Field"  telescope 
is  to  afford  spectacular  views  of  the  heavens,  and  to  show  the  greatest 
possible  number  of  Milky  Way  stars  at  one  view. 

The  intention  of  the  design  of  the  Richest-Field  telescope  is  to 
provide  an  optimum  combination  of  a  wide  field  of  view  and  reaching  to 
fairly  faint  stars.  This  is  just  about  what  is  needed  for  a  satellite 
tracking  instrument,  so  it  is  not  at  all  surprising  that  the  Richest-Field 
design  was  selected  for  the  first  satellite  telescopes. 

Standard  specifications  for  Richest-Field  design  are  an  objective  of 
f/4  to  f/6  focal  ratio  and  a  wide-angle  eyepiece  to  give  3.39  power  per 
inch  of  aperture  of  the  objective.  Arguments  put  forth  in  the  earliest 
published  articles  on  this  design  indicate  that  these  specifications  were 
derived  almost  entirely  from  theoretical  considerations  and  assumptions 
which  seemed  reasonable  at  the  time. 

Many  Richest-Field  telescopes  have  been  built,  and  their  builders 
have  attested  to  tne  spectacular  views  which  they  afford;  but  there  appears 
to  have  been  little  systematic  experimental  testinq  to  ascertain  whether 
or  not  these  specifications  lead  to  a  truly  optimum  design.  It  is  the 
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writer's  opinion  that  they  are  in  rather  serious  error  on  two  points.  First, 
the  optimum  specific  magnification  (magnifying  power  per  inch  of  aperture 
of  the  objective)  should  be  somewhere  between  5-  and  10-  power  per  inch, 
depending  on  the  individual's  feelings  for  what  he  sees  in  the  eyepiece  of 
his  telescope,  rather  than  at  the  earlier  published  value  of  3.39;  and 
second,  an  increase  in  the  focal  ratio  of  the  objective  system  from  the 
f/4  -  f/6  range  to  the  f/12  -  f / 1 7  range  will  give  a  significant  improvement 
in  the  sharpness  of  star  images  over  a  large  fraction  of  the  field  of  view. 

Limiting  Magnitude  versus  Specific  Magnification 

One  of  the  assumptions  made  in  setting  up  the  earlier  specifications 

was  that  the  limiting  magnitude  of  a  telescope  is  determined  by  aperture 

alone  and  is  entirely  independent  of  its  magnifying  power.  A  more  recent 

(?) 

experimental  investigation  into  this  matter'  '  has  shown  that  limiting 
magnitude  is  very  much  a  function  of  magnifying  power  as  well  as  light 
gathering  power.  For  example,  an  increase  in  magnification  from  3.39  to 
10  power  per  inch  will  increase  the  limiting  brightness  of  a  small  telescope 
by  1.9  stellar  magnitudes.  The  reason  for  this  is  that  the  ability  of  the 
eye  to  detect  a  faint  star  image  depends  not  only  on  the  amount  of  light  in 
that  image,  but  also  on  the  apparent  brightness  of  the  surrounding  background. 
The  brightness  of  the  star  image  is  dependent  almost  entirely  on  the  light 
gathering  power  of  the  objective  while  the  brightness  of  the  surrounding 
background  presented  to  the  eye  varies  inversely  as  the  square  of  the 
specific  magnification. 

Although  the  impression  which  a  given  view  will  make  on  a  person  will 
probably  be  different  for  different  observers,  it  is  the  writer's  feeling 


that  the  relatively  milky  background  presented  by  a  specific  magnification 
of  3.39  power  per  inch  generally  detracts  from  the  beauty  of  the  view. 
Eight-to-ten-power  per  inch  gives  a  much  darker  background  and  greater 
contrast  between  the  stars  and  background.  Even  at  11.3  power  per  inch, 
which  the  writer  is  using  for  satellite  tracking,  the  background  illumina- 
ation  offered  by  the  darkest  part  of  the  sky  on  the  darkest  night  is  still 
sufficient  to  make  cross-wires  of  0.006-inch  diameter  easily  visible  in  the 
eyepiece  without  any  supplemental  illumination. 

One  of  the  implied  assumptions  used  in  writing  the  specifications  for 
the  Richest-Field  telescope  is  that  it  will  be  used  mainly  on  clear  dark 
nights.  Satellite  tracking,  on  the  other  hand,  must  be  carried  out  on  as 
many  occasions  as  possible;  and,  as  any  experienced  satellite  observer  knows, 
this  means  working  in  full  moon  light,  strong  twilight,  and  with  light  fog 
or  clouds  in  the  sky.  Thus  the  specifications  which  yield  the  best  Richest- 
Field  instrument  are  not  necessarily  the  best  prescription  for  a  satellite 
tracking  telescope.  As  the  background  illumination  is  increased,  higher 
magnification  is  of  increasing  value  in  revealing  faint  objects.  Also, 
averted  vision  is  very  helpful  in  attracting  an  observer's  attention  to  a 
faint  moving  object,  but  this  works  only  if  the  apparent  background  illumina¬ 
tion  is  fairly  low.  It  won't  be  of  any  help  in  broad  daylight.  For  these 
reasons  a  satellite  telescope  should  be  designed  with  a  higher  magnifying  power 
than  a  Richest-Field. 

If  a  satellite  appears  well  above  threshold  brightness  in  the  field  of 
view,  high  angular  velocity  does  not  interfere  appreciably  with  its  being 
seen.  If  it  is  near  limiting  brightness,  however,  a  high  apparent  velocity 
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increases  the  probability  of  its  beinq  missed  by  the  observer.  For  this 
reason  as  magnifying  power  is  increased  a  value  is  soon  reached  above  which 
there  is  little  net  gain  to  be  realized  by  going  higher.  From  long  experience 
observing  faint  satellites  with  a  6-inch  refractor  and  other  smaller  telescopes, 
it  is  the  writer's  opinion  that  10  power  per  inch  is  just  about  optimum 
specific  magnification  for  telescopes  up  to  6-inches  aperture.  Experiments 
with  higher  magnification  have  indicated  that  there  is  little  to  be  gained  by 
going  above  about  15  power  per  inch.  On  the  other  hand,  if  the  satellite  is 
so  bright  that  high  magnification  is  unnecessary  to  make  it  easily  visible 
under  the  existing  sky  conditions  and  if  its  angular  velocity  is  fairly  high, 
a  lower  power  is  to  be  preferred.  If  the  angular  velocity  of  the  satellite 
is  so  high  that  it  will  cross  the  field  of  view  in  less  than  about  two 
seconds,  some  loss  in  accuracy  of  the  observations  may  result.  Linder  such 
circumstances  a  lower  magnification  should  be  used,  assuming, of  course, 
that  the  lower  magnification  will  not  render  the  satellite  invisible  to  the 
observer's  eyes.  On  the  other  hand,  if  the  satellite  is  very  high  and 
moving  very  slowly,  an  increase  in  magnification  can  easily  lead  to  a  more 
accurate  observation. 

Eyepiece  Aberrations  versus  Focal  Ratio  of  Objective 

The  specification  for  the  focal  ratio  of  the  objective  of  the  Richest- 
Field  telescope  was  probably  set  as  a  matter  of  convenience  rather  than  from 
any  consideration  of  what  effect  it  might  have  on  the  over-all  performance 
of  the  instrument.  To  obtain  a  specific  magnification  of  3.39  ;  ower  per 
inch  with  a  conventional  f/15  refractor  objective  would  require  an  eyepiece 
of  4.42-inch  focal  length.  For  the  amateur  telescope  makers  who  built  the 
first  Richest-Field  telescopes,  such  an  eyepiece  was  --  and  still  is  -- 
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difficult  to  obtain;  and,  if  available,  would  be  a  relatively  costly  item. 
So,  they  followed  the  route  of  least  resistance  and  made  or  purchased 
objectives  to  match  their  eyepieces.  An  eyepiece  of  15-inch  focal  length  -- 
which  is  about  as  large  as  can  be  readily  obtained  --  requires  an  f/5.085 
objective  to  give  a  specific  magnification  of  3.39  power  per  inch. 

In  the  center  of  field  of  view,  some  of  our  best  eyepiece  designs  can 
handle  an  f/5  cone  of  rays  without  serious  aberration;  but  as  we  go  away 
from  the  center,  even  the  best  designs  get  into  rather  serious  trouble.  The 
very  low  power  employed  in  the  Richest-Field  telescope  helps  to  reduce  the 
apparent  magnitude  of  these  aberrations,  but,  as  the  original  designer  says: 
"Marginal  definition  must  be  a  little  inferior  to  central  in  every  Richest- 
Field  telescope, 

This  loss  in  performance  toward  the  outer  part  of  the  field  of  view  can 
be  reduced  to  negligible  proportions  by  simply  increasing  the  focal  length 
of  both  the  objective  and  eyepiece.  Although  this  will  require  larger  and 
more  costly  eyepieces,  to  obtain  a  magnification  of  10  power  per  inch  with 
a  standard  f/15  objective  requires  an  eyepiece  of  only  1.5-inch  focal  length. 

The  following  example  is  presented  in  order  to  explain  in  detail  why 
this  improvement  in  optical  performance  comes  about  when  we  increase  the 
focal  lengths  of  the  objective  and  eyepiece.  Consider  two  telescopes,  1 
and  2,  which  are  of  identical  optical  design  throughout.  They  are  of  the 
same  aperture  and  magnification;  but  2  has  a  longer  focal  length  than  1, 
and  the  radii  of  curvature  of  all  its  surfaces  are  greater  than  those  in  1 
by  the  same  ratio.  Both  telescopes  are  pointed  in  the  same  direction,  but 
we  are  looking  at  a  star  some  distance  off  the  axes  of  the  telescopes.  At 


corresponding  points  in  each  telescope  the  cones  of  rays  from  the  star  will 
be  of  the  same  diameter;  and  the  angles  of  incidence  of  the  central  rays  at 
corresponding  optical  surfaces  will  be  the  same  in  both  instruments. 

We  will  assume  that  the  design  of  the  objective  is  such  that  virtually 
no  aberration  is  produced  by  either  objective  and  all  the  observed  errors 
are  produced  in  the  eyepieces.  Consider  the  various  aberrations  produced 
at  the  first  surface  of  the  eyepiece  in  each  instrument.  Spherical 
aberration  is  given  by  the  following  expression: 


where  w$  is  the  maximum  wave-front  aberration  in  the  same  units  as  the  other 
linear  dimensions,  C$  is  a  number  which  depends  on  the  geometry  of  the 
rays  and  the  index  of  refraction  of  the  glass,  d  is  the  diameter  of  the  cone  of 
rays,  and  R  is  the  radius  of  curvature  of  the  lens  surface. 

Since  C$  and  d  will  be  the  same  for  both  instruments,  the  ratio  of 
spherical  aberration  produced  by  this  first  surface  in  the  two  instruments 
will  be: 


w 
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This  indicates  that  the  ratio  of  the  spherical  aberration  produced  by 
the  two  eyepieces  will  be  inversely  proportional  to  the  cubes  of  the  radii 
of  curvature  of  the  lenses  and  this  will  be  inversely  proportional  to  the 
cubes  of  the  focal  lengths  of  their  objectives. 


First-order  coma 


will  be  given  by : 
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w  =  C  — ■*- 
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where  w  is  the  maximum  wave-front  aberration  due  to  coma,  and  C  is  a 
c  c 

constant  which  depends  on  the  geometry  of  the  rays  and  the  index  of 
refraction.  Therefore,  the  ratio  of  coma  produced  by  the  two  eyepieces  will 
be  inversely  proportional  to  the  squares  of  the  focal  lengths  of  the  objec¬ 
tives. 

First-order  astigmatism  is  given  by: 


-s  ‘  cs  r 


(4) 


This  means  that  the  ratio  of  astigmatism  produced  by  the  two  eyepieces 
will  be  inversely  proportional  to  the  first  power  of  the  ratio  of  the  focal 
lengths. 

Chromatic  aberration  varies  in  a  way  similar  to  astigmatism,  so  it 
too  will  be  inversely  proportional  to  the  focal  lengths  of  the  objectives. 

Thus  in  telescope  2,  with  its  longer  focal  length,  all  of  these  eye¬ 
piece  aberrations  will  be  reduced.  If  we  go  to  an  objective  of  long 
enough  effective  focal  length,  all  of  these  aberrations  can  be  reduced  to 
negligible  amounts.  The  only  eyepiece  aberration  not  affected  by  this 
increase  in  focal  length  will  be  field  curvature.  In  a  well-designed  eye¬ 
piece  this  defect  is  not  serious,  so  we  can  practically  eliminate  all 
eyepiece  aberration  by  simply  going  to  an  objective  and  eyepiece  of  long 
enough  focal  length. 


Conclusions 


In  designing  a  telescope  for  satellite  tracking  the  focal  ratio  of  the 
objective  system  should  be  made  as  high  as  f/15  if  at  all  possible.  A 
refractor  is  to  be  preferred  over  a  reflector  because  with  the  former  it  is 
usually  possible  to  design  an  internal  system  of  light  baffles  which  will 
give  a  darker,  more  contrasty  field  than  with  a  reflector.  In  order  that 
the  instrument  should  be  able  to  operate  near  its  maximum  performance 
over  a  wide  range  of  observing  conditions,  and  satellite  brightness  and 
speed,  a  good  selection  of  eyepieces  should  be  provided  which  will  provide 
good  coverage  over  a  range  of  specific  magnification  from  about  5-  to  20- 
power  per  inch  of  aperture.  For  telescopes  of  more  than  6-inches  aperture 
the  range  of  specific  magnification  should  be  lowered  somewhat  from  these 
figures. 
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Abstract 

New  theoretical  results  are  derived  which  predict 
the  princi|>al  effect  of  the  gravitational  and  magnetic 
torqueiymjjjc  rotatajj^il  nyjtiqnjjf  an  orbiting  triaxial 
rigid  body.  Three  integrals  are  derived  by  the  averaging 
principle  of  Kryloff  !  igoliuboff  to  determine  the  motion 
of  the  rotational  angular  momentum  vector.  The  results 
are  compared  with  the  observations  of  the  rotational  mo 
tion  of  the  Pegasus  satellites. 


is  treated  here,  reduces  to  the  aforementioned  symmet¬ 
rical  problems  if  the  rigid  body  is  assumed  to  be  a  figure 
of  revolution  in  a  fixed  orbit  about  the  primary  mass. 

♦  <».  ♦  v  Theory 

Formulation  of_Di fferential  Equations  for  Rotational  Motion 
The  Newtonian  gravitational  torque  about  the  center 
of  mass  of  a  small  triaxial  rigid  body  orbiting  a  spherical 

(3) 

primary  mass  mav  be  written 


Introduction 

The  rotational  and  translational  motion  of  a  rigid 
body  are  generally  coupled.  Because  the  coupling  effect 
is  negligible  in  most  applications  in  celestial  mechanics, 
the  two  problems  may  be  treated  se|>arately.  They  are, 
however,  analogous  in  many  respects.  For  example,  if 
there  are  no  external  torques  acting  in  the  two -body  pro¬ 
blem  of  two-point  masses,  then  the  orbital  angular  mo¬ 
mentum  vector  is  constant,  and  the  motion  of  the  bodies 
is  called  "free  Keplerian.  "  Likewise,  if  there  are  no 
external  torques  acting  in  the  rotational  problem  of  a  rig¬ 
id  body,  then  the  rotational  angular  momentum  vector  is 
constant  and  the  motion  of  the  rigid  body  s  called  'tree 
Eulerian.  "  If  the  primary  body  in  the  translational  two- 
body  problem  i3  triaxial  or  oblate,  then  the  principal  ef 
feet  is  that  the  uim-  averaged  orbital  angular  momentum 
vector  follows  periodic  paths  of  motion.  Likewise,  it  is 
shown  in  this  paper  that  if  the  orbiting  rigid  body  is  tri¬ 
axial  or  oblate,  then  again  the  principal  effect  is  that  the 
time-averaged  rotational  angular  momentum  vector  fol 
lows  periodic  paths  of  motion. 

Beletakii^  and  Colombo^  have  studied  the  symmet¬ 
rical  rigid  body  rotational  problem  ( i.  e.  .  two  moments 
of  inertia  are  equal)  in  a  circular  orbit.  Colombo  pub 
lished  two  integrals  of  this  problem  analogous  to  the  first 
two  derived  in  this  paper  except  that  they  refer  to  a  dif¬ 
ferent  coordinate  system. 

This  |uiper  extends  the  work  to  asymmetrical  bodies; 
in  addition,  a  method  is  also  indicated,  but  not  developed 
for  an  extension  to  eccentric  orbits. 

The  more  specialized  problem  of  the  rotational  mo 
tion  of  a  symmetrical  rigid  body  moving  in  :  fixed  orbit 
about  another  body  has  been  treated  successfully  in  the 
astronomical  literature  for  the  case  of  the  rotational  mo¬ 
tion  of  the  earth  under  the  .nfluencc  of  the  gravitational 
attraction  of  the  sun  and  moon;  this  is  the  luni  solar  pre¬ 
cession  problem. 

The  more  general  problem  of  the  motion  of  a  triaxial 
rigid  body  in  a  regressing  orbit  about  another  body,  which 


N  3-5jr  ltv(S)  R  (1) 

where  m  is  the  mass  of  the  primary,  O  is  the  universal 
gravitational  constant,  R  is  the  vector  from  the  center  of 
mass  of  the  small  rigid  body ,  m1,  to  the  center  of  mass 
of  the  primary  mass,  m  (  Fig.  1) . 


EARTHS  POL AR  AXIS 
I  DIRECTION 


IO.X.Y.J)  ARE  CONSIDERED  PARALLEL  TO  INERTIAL  COORDINATES 
0  IS  THE  PRIMARY  m  CENTER  DF  MASS 

ARE  THE  ROTATING  ORBlT  PLANE  COORDINATES 


Figure  1.  Geometry  of  Equatorial  and  Orbit 
Plane  Coordinate  Systems 


The  symbol  R  is  the  magnitude  of  R.  and  (S)  is  the  mo¬ 
ment  of  inertia  matrix  of  m'  with  respect  to  any  particu 
lar  coordinate  system  which  may  be  arbitarily  chosen 
for  the  convenience  of  the  particular  problem  or  analyst. 


From  Newton's  second  law  for  rotational  motion,  the 
torque  given  in  equation  (  1)  may  be  equated  to  the  time 
rate  of  change  of  the  angular  momentum  vector  ( i.  e.  . 


in  an  inertial  coordinate  system. 


If  it  is 


desirable  to  know  the  motion  of  the  angular  momentum 
vector,  in  a  rotating  system,  the  equation  of  motion 
liecomes 


(2) 


V  +  uxL  =  N 
at 

where  u  is  the  angular  velocity  vector  of  the  rotating 
coordinates. 

If  the  rotating  coordinates  are  chosen  to  coincide  with 
the  principal  moment  of  inertia  axes  of  the  small  rigid 
body,  equation  (2)  becomes  Euler's  dynamical  equations. 
These  equations  have  known  periodic  solutions  if  the 
torque  ft  is  zero. 

This  free  Eulerian  motion  is  usually  defined  by  the 

periodic  Euler  angles  ,6  ,  <t>  ,  it  ,  in  a  coordinate 
m  m  m 

system  which  has  its  polar  axis  along  the  constant  vector 
L  (Figs.  2  and  3) . 
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Figure  3.  Geometry  of  Principal  Body  and  Momentum 
Coordinate  Systems 

Now  equation  ( 2)  may  be  interpreted  as  the  pertur¬ 
bation  equation  for  the  variation  of  L  which  is  a  constant 
of  the  motion  in  the  torque-free  Eulerian  motion. 

From  equations  ( i)  and  (  2) 


10,1,1,11  ARE  SAME  AS  DEFINED  III  FISURE  I 

10,  i«.  1.4ml  ARE  COORDINATES  OCFINCO  WITH  |„ 
ALOMO  THE  I  VECTOR  THIS  If 
REFERRED  TO  AS  THE  MOMENTUM 
COORDINATE  STSTEM 


Figure  2.  Geometry  of  Momentum  and  Orbit 
Plane  Coordinate  Systems 

In  the  present  problem,  the  torque  is  not  zero;  hence 
L  is  not  a  constant  of  the  motion. 

It  is  convenient  for  the  purposes  of  this  study  to 
choose  rotating  coordinates  so  that  the  polar  axis  z  is 
normal  to  the  orbit  plane  of  m'  about  m  and  the  x  axis  is 
in  a  direction  parallel  to  the  line  of  nodes  (  Fig.  1) .  In 
this  coordinate  system,  the  angular  velocity  vector  in 
equation  ( 2)  is 


where  fl  is  the  rate  of  change  of  the  right  ascension  of 
the  ascending  node  and  i  is  the  Inclination  of  the  orbit 
plane  to  the  celestial  equator.  These  quantities  will  be 
assumed  to  be  constant  so  that  A  and  B  are  constants. 

A  circular  orbit  will  also  be  assumed  so  that  R  will  be  a 
constant. 


^  =  yHxIS)  R-wxL  (4) 

where  all  the  vectors  and  (S)  are  now  referred  to  the 
rotating  orbit  plane  coordinates  defined  in  Figure  1. 

The  matrix  ( S)  may  be  obtained  by  a  transformation 
from  the  principal  axis  at  m'  to  the  rotating  coordinates 
(  Figs.  2  and  3) 

(S)  =  (T)  (F)  (T) T  (5) 

where  (T)  is  the  transformation  matrix  from  the  prlnci- 

T 

pal  axes  to  the  rotating  orbit  plane  coordinates.  (T)  is 
the  transpose  of  (T),  and  (I')  is  the  diagonal  moment  of 
inertia  matrix  for  m1 

/l|  0  o\ 

(!')  -  I  o  1|  0  I 

\o  0  IS  / 

The  matrix  (T)  may  be  written  as  the  product  of  two 
matrices  (M)  and  (B) .  The  matrix  (M)  is  the  transfor¬ 
mation  matrix  from  the  momentum  vector  coordinate 
system  to  the  rotating  orbit  plane  coordinate  system 
(  Fig.  2) .  The  matrix  (B)  is  the  transformation  matrix 
from  the  rigid  body  principal  axes  to  the  momentum  co¬ 
ordinate  system  (  Fig.  3) . 


k)‘> 


The  matrices  (M)  and  (B)  are 


cos  0  -  cos  6  sin  < 
(M)  sin0  cos  6  cos  i 
0  sin  6 


sin  8  sin  0  \ 
-sls8cos0  j  (6) 
cos  8  / 


And  (B)  - 


•  I  rv*  MOlf4  «OS0  U«f 

in  ni  rr  rf 

11*6  cm 

m  m 


(7) 


The  reason  for  splitting  the  transformation  into  two  ma¬ 
trices  is  to  separate  the  Euler  angles  ,6  ,0  and  0  , 

m  m  m 

which  define  the  orientation  and  Eulerian  motion  of  the 
rigid  body  with  respect  to  the  momentum  coordinate  sys¬ 
tem,  from  8  and  0  which  define  the  position  and  motion 
of  L  with  respect  to  the  rotating  orbit  plane  coordinate 
system. 


Under  the  above  assumption,  the  first  approximation 
of  Kryloff-Bogoliuboff  states  that 

dL  _  L  (t  ♦  T)  -  L  (t)  1  r+  T  dL  . 

dt  3  t -  *7  J  !Tdt 

( 10) 

The  previous  section  was  devoted  to  formulating  the 

differential  equations  in  terms  of  periodic  parameters, 

6_  ,0  ,  0  ,  and  u  from  the  unperturbed  problem  so  that 
m  m  m 

the  averaging  operations  Indicated  by  equation  ( 10)  could 
be  applied. 


The  averages  were  performed  Independently  with 

respect  to  the  angles  ip  ,  0  and  u;  therefore,  no 
m  m 

resonance  conditions  must  exist  and  the  Eulerian 
motion  must  be  nearly  uniform. 


Then, 


(ID 


From  the  preceding  definitions  it  Is  clear  that 

(T)  -  (M)  (B)  (8) 

Combining  equation  ( 4)  with  equations  (  5)  and  ( 8) ,  the 
equation  of  motion  may  be  written 

^  -  ^Rx|M)  |B)  (I1)  <B)T(M)TR-wxL 

(9) 

This  equation  is  now  in  the  form  desired  far  the  applica¬ 
tion  of  the  "averaging  principle"  of  Kiyloff  and 

Bogoliuboff^4  . 

Formulation  of  the  Differential  Equations  for  the  Averag¬ 
ed  Rotational  Motion 

For  the  application  of  the  Kryloff-Bogoliuboff  theory 
to  equation  ( 9) ,  It  Is  necessary  to  assume  that  the  rota¬ 
tional  kinetic  energy  Is  large  compared  to  the  gravita¬ 
tional  potential  energy.  This  is  necessary  so  that  the 

chosen  dependent  variables  (L  ,  L  ,  L  )  =  L  of  the  per- 

x  y  z 

turbed  problem  will  vary  slowly  during  a  period  T  of  any  of 

the  periodic  parameters,  8  .0  .0  ,  or  u  of  the  unper- 
m  m  m 

turbed  problem.  The  symbol  u  is  the  argument  of  latitude 
of  the  orbiting  rigid  body.  This  assumption  is  very  im¬ 
portant  from  a  theoretical  point  of  view  because  it  has 
the  effect  of  separating  librational  theories  for  bodies 
such  as  the  moon  from  rotational  theories  such  as  for 
the  earth.  The  problem  considered  here  is  in  the  rota¬ 
tional  category. 


I  i  III  1  •  I  -S||K,css,i*  •  lHl'i 


sin  n  cos  u 

sin  u  tea  «  •  I,,**#**  w  cos  u  •  u 


(12) 


where  are  elements  of  (8). 
The  average  over  u  Is 


2* 


S,, 


-^  /  R  x  (S)  K  du  -  i-R*[  -SM  J  .  (13) 


The  last  component  is  zero  because  of  the  symmetry  of 

(S). 


Now  equation  (9)  becomes 


The  averages  are  carried  out  over  0  and  0  in  the  same 

ID  ZD 

manrer,  that  is,  by  expressing  8n  and  Sjj  in  terms  of 

.  0„ ,  0„ .  8  and  0  by  equations  ( 5)  ,  ( 6)  ,  (7)  and 
rn  m  in 

( 8)  ,  then  averaging  over  0  and  0 

m  m 


The  result  is 


dL 

dt 


A 

(n  • 


_  A  _  _  _ 

L)  ( 0  x  L)  -u  x  L 
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Where 


K  =  -T5T  (*t  +  4  -  2IJ)  ( t  -  3  C08!Q  )  (15) 

oi\  m 

_  A 

L  Is  the  magnitude  of  L  and  0  is  a  unit  vector  normal  to 

A 

the  orbit  plane,  that  is,  0  =  (0,  0,  1) .  This  is  the  form 
of  the  differential  equation  for  describing  the  averaged 
motion  of  L  (to  be  discussed  In  the  next  section) . 

The  preceding  formulation  may  also  be  used  to  deter¬ 
mine  the  effect  of  an  eccentric  orbit.  The  trigonometric 
terms  and  R  which  involve  u  must  be  written  as  a  Fourier 
expansion  of  the  mean  anomaly  M  and  the  averages  car¬ 
ried  out  over  a  period  of  M. 

The  spin  and  precession  are  assumed  to  be  about  I}* 
in  the  above  definition  of  K.  This  choice  places  no  re¬ 
strictions  on  the  generality  because  no  assumption  has 
been  made  regarding  the  relative  magnitudes  of  the  mo¬ 
ment  of  inertia . 


y 


2AL 


C» 


(20) 


This  is  a  second  integral  of  the  motion.  The  first  inte¬ 
gral  states  that  the  motion  of  the  end  point  of  L  takes 
place  on  the  surface  of  a  sphere.  This  second  integral 
states  that  the  motion  takes  place  on  the  surface  of  a 
parabolic  cylinder.  Therefore,  the  motion  takes  place  on 
the  curve  ( or  path)  formed  by  the  intersection  of  the 
sphere,  equation  ( 18) ,  with  the  parabolic  cylinder,  equa¬ 
tion  ( 20) . 


This  result  is  sufficient  to  infer  that  the  motion  is 
periodic;  however,  the  third  integral  given  below  also 
implies  periodicity. 


From  the  last  equation  in  ( 19)  combined  with  equa¬ 
tions  ( 18)  and  (20) 


z*  -  z  +  C 
A 


Periodic  Solutions 

From  equation  ( 15)  ,  it  is  evident  that 


or 


-Adt 


dz 


(21) 


Therefore 


Integrating  equation  (21)  gives 


L2  =  L2  +  L2  +  L2  =  const,  s  C,  ( 16) 
x  y  z 

which  is  a  first  integral  of  the  motion  and  states  that  the 
averaged  magnitude  of  the  angular  momentum  is  constant 


-At  +  C,  -  / 


dz 


(22) 


Now  equation  ( 16)  may  be  used  to  write  equation 
( 15)  in  terms  of  the  direction  cosines  of  L.  Denoting 
these  by 

£  =  (x,  y.  z)  (17) 


which  is  immediately  recognized  as  an  elliptic  integral 
and  implies  that  z  is  a  periodic  function  of  time. 

This  concludes  the  proof  that  the  averaged  motion  of 
the  angular  momentum  vector  of  a  rigid  body  under  the 
gravitational  perturbation  of  an  oblate  body  is  periodic, 
described  by  the  above  three  integrals. 


Then 


x2  +  y2  +  zl  =  1 


<  18) 


Combining  equations  (17)  ,  (18)  and  (3)  with  equation 
( 15)  ,  the  differential  equations  ( 15)  become 


K  ^ 

*  =  —  zy  +  Az  -  By 
L 

y  =  -  —  zx  +  Bx 


z  =  -  Ax 


(19) 


The  Permanent  Magnetic  Torque.  When  a  rigid  body 
in  orbit  in  the  earth's  magnetic  field  possesses  a  magne¬ 
tic  moment  M  along  L,  a  torque,  given  by 

Li 

_  |~  |  A  a  o  a  a  a  a  I 

Nm  *  ML  Bn  j^2  k*  *  L*  -  ~  (k*  •  0*)  ( 0*  xL*)J 

(  23) 

is  experienced  [5] ,  where  is  the  constant 

magnetic  couple  and  is  a  unit  vector  along  the  Z*  iner¬ 
tial  axis.  An  earth -centered  dipole  representation  of  the 
magnetic  field  is  assumed. 


From  the  last  two  equations  u  Wf'en  < 23)  is  resdved  °°to  the  rotati“« 

orbit  plane  coordinates,  written  out  in  component  form, 

^  B  and  combined  with  equation  ( 15) ,  the  result  is  merely  to 

y  =  zt  -  —  i  change  the  definitions  of  the  constants  A  and  B.  That  is, 

instead  of  equations  (19)  we  have 

which  may  be  integrated  to 
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Observations  of  the  Pegasus  Satellites 


Figure  4  is  a  plot  of  y  versus  z  obtained  from  nu¬ 
merical  integration  of  equation  ( 15)  ,  the  second  integral 
of  motion  given  by  equation  (20)  and  experimental  obser¬ 
vations  obtained  from  the  artificial  earth  satellite  Pegasus 
I  during  the  month  of  May,  1965. 


o  »  «  t  o  *  •  »  •  mi/ 


Figure  4.  Motion  of  Rotational  Angular  Momentum 
Vector  -  Pegasus  1  Satellite  -  Rectangular 
Coordinates  y  vs.  z 
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Figure  5,  Motion  of  Rotational  Angular  Momen¬ 
tum  Vector  —  Pegasus  I  Satellite  —  Spherical 
Coordinates  4>  vs.  0 

Figure  6  is  a  conceptual  sketch  of  the  various  |>aths 
which  the  momentum  vector  may  follow  in  the  orbit  plane 
coordinates.  The  particular  path  which  the  momentum 
vector  takes  depends  on  the  physical  characteristics  of 
the  bod)’  and  the  initial  conditions.  It  may  be  of  practical 
astronautical  significance  to  note  that  both  of  these  condi¬ 
tions  may  be  controlled  passively. 

All  three  of  the  Pegasus  satellites  followed  different 
types  of  paths. 


i 


Figure  5  is  a  plot  of  the  spherical  coordinates  6 
versus  O,  again  obtained  from  the  present  theory,  numer¬ 
ical  integration  and  Pegasus  1, 


Figure  6.  Conceptual  Sketch  of  Possible 
Paths  of  Motion 


The  agreement  with  the  experimental  observations  This  method  has  been  applied  to  the  lum-solur  pro- 

and  numerical  integrations  is  further  evidence  that  the  cession  problem  and  the  results  agreed  with  the  wi  II 

first  order  theory  developed  here  describes  the  physical  known  classical  theory, 
situation  quite  adequately. 
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Attitude  determination  of  the  Ariel  III  satellite 


by 

R.  B.  Bent 

S.R.C.,  Radio  and  Space  Research  Station,  Ditton  Park,  Slouch,  Bucks, 


Ariel  III,  the  third  U.K./U.S.  satellite,  was  successfully  placed  into  orbit 

on  May  5th  1 567  with  apogee  and  perigee  approximately  600  km  and  500  km  respectively 

and  inclination  60°.  It  had  been  suggested  some  time  before  launch  that  the  spin 

axis  direction  required  for  the  interpretation  of  data  from  some  of  the  experiments 

could  be  monitored  by  observing  fluctuations  in  the  telemetry  A.&.C.  recordings  as 

the  satellite  rotated;  this  method  had  been  adopted  for  the  Ariel  II  satellite  ^ 

but  the  analysis  proved  difficult  and  as  the  telemetry  polar  diagrams  for  Ariel  III 

were  found  to  be  nearly  olroular  this  technique  was  regarded  as  unpromising.  An 

(2) 

alternative  method  was  therefore  proposed  similar  to  that  used  on  Telstar  '  '  , 
namely  observing  the  glints  of  solar  reflections  from  the  solar  oell  panels  and 
from  special  mirrors  plaoed  on  the  spacecraft.  Solar  aspect  sensors  already 
inoluded  in  the  payload,  would  determine  the  angle  between  the  spin  axis  of  the 
satellite  and  the  direction  of  the  sun  thus  defining  a  conical  surface  containing 
the  spin  axis.  To  position  the  spin  axis  on  this  oone  another  frame  of  reference 
would  be  required  and  the  mirror  glints  would  define  this.  The  method  was  adopted 
suooessfully  for  Ariel  III  and  spin  axis  direotion  to  better  than  +  2°  has  been 
oaloulated  together  with  coning  angle  and  spin  rate. 

o 

Six  echelon  mirrors  eaoh  with  an  effective  reflecting  surfaoe  at  30  to  the 
spin  axis  direotion  were  placed  around  the  body  of  the  spacecraft  which  was  launched 
with  a  spin  rate  of  30  r.p.sw  When  the  spin  axis  is  correctly  related  to  the 
directions  of  the  sun  and  the  observer  there  are  three  reflections  per  seocnd  from 
the  mirrors,  similarly  the  body  solar  oells  refleot  glints  at  a  rate  of  six  par 
second  with  the  appropriate  directions  to  the  spin  axis  and  the  double  sided  boom 
solar  oells  at  a  rate  of  two  per  seoond.  As  an  example,  with  the  sun's  direotion 

normal  to  that  of  the  spin  axis  there  are  glints  at  angles  of  25,  60,  90  and  155 
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degrees  to  the  scin  ajcis  direction.  The  mirrors  are  rhodium  plated  and  have  a 
reflection  for  normally  inoident  light  in  the  visible  spectrum  whilst  the  solar 
oells  have  a  reflectivity  of  20 ^  for  normal  light.  If  the  mirrors  had  been 
optioally  flat  only  2  or  3  glints  would  have  been  detected  during  a  pass  over  an 
observer  on  the  earth  since  the  sun's  reflection  would  illuminate,  at  any  instant, 
an  area  of  only  3  Mm  diameter,  so  it  was  decided  to  manufacture  them  flat  to  only 
1°  thereby  giving  a  2?'  dispersion  of  the  sun's  image.  According  to  calculations, 
about  ten  7  m.seo.  glints  of  first  magnitude  would  be  observed  at  a  ground  station 
in  a  period  of  about  four  seconds  from  mirror  reflections  whilst  glints  from  the 
solar  cell  panels,  whioh  had  greater  surface  irregularities,  would  eaoh  have  about 
25  m.seo.  duration  with  their  glint  train  lasting  for  at  least  15  sec  and  the  most 
intense  magnitude  glint  being  also  first  magnitude.  As  the  eye  has  an  integration 
period  longer  than  the  individual  glint  duration  these  glints  appear  to  be  about 
two  magnitudes  lower  to  the  eye,  although  photographic  and  photo-electric  devices 
record  the  magnitude  correctly.  An  experiment  carried  out  in  the  laboratory 
showed  that  the  mirror  glints  would  be  easily  visible  to  the  human  eye  through 
binoculars  of  magnification  10  although  reasonably  accurate  predictions  of  the 
azimuth  and  elevation  of  observation  would  be  required.  Owing  to  the  number  of 
reflecting  surft.oe  angles  some  glints  would  be  visible  on  eaoh  pass  over  any  station 
provided  the  sun  »7as  at  least  about  10°  below  the  horizon  and  the  satellite  was 
illuminated  by  the  sun. 

A  number  of  organisations  agreed  to  oo-operate  in  the  observations  and  thqy  are 
listed  below  along  with  the  type  of  equipment  employed. 

1.  National  Aeronautios  and  Space  Administration  (NASA)  -  Worldwide  network  of 

MOTS*  cameras 

2.  Smithsonian  Astrophyaical  Observatory  (SAO)  -  Worldwide  network  of  Baker-Nunn 

cameras 

3.  United  States  Air  Force  (USAF)  -  Baker-Nunn  cameras 

4.  Royal  Observatoiy  Edinburgh  -  Kinetheodolites  at  Edinburgh  and 

Mai  t.a 

5.  Aerospace  Fesearch  Laboratories  U.S.A.F.  -  Photoelectric  camera  at  Ohio 

(ARL) 

6.  N.A.S.A.  Optical  Facility  -  Photoelectric  camera  at  Washington 

•M.O.T.S.  -  Mini  track  Optioal  Tracking  System  17  ** 


Two  teams  of  amateur  observers  controlled  by  the  S.A.O.  Moon  watch  group  and 
R.S.R.S.  were  also  alerted.  As  it  is  difficult  to  prediot  the  position  along  the 
orbit  where  glints  are  expected  at  a  particular  station  the  photoelectric  cameras 
that  are  capable  of  traoking  the  spacecraft  from  h orison  to  horizon  are  the  most 
satisfactory  recording  system. 

echoic’  nincor  rt:ni  Fi  •  ;:'o  J  .-.hov.v  them 

Figure  1  is  a  photo  :mph  of  n  /  mounted  on  the  satellite.  They  oould  be 

adjusted  in  angular  position,  and  alignment  teats  were  carried  out  an  a  flight 
satellite.  The  spacecraft  was  vibrated  to  the  level  expected  during  lift  off 
and  the  resulting  change  in  alignment  was  only  1/6°.  The  final  flight  model 
alignment  was  oarried  out  just  prior  to  launch  when  the  spaoeoraft  doors  had  been 
sorewed  down  for  the  last  time  and  the  full  dispersion  characteristics  of  eaoh 
mirror  were  determined.  The  mirrors  were  found  to  exhibit  different  dispersion 
characteristics,  and  it  is  possible  on  analysis  of  ground  photoelectrlo  photometry 
records  to  indicate  which  mirror  has  led  to  a  particular  flash. 

Slnoe  launch  there  have  been  many  reports  of  glint  observations  from  satellite 
traoking  cameras  and  amateur  observers  all  over  the  world.  A  sample  record  from 
a  Baker-Nunn  camera  in  stationary  mode  is  shown  in  Fig.  3  in  which  the  satellite 
travels  from  right  to  left.  The  camera  shutter  was  opened  just  after  the  satellite 
entered  the  field  of  view.  Some  remarkable  photoelectrlo  records  from  the 
Aerospace  Research  Laboratories,  Ohio  have  led  to  the  oaloulation  of  the  satellite's 
ooning  angle,  which  is  the  angle  between  the  geometrlo  and  dynando  axes,  and  have 

o 

o  _ 

shown  it  to  be  less  than  1  .  There  is  an  r.m.s.  error  in  analysed  spin  axis 
direction  of  about  2°  which  is  no  doubt  caused  by  an  inability  to  identify  the 
peak  magnitude  in  the  glint  train  either  on  photographs  or  during  visual  observations 
by  amateurs.  The  photoeleotrio  data  from  Ohio  however  have  allowed  this  peak 
magnitude  position  to  bs  determined  to  an  acouracy  of  about  1/6  seo.  of  time. 

Thia  type  of  reoording  is  therefore  to  be  reoommended  in  any  future  use  of  the 
optioal  glint  technique.  Such  a  reoorder,  with  azimuth  and  elevatioo  angle 
indicates  oorrect  to  0.1°  enables  the  spin  axis  direction  to  be  measured  to  better 

than  +  1°  and  the  coning  angle  to  0.1°. 
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Furthermore,  a  series  of  measurements 


of  the  glints  is  suiTioisnt  to  measure  thess  angles  without  the  use  of  solar 
aspect  sensors,  particularly  when  more  than  one  set  of  glints  can  be  oba erred 
on  a  single  pass*  By  this  means  the  use  of  telemetry  channels  for  transmitting 
solar  aspect  data  is  avoided. 

The  results  of  spin  axis  direction  during  the  first  115  days  are  shown  in 
fig*  4,  the  day  of  the  year  being  narked  periodically  on  the  curve  with  launoh 
ooourrlng  on  day  125* 

The  eoning  angle  of  the  satellite  which  is  the  angle  between  the  geometric 
axis  and  the  dynamio  axis  around  whioh  the  satellite  spine  has  been  deter  nine  d 
to  be  0.6°  ♦  0.1  °.  As  the  satellite  passes  over  an  observer  it  moves  in  a 
aanner  displayed  in  Fig,  5  where  6  is  the  eoning  angle*  If  the  observer  was 
monitoring  glints  phot oeleetrioally  his  reoonl  should  show  some  modulation  in 
magnitude  due  to  the  faot  that  the  refleoting  surfaces  are  aligned  to  the 
geometric  axis*  Considering  first  the  simple  case  of  reflections  from  the  solar 
osll  pansls  whioh  are  fixed  parallel  to  the  ge  cate  trio  axis  we  have  the  following 
sequemoe  of  events  illus  :  rated  by  Fig*  5*  Whenever  the  satellite  is  in  position 
A  (nose  up),  the  same  panel,  say  P^,  will  be  in  the  position  to  give  specular 
reflections*  The  suoaesslon  of  reflations  from  will  pass  through  a  maximum 
intensity  with  a  given  geomstrioal  arrangement*  Owing  to  the  ocrdng  aotlon 
other  panels  will  give  maxima  at  different  times*  At  the  other  extreme  for 
example,  P^  will  give  maxima  intensities  later  than  Rariel  P^,  and  the  time 
displacement  between  the  conditions  of  maximum  intensity  is  a  measure  of  the 
eenlng  angle  6 . 

Figure  6  shows  the  intensities  of  the  suooessive  reflections  from  all  12 
panels*  From  the  reoord,  taking  every  12th  reflection,  the  intensity  variation 
from  a  single  panel  a  an  be  plotted,  and  Fig*  7  shows  these  variations  In  each 
ease*  The  ourves  marked  6  and  12  relate  to  the  extreme  positions  of  the 
satellite  as  indioated  in  Fig*  5, and  the  time  difference  between  them  enables 
the  ooning  angle  to  be  oaloulated. 

The  dispersion  from  ths  solar  oslls  is  large  ooTpared  with  that  from  the 
60°  mirrors  fixed  to  the  spaoeoraft  and  ths  corresponding  period  of  time  when 
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refleotions  oo our  from  the  mirrors  is  therefor*  shorter*  Similar  but  more 
acourate  results  from  these  slirors  are  displayed  in  Figs.  8  and  9. 

By  considering  reflections  fro*  the  body  solar  oell  panels  it  is  expeoted 
that  monitoring  the  time  at  vhiah  eaoh  panel  refleots  the  sun  with  maximum 
intensity  against  its  solar  oell  panel  nuaber  a  resulting  graph  will  follow  a 
sine  ourre*  Such  a  curve  has  been  plotted  in  Fig.  10  from  the  results  shown  in 
Fig.  7  and  it  is  quits  remarkable  how  aoourately  these  points  follew  a  sins  ourre. 
This  would  indloats  that  the  short  term  scintillation  of  the  solar  reflections  are 
extremely  snail.  Another  example  illustrating  this  point  is  shown  in  Fig.  11 
where  the  magnitude  ▼  time  plot  of  two  adjacent  solar  oell  panels  shows  almost 
perfeot  individual  symmetry.  The  difference  in  magnitudes  between  the  adjaoent 
panels  is  oaussd  by  the  different  number  of  solar  oells  on  alternate  panels  af 
the  spaoearaft  as  o*n  be  seen  in  Fig.  2. 

Thanks  art  due  to  N.A.S.A.,  3.A.O.,  U.S.A.?.  and  the  Royal  Observatory, 
Edinburgh,  for  their  oontlmuing  help  In  monitoring  glints  from  the  satellite,  In 
particular  Drs,  Kissel  and  Vanderburgh  of  the  Aeroepnoe  Research  Laboratories  of 
the  Wright  Patterson  Air  Fores  Base,  Ohio  and  also  to  the  many  amateur  observers 
for  their  valued  asslstanoe.  The  British  Airoraft  Corporation  also  deserve 
special  mention  far  their  s ©-operation  in  oarrying  out  the  installation  end 
meohanioal  tests  or  the  mirrors  in  a  period  of  only  two  months  Just  prior  to  launch. 
Acknowledgements  are  due  to  the  Soienoe  Reaearoh  Counoil  who  have  been  responsible 
for  the  management  and  financial  support  of  the  Ariel  III  international 
oo-operative  projeot. 

Dr.  D.  E.  Smith  of  R.S.R.S,  co-operated  in  the  planning  of  the  experiment  in 
detail  and  developed  t  w  mathematical  techniques  needsd  in  ths  data  analysis.  He 
has  also  organised  the  visual  observations  made  by  many  of  the  amateur  observers. 

The  work  described  was  carried  out  at  ths  Radio  and  Spaea  Reaearoh  Station 
of  the  8olenoe  Reaearoh  Counoil  and  this  pspsr  is  published  with  the  permission  of 
ths  Director. 
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(copy  of  figure  2  not  available) 


Figure  2.  Echelon  mirrors  mounted  on  satellite. 


18  2 


Figure  3  -  Record  of  satellite  traveling  from  right  to  left. 
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Photographed  by  S.A.O.  Baker- Nunn  camera  at  Island 
Laqoon,  Australia,  14  May  1967 


Figure  4  -  Results  of  spin  axis  direction  during  the  first 
115  days 
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^ght  ascension 


Figure  5  -  Illustration  of  a  Sequence  of  Events 
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Figure  6  -  Intensities  of  Successive  Reflections 
from  all  12  panels 
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Figure  7 


-  Plotting  of  intensity  variations 


(Copies  of  Figures  8  &  9  not  available) 


Figures  8  and  9  -  Reflections  from  60°  Mirrors 


SOLAR  CELL  PANEL  NUMBER 


TIME  WHEN  SOLAR  CELL  REFLECTIONS 
REACH  PEAK  MAGNITUDE 


Figure  10  -  Sine  Curve  of  Sun  Reflections 
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REFLECTIONS  FROM  TWO  ADJACENT  SOLAR  CELL 

PANELS  ON  ARIEL  m 

(MAGNITUDE  v  TIME  PLOT  FROM  A.R.L.  PHOTOELECTRIC 
PHOTOMETRY  RECORD  JULY  17  th  1967.) 


Figure  11  -  Curve  showing  magnitude  v  time  plot 
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CZECHOSLOVAKIA  CONTRIBUTION  ON  THE 
SATELLITE  PHOTOMETRY 


F.  Link 

Astronomical  Institute,  Ondrejov,  Czechoslovakia 


Czechoslovakia  Contribution  on  the  Satellite  Photometry 


F.  Link 

Astronomical  Institute,  Ondrejov,  Czechoslovakia 

1.  INTRODUCTION 

Our  work  at  the  Ondrejov  Observatory  (Dept,  of  the  Upper  Atmosphere) 
has  been  concerned  with  the  balloon  satellites,  and  the  photometry  of 
satellite  eclipses  in  order  to  explore  the  upper  atmosphere. 

A  by-product  of  this  work  is  the  study  of  satellite  light  curves  outside 
the  eclipses.  The  form  of  the  light  curves  depends  on  surface  properties 
of  these  satellites  and,  therefore,  I  think  that  some  of  our  results  may 
be  of  interest  for  the  participants  of  this  symposium. 

2.  PHOTOMETRY  OF  SATELLITES 

We  have  used  a  two-color  photoelectric  photometer  in  recording  the 
light  curves  of  the  satellite  transits.  The  scheme  of  each  half  of  the 
photometer  consists  (Fig  1)  of: 

a.  an  achromatic  objective:  50  mm  diameter  and  40  cm  focal  length. 

b.  a  field  stop  C  giving  0.75°  of  the  sky  followed  by  a  Fabry  lens  L 
which  forms  the  image  of  the  objective  on  the  photocathode  FK  of  the 
multiplier  EMI  6094. 

c.  a  wheel  F  containing  the  color  filters. 

d.  a  radioactive  luminiscent  capsule  R  as  the  light  standard. 

e.  a  pen  recorder  with  paper  speed  of  3  mm  per  sec. 

In  the  latest  version  we  have  adopted,  instead  of  two  separate  recorders, 
one  high  speed  Ripp  recorder  micrograph.  The  photocurrents  of  both 
photometers  are  connected  alternately  by  an  astronomical  clock  (each  three 
seconds)  to  this  recorder;  in  this  manner  we  obtain  two  superposed  light 
curves  whose  interruptions  are  exactly  timed. 

The  whole  photometer  has  an  altitude/  azimuth  mounting.  The  tracking 
is  performed  manually  with  the  observer  using  both  hands  -  one  hand  acting 
to  push  or  pull  in  azimuth  and  the  other  to  set  altitude  by  soft  striking  of 
a  guiding  bar.  With  some  training  it  is  possible  to  track  all  the  balloon 
satellites  by  this  method. 
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3.  ROTATION  OF  SATELLITES 


All  three  balloon  satellites  Echo  I,  II  and  Pageos  I  display  light 
fluctuations  which  may  be  interpreted  in  terms  of  rotation.  On  the 
light  curves  one  can  find  the  repetition  of  many  details.  Their  time 
separation  gives  directly  the  length  of  the  photometric  period,  Tph» 
which  is  of  course  slightly  different  from  the  sidereal  period,  Tsid» 
because  of  the  combination  of  rotation  with  the  topocentric  revolution 
of  the  satellite. 

In  the  most  simple  case,  when  the  orbit  of  the  satellite  is  placed 
in  the  plane  containing  also  the  sun  and  the  observer  and  the  spin  axis 
is  perpendicular  to  this  plane,  we  have  the  following  relation  (2,3): 


W1  ± 


A 

720° 


) 


where  A  is  the  topocentric  angle  described  by  the  satellite  during  one 
photometric  rotation.  It  is  assumed  that  the  balloon  surface  produces 
a  specular  reflection  and  that  the  brightness  fluctuations  result  from 
local  irregularities  in  the  surface  figure. 

In  the  above  relation  the  topocentric  angle  A  increases  with  decreasing 
topocentric  distance  of  the  satellite,  i.e.  ,  toward  the  zenith.  Our  formula 
gives,  therefore,  the  possibility  to  fix  the  sense  of  the  rotation  by  the  use 
of  the  following  rule: 

"If  the  photometric  period  of  rotation  ir. creases  toward  the  culmination, 
the  sense  of  rotation  is  direct,  i.e.,  identical  with  the  sense  of  revolution, 
and  vice  versa.  " 


4.  EXPERIMENTAL  RESULTS 


The  main  results  obtained  for  the  three  balloon  satellites  are  as 
follows: 

ECHO  I  (I960  Iota  I) 

The  light  curves  have  a  very  irregular  shape  with  many  periodic  details, 
enabling  good  determinations  of  the  rotation  period.  From  November  1964 
to  February  1967  the  length  of  the  period  has  dropped  from  0.  93m  to  0.  66m 
(2,  4).  The  internal  accuracy  of  these  determinations  is  as  high  as  +0.  003m  + 
0.  2s;  the  main  uncertainty  is  due  to  the  systematic  differences  between  the 
photometric  and  sidereal  periods  which,  of  course,  change  from  one  transit 
to  another. 

The  form  of  the  light  curves  seems  to  be  determined  mainly  by  the 
relative  position  of  the  three  bodies:  sun,  satellite,  observer,  since  we 
have  found  nearly  the  same  form  of  the  light  curve  (Fig  2)  after  an  interval 
of  5  months.  This  shows  also  the  constancy  of  the  satellite  form  during 
this  time  interval.  The  sense  of  the  rotation  is  found  to  be  retrograde. 

In  fact,  we  found  for  the  differences  of  the  photometric  periods  (horizon  - 
zenith)  the  following  results: 
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in  41  cases 
in  29  cases 
in  16  cases 


+  0.  414m 
-  0.  129m 
zero 

ECHO  II  (1964  -  04A) 

In  1964  we  observed  many  interesting  light  curves  showing  (as  in 
Fig  3)  the  existence  of  discontinuities  due  to  antennae  or  similar  systems 
on  the  opposite  sides  of  the  balloon.  From  September  1964  to  August  1965, 
the  rotation  period  had  slowed  from  1.  52m  to  about  2.  3m.  After  this  date 
a  lack  of  details  on  the  light  curves  makes  the  determination  of  the  rotation 
quite  difficult. 

PAGEOS  (1966  -  56 A) 

The  light  curves  always  show  many  details,  but  their  repetition  is  not 
obvious,  in  contrast  to  Echo  I.  In  the  summer  of  1966,  the  sidereal  rota¬ 
tion  period  was  between  12-13m(4).  Some  of  the  autumn  1966  curves  dis¬ 
play  a  series  of  light  flashes  of  great  amplitude  (Fig  4).  At  the  end  of 
March  1967,  the  period  seems  to  have  fallen  to  about  6m. 

In  every  case  the  use  of  Pageos  as  a  geodetic  satellite  would  not 
have  been  too  easy  because  of  the  large  fluctuation  of  the  light. 

5.  ECLIPSES  OF  ECHO  2 

As  the  eclipses  of  the  satellites  lie  outside  the  scope  of  this 
symposium,  let  me  mention  only  briefly  some  work  accomplished  at 
Ondrejov  Observatory. 

The  photometric  theory  of  satellite  eclipses  had  been  already  developed 
in  1961  (5,6)  on  the  basis  of  our  30-year-old  lunar  eclipse  theory.  The 
computation  of  an  ephemeris  supplements  this  work  (7).  Finally,  to  aid 
in  solutions  of  satellite  eclipse  problems,  we  have  published  a  set  of  very 
convenient  Tables  (8)  giving  the  density  of  the  shadow  for  Echo  II  eclipses 
in  every  possible  case  between  30°  N  and  70°  N  latitude. 

Up  to  the  present  time  we  have  measured  several  tens  of  eclipse 
curves  of  Echo  2.  Their  final  reduction  depends  on  the  exact  orbital 
elements  which  will  soon  be  (we  hope)  available  from  the  Smithsonian 
Astrophysical  Observatory.  Some  preliminary  attempts  in  this  direction 
on  two  eclipses  are  very  promising  (9). 
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FIGURES 


Figure  1 
Figure  2 

Figure  3 

Figure  4 


Two-color  photoelectric  photometer. 

Light  curves  of  Echo  1  in  5-month  interval.  Similarity 
of  both  curves  but  change  of  the  period  (0.  93m-0.  82m). 

Light  curves  of  Echo  2 

3-4Mid  September  1964  (rot.  per.  1.  57m) 

47-49  end  of  October  1964  (rot.  per.  1.  68m) 

72-79  beginning  of  May  1965  (rot.  per.  2.  0m) 

Portion  of  Pageos  light  curve  (1966  X.  22)  showing  several 
flashes. 
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Figure  2.  Light  curves  of  Echo  1  in  5-month  interval.  Similarity 
of  both  curves  but  change  of  the  period  (0.  93m -0.  82m). 
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Figure  3  -  Light  curves  of  Echo  2 

3-5  Mid-September  1964  (rot.  per.  1.57m) 
47-49  end  of  October  19o4  (rot.  per.  1 . 68™) 
72-79  beginning  of  May  1965  (rot.  per.  2.0™) 
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Figure  4.  Portion  of  PAGEOS  light  curve  (1966  X.  22)  showii.g 
several  flashes. 
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ABSTRACT 


Photoelectric  photometry,  displayed  as  continuous  lightcurves  on  a 
logarithmic  scale,  are  presented  with  sufficient  detail  to  enable  the  reader 
to  ascertain  their  phenomenological  character  and  general  quality.  Selected 
portions  of  these  lightcurves  are  presented  and  discussed  in  greater  detail 
to  enhance  the  understanding  of  the  physical  processes  important  in  the 
optical  scattering  of  spacecraft  and  in  the  collection  of  data.  Calibration 
techniques  and  data  reduction  methods  are  outlined,  with  commentary  on 
the  significar  f  short  and  long  term  oiiguiness  liuctuations.  A  portion 
of  the  data  is  analyzed  and  attention  is  drawn  to  the  types  of  information 
which  is  extractable  and  which  may  be  of  interest  to  geodecists,  PAGEOS 
project  scientists  and  engineers,  fabricators  of  future  balloon-satellites, 
and  investigators  of  satellite  dynamics.  The  earliest  records  indicate 
that  a  ^pheroid  of  about  fifteen  meters  radius  was  deployed,  the  radius 
being  nearly  uniform.  The  most  recent  record  strongly  suggests  a 
shrinking  of  the  original  sphere  to  a  prolate  spheroidal  shape  during  its 
fifteen  months  in  orbit.  Intermediate  records  indicate  varying  but 
generally  increasing  degrees  of  surface  anomaly.  These  data  compirise 
the  most  extensive  and  definitive  information  at  optical  wavelenglns 
known  to  exist  for  the  PAGEOS  satellite,  and  were  collected  at  the 
Aerospace  Research  Laboratories'  Sulphur  Grove  (Ohio)  Observatory. 
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INTRODUCTION 


As  part  of  the  program  on  photoelectric  photometry  of  orbiting 
spacecraft  carried  out  at  the  Aerospace  Research  Laboratories  since 
early  1962,  the  PAGEOS  spacecraft  was  singled  out  for  study  for 
several  reasons.  These  include  the  intent  (1)  to  study  the  spin  history 
from  deployment  until  a  condition  of  steady  state  was  reached,  (2)  to 
study  the  scattering  properties  of  the  object,  both  diffuse  and  speculer 
if  possible,  for  comparison  with  those  of  Echo  I  (its  prototype)  and 
Echo  II,  (3)  to  study  the  buildup  of  brightness  fluctuations  following 
relaxation  of  the  internal  inflation  pressure,  and  (4)  to  use  this  object 
for  measurement  of  eclipse  phenomena  as  measured  near  its  apogee  for 
comparison  with  eclipse  measurements  of  the  nearer  Echos  I  and  II. 
Another  inflation-deployed  spacecraft,  OVl-8  (1966-63A)  was  scheduled 
for  launch  the  following  month  and  was  also  to  be  observed  from  soon 
after  initial  injection  for  the  same  reasons. 

The  anticipation  of  a  spin  build-up  was  predicated  on  the  previous 
optical  measurements  at  ARL  and  at  Ondrejov,  Czechoslovakia  by  F.  Link 
(Ref.  1)  of  inflatable  spacecraft,  namely  Echo  I  (19b0  Ioto  1),  Explorer  9 
(1961  Delta  1),  DASH  2  (1963-30D),  Explorer  19  (1963-53A),  Echo  II  (1964- 
04A),  and  Explorer  24  (1964-76A),  all  indicated  that  a  state  of  rotational 
motion  at  surprisingly  rapid  rates  was  the  eventual  equilibrium  state  for 
all  of  these  vehicles  although  the  rotation  of  Echo  II  has  dropped  to  a  very 
low  level.  Table  I  gives  rotation  rates  inferred  from  optical  photometry 
on  these  spacecraft  on  various  dates  of  observation  when  definite  periodic 
brightness  variations  were  discerned.  Included  are  data  on  OVl-8.  Early 
Echo  II  data  based  on  radio  measurements  of  the  Goddard  Space  Flight 
Center  are  also  included,  as  reported  by  Mar  and  Vigneron  of  the  Defense 
Research  Telecommunications  Establishment  of  Canada. 

Photoelectric  data  on  the  PAGEOS  spacecraft  were  obtained  as  early 
as  6620m  after  inflation,  i.  e.  ,  after  2.  1  orbits  of  the  earth  near  the 
beginning  of  Rev  3.  Additional  early  data  were  obtained  on  Rev  11  and 
Rev  27.  These  data  were  obtained  at  the  ARL  OPOS  Observatory  at 
Sulphur  Grove,  Ohio.  Observing  equipment2  consisted  of  a  4-axis 
mounted  24"  Cassegrain  telescope  shown  in  Figure  1  employing  a  loga¬ 
rithmic-feedback  photometer3  using  a  specially  selected  RCA  7029  photo¬ 
multiplier  which  has  an  S-17  photocathode.  The  detector  was  used  without 
filter;  the  resulting  broad-band  photometry  thus  covered  the  wavelength 
range  from  310  to  600  millimicrons.  Data  were  recorded  on  an  FM  tape 
recorder.  Aided  manual  tracking  maintained  the  target  in  a  two  arc- 
minute  diameter  aperture.  Stars  were  recorded  for  brightness  compari¬ 
sons,  and  reference  voltages  recorded  as  a  brightness  scale. 

Selected  oscillograms  are  shown  in  Figures  2-7.  These  are  time- 
compressed  copies  produced  by  a  light-beam  galvanometer  recorder,  and 
read  from  right  to  left.  One  data  channel  records  photometer  output, 
another  time  pulses  recorded  directly  from  the  WWV  radio  transmissions. 


203 


Table  1.  Rotation  Rates  of  Inflation -Deployed  Spacecraft  Data  are  apparent  rotational  periods 
derived  by  optical  photometry  at  the  ARL  OPOS  Observatory  unless  otherwise  noted. 
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Link,  B.A.C.  Vol.  12  (1966)  No.  1  (Reference  1). 

Mar  and  F.  R.  Vigneron,  COSPAR  - 1 AU -IUT  AM  Symposium  of  20-23  April  1965  (Reference  9). 
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Figure  2 1  PAQEOS  Lightcurre*  21i  June  1966  -  11  July  1966 
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Flgnr*  3i  PAG£OS  LightcurvM  16  July  1966  -  21  July  1966 
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Figure  5 1  PAQBOS  Lightcurrea  22  February  1967  -  25  July  1967 


U  August  1967  Oh?7 

Figure  6»  PAGEOS  Lightcurree  1  July  1967  -  1*  August  1967 


Figurs  7<  PAGEOS  Lightcurvss  5  August  1967  -  25  August  1967 


Oci  .iyion.il  brightness  drops  between  satellite  and  dark  indicate  sky 
brightness  as  accepted  by  the  field  aperture.  Deeper  drops  indicate 
the  dark  current  level  of  the  photomultiplier.  In  the  earlier  records 
these  are  effected  at  .10- second  intervals  by  a  timer. 

Brightness  Calibration 

A  star  calibration  technique  is  used  to  measure  photometer  response 
and  atmospheric  transmission.  Selected  stars  are  observed  photometri¬ 
cally,  and  the  corresponding  output  voltage  (deflections)  recorded. 

Figure  8  is  a  plot  showing  the  results  of  star  calibration  conducted  on 
27  June  19t>6.  The  groups  of  stars  defining  the  zenith  curve  were  observed 
near  the  zenith,  and  their  published  visual  magnitude  is  plotted  versus 
photometer  output  voltage,  with  corresponding  corrections  applied  to 
the  measured  output  voltages  to  account  for  these  small  zenith  distances. 
This  curve  thus  represents  output  voltage  as  inferred  through  one  air- 
mass,  as  a  function  of  extra -atmospheric  visual  stellar  magnitude. 

Values  for  the  visual  stellar  magnitudes  have  been  taken  from  recently 
published  catalogs4’5  prepared  from  modern  conventional  photoelectric 
photometry.  G-type  stars  (matching  the  sun's  energy  distribution)  were 
used,  although  for  the  bright  region  only  the  blue  star  Vega  was  accessible 
at  this  sidereal  time.  In  this  case,  a  color  correction  was  applied  to 
obtain  a  simulation  of  a  G  star  of  zero  magnitude.  From  this  approach, 
absolute  brightness  is  measurable  to  .  05  to  .  5  magnitudes  depending 
upon  observed  brightness  and  sky  quality.  The  observed  brightness 
range  for  all  the  recorded  PAGEOS  transits  lies  along  the  linear  portion 
of  the  response  curve,  so  short  term  fluctuations  were  subject  to  mini¬ 
mum  measurement  error.  On  Revs.  27  and  3566  conditions  were  excellent 
and  the  absolute  errors  believed  less  than  0.1  magnitude. 

Analysis  of  the  Brightness  Data 

Physical  interpretation  of  brightness  data  requires  separation  of  the 
received  brightness  into  specular  and  diffuse  components.  While  pre¬ 
launch  laboratory  measurements  by  NASA6  estaolished  the  initial  values 
of  the  scattering  components  quite  precisely  (See  Table  II), 


Table  II.  Initial  Reflectivity  of  PAGEOS 


R  ,  .  0.  862  +  0.  0046 

specular  — 

R  .  0.  029  +  0.  0046 

uiffuse  — 

R  .  0. 891  +  0. 0032 

total  — 

These  values  are  weighted  to  represent  the  solar  spectrum  within 
the  range  3,  600  A  to  7, 000  A.  (Reference  6) 


there  is  sufficient  uncertainty  concerning  effects  of  the  orbital  environ¬ 
ment  to  warrant  verification,  especially  after  the  satellite  has  been  in 
orbit  for  some  time.  This  is  possible  for  an  approximately  spherical 
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satellite  by  comparison  of  selected  average  brightness  values  over  a 
wide  range  of  phase  angle  (the  angle  at  the  satellite  formed  by  sun  and 
observer  vectors).  Given  a  suitable  change  in  illumination  geometry, 
a  diffuse  contribution  exceeding  5%  should  be  detectable  from  the  re¬ 
corded  brightness  data.  Because  of  the  short-term  brightness  fluctua¬ 
tions,  primarily  random  in  nature,  characteristic  of  most  of  the  data, 
it  is  difficult  to  determine  localized  average  brightness  values.  How¬ 
ever,  experimental  conditions  for  the  transits  yielding  the  27  June  19 1> 6 
and  30  July  1967  records  cover  a  large  range  of  phase  angle  and  occurred 
on  nights  of  goo A  photometric  quality.  Reduction  of  these  data  shows 
that  the  diffuse  contribution  was  less  than  5  percent  on  Rev.  27  and 
has  remained  at  least  that  low  during  the  first  13  months  in  orbit.  It 
seems  reasonable  to  infer  that  since  specularity  has  remained  high, 
there  has  been  little  microscopic  surface  degradation,  i.e.,  a  dulling 
of  the  shiny  surface  would  have  increased  the  diffuse  component.  This 
bound  on  the  constancy  in  reflectance  properties  is  in  accord  with  the 
conclusions  reached  by  Emmons,  et.al.  ,  based  on  strongly  averaged 
data  after  2  months  in  orbit  and  inferred  by  defining  a  specularity  index 
as  the  ratio  of  specular  to  specular  plus  diffuse  scattering. 


Table  III  presents  positional  parameters,  absolute  brightness,  and 
inferred  radii  of  curvature  as  functions  of  time,  for  selected  transits. 
The  positional  parameters  were  computed  from  Smithsonian  orbital 
data  and  the  mean  anomaly  adjusted  by  local  position  checks  generally 
by  photoelectrically  observed  near  occultations .  Observed  brightness 
was  converted  to  absolute  brightness  from  the  usual  relation10 
R  2 

M  =  m  -2-5  l°g(-[QQQ  ~)  where  M  is  the  absolute  magnitude  for  1000 


kilometers  slant  range,  and  m  the  observed  magnitude  (corrected  for 
extinction).  Radii  of  curvature  were  calculated  by  simultaneous 
application  of  the  specular  reflection  and  diffuse  scattering  expressions:1  1 

E  r  b2 


E 


sp 


s  sp 
4R  2 


(3) 


2E  r  b2 

E  ,  =  - - -  [  sin  0  +  (tt-0)  cos  0]  (4) 

d  3 it  R2 

where  E  is  the  incident  solar  flux1  2  13.  21  x  104  (lumens/ m2),  b  the 
radius  of  the  satellite  (meters),  R  the  distance  from  observer  to 
satellite  (meters),  0  the  phase  angle  (formed  at  the  satellite  by 
observer  and  the  sun  vectors),  Egp  the  specular  contribution,  E^  the 
diffuse  contribution  of  satellite  brightness  (lumens/ m2),  and  rgp 
and  r^  the  specular  and  diffuse  reflectivity,  respectively.  Using 
the  laboratory  data  of  Table  II  for  the  specular  and  diffuse  reflectivity: 
i.e.  ,  assuming  only  slight  departures  from  the  predeployment  values: 
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Table  III.  Reduction  of  PAGEOS  Photometric  Data  to  Absolute  Magnitude  and  Radius  of  Curvature, 


CO 

> 

<u 

Pi 


X 

x 

in 

CO 


T3 

s 

IT) 


r- 

INI 


W 


> 

0) 


X 

X 

D 

3 

3 

M* 

OJ 


m 

3 

T> 

(T) 

OJ  ' 

4)  — 

Ih  (0 

oj  ^ 

a 

CO  " 

c 

D 

hi 

Ifl 

a 
a 
< 


x 

IT) 

s 


e 

o  —  .5 
§ 


-  s 

5  w 

So 

it!  5- 

3 

in 


>  s 


4) 

'3) 

3 

< 

<U  1 

10 

it) 

X 

Hi 


e 

<u 

M 

rt  § 


N 

< 


H 

D 


o 

a 


t  M)  r-  i/i  co 

— *  O 

Is-  Is-  oo  oo  Is-  i- 

in!  in  mi  m3 

M*  O'  oo  oo  in 

m  so  in  co 

cO  O  --  O  O'  O' 

V  r-  M3  X 

<\]  CM  f\J 

ro  — <  tn  (M  O 

r-  o  o  im 

CO  CO  V  CO  (M 

in  —  o  — ■ 

II  till 

o  — *  — * 

i  i  i  i 

tf  o  0  oo  r- 

O  IM  x  O 

O'  m3  n-  m3  in  in 

O'  co  — •  M* 

1  1  1  1  1  1 

O  ^  ^ 

1  1  1  1 

— < 

> 

0  0  O  0  0  0 

«  0  O  0  0 

X  — •  X  O  M*  Is- 

K  o-i  oo 

O  — i  — <  IM  IM  CM 

o  — •  — 1  (M 

O  O  O  O  O  O 

o  o  o  o 

N  I-  Is  O  O  ^ 

O  X  M3  O 

m  m*  tj<  m  md  r- 

co  —>  (M  co 

m  m  m  m  m  m 

in  m  x  in 

0  0  O  0  0  0 

0  0  0  0 

Tf  m  in  m  co  o 

—1  CO  if  C 

oo  Is-  x  m  i** 

0s  r»  LO  If 

o  o  o  o  o  o 

o  o  o  o 

o 

o 

o 

o 

o 

o 

0  0  0  0 

'O  t--  n-  t—  m  co 

— 1  CO  (M  — < 

CO  CO  CO  CO  CO  CO 

It*  M*  ^ 

in  cm  m  o  oo  m 

N  ^  ro 

m  m  m  ro  m 

co  in  ^  m 

on  <v>  in 

x  — i  in  x 

m  m  o  o  o  o 

M3  o  — 1  — 1  — 1 

n-  r-  oo  oo  oo  oo 

m3  X  X  X  X 

o  o  o  o  o  o 

o  o  o  o 

0 

3 

3 

— i  CM  co  M*  X  X 

— 1  (M  CO  Iji 

m 

#***•«•* 

IM  *  #  *  ft 

Is- 

fM 

> 

0) 

Pi 


X 

X 

V 

3 

3 

h-J 

(M 


00  O  O'  —  X 


I 


in  co  'r  m 


xxxo'xxxxxco 

i^co'r^,r^a0'O'OTt<r- 


X  O'  (M  M"  M3 

O'  O'  00  00  O 

•  •  •  •  • 

o  o  o  o  — 

I  I  I  I  I  I  I  I  I  I 


00  O  CO  00  M3 
CO  00  O'  O'  CO 


-•  o  o  o  -• 

I  I  I  I  I 


in  O'  in  m  oo 

TJ<  — I  (M  O'  IM 


I  I  I  I  I 


oooooooooo 

M'mo'ooo'in— <M30'— > 
rs-rs-r-ooo'0— •— <  —  fM 


o 

o 

o 

o 

o 

o 

O 

O  O 

o 

in 

IM 

M* 

00 

— 4 

00 

M* 

r-  rt* 

PO 

00 

r- 

— M 

00 

r- 

00 

O'  o 

X 

X 

X 

X 

M* 

M* 

Tf  X 

m 

oooooooooo 

O'ooNO'Ncciintom 

cococo— i  o'  r-  m  ij1  iji  co 

__i  — i—iOOOOOO 


oooooooooo 

coinNcor^^iMO't'in 

NiMo^inmin^ii'i)1 


X 

M*  O  X 

O 

o 

M* 

M* 

o 

-H 

o 

o 

— * 

M* 

o 

O 

o 

m 

6 

M* 

X  X 

X 

X 

O 

(M 

CM 

(M 

CM 

CO 

ro 

rn 

rr 

M* 

00 

00 

X 

X 

CO 

00 

X 

X 

00 

00 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CM 

CO 

in 

mO 

X 

O' 

«■ 

* 

* 

* 

* 


214 


Indicates  brightness  enhancement  point 


Table  III  (continued) 
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Table  III  (continued) 
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.  029  and  recalling 


r 

sp 


rd 


that  absolute  stellar  magnitudes  are  extra  atmospheric  values  at  1000 
km  distance, 


E  b2 

E  =  - - -  [  0,  216  +  0.0062  (sin  0  +  [ tt - 0 J  cos  0)]  (5) 

101  2 


where  E  is  the  total  brightness  contribution. 


or 


_ 7.  56  x  106  E _ 

0.  216  +  0.0062  [  sin  0  +  (tr-0)  cos  0] 

The  conversion  of  absolute  stellar  magnitude  to  lumens/  m2  is 
accomplished  by  the  expression: 


2.  65  x  10~6 
2.  512M 


(7) 


Thus  from  the  measured  value  E  and  the  computed  value  0  the  radius 
of  curvature  could  be  calculated,  assuming  constant  reflectivity.  The 
absolute  values  of  the  calculated  radii  of  curvature  are,  of  course, 
only  as  good  as  the  brightness  measurements  and  the  degree  to  which 
the  assumed  reflectivity  coefficients  represent  the  optical  properties 
of  the  satellite  surface.  Figure  9  gives  a  graphical  solution  of  Eqns. 

(6)  and  (7)  in  the  range  of  brightness  of  importance  in  this  paper,  where 
0  was  essentially  constant.  From  Figure  9,  it  may  be  noted  that  an 
absolute  brightness  uncertainty  of  0.05  magnitudes  produces  a  radius 
of  curvature  uncertainty  of  0.  3  meters  in  the  neighborhood  of  the 
nominal  15  meter  radius.  In  this  same  region,  a  brightness  uncertainty 
as  large  as  0.  5  magnitudes  would  produce  an  apparent  radius  of  curva¬ 
ture  uncertainty  of  3.7  meters.  However,  as  noted  above,  short  term 
relative  brightness  fluctuations  are  independent  of  absolute  values,  and 
thus  their  measurement  carries  the  minimum  error.  The  absolute 
brightness  values  in  Table  III  for  24  and  25  June  1966  are  probably 
subject  to  errors  up  to  +  0.  5  magnitudes  as  the  result  of  poor  skv 
conditions;  for  the  other  three  dates,  tabulated  absolute  magnitudes  are 
almost  certainly  good  to  +_  0. 1  magnitudes  or  better. 

Rotational  History  of  PAGEOS 

Rotation  of  the  satellite  about  its  center  of  mass  is  detectable  from 
some  of  the  brightness  data.  Figure  10  is  a  plot  of  apparent  spin  rates,  and 
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indicates  that  PAGEOS  has  spun  up  from  a  low  rate  of  about  0.  083 
revolutions  per  minute  in  the  fall  of  1966  to  a  high  of  about  0.  325  rpm 
in  the  summer  of  1967.  The  first  three  points,  showing  an  apparent 
fast  spin  at  launch  are  based  on  brightness  fluctuations  whose  periodic¬ 
ity  may  be  multiples  of  the  true  period  although  these  points  join  the 
others  smoothly.  Periodic  patterns  on  later  records  were  sufficiently 
detailed  to  allow  determination  of  spin  rates  with  a  high  degree  of 
confidence.  Two  additional  points  were  obtained  from  independent 
optical  measurements,8  and  as  the  plot  shows,  one  is  in  close  agreement 
with  the  ARL  data  while  the  other  is  slightly  off  the  apparent  curve. 

Since  the  discrepancy  occurs  at  a  time  when  the  period  was  longest,  it  is 
most  likely  due  to  differences  in  observing  geometry,  i.e.,  use  of 
apparent  period  to  infer  sidereal  period  plus  contributing  shifts  of 
reference  points  for  the  apparent  rotation.  In  any  case,  it  is  evident  that 
the  rotational  period  did  decrease  from  a  maximum  of  about  7  30  seconds 
per  turn  to  a  minimum  of  about  195  seconds  per  turn  in  245  days.  This 
represents  a  change  in  kinetic  energy: 

ltd*  IW; 

AK  =  - —  -  - - -  (1) 

2  2 

where  I  is  the  moment  oi  inertia,  ubthe  minimum  angular  velocity  and 
U>2  the  maximum  angular  velocity.  Tor  PAGEOS  I  =  2/3  (mass)' 
(radius)2  =  8375  Kg-m2;  cd,  =  8.6  x  10"3  radians/  second  =  3.  22  x  10  2 
radians/  second.  Thus  AK  =  4.04  joules.  The  average  torque  required 
to  produce  this  energy  change  is 

T  =  I  a 

2.  36  x  10  2  radians/  sec  (2) 

2.  12  x  107  seconds 

=  9.  30  x  10  6  newton  -meters 


u>2  ' 

where  a  =  - 

245  days 


These  indicated  small  values  for  energy  change  and  torque  suggest 
that  rotation  of  PAGEOS  about  its  center  of  mass  and  the  observed 
rate  changes  could  have  been  caused  by  solar  radiation  pressure, 
the  earth's  magnetic  field,  or  both.  The  analysis  of  Echo  II  rotation 
in  Reference  9  anticipates  a  radiation  pressure  torque  of  8  x  10"6 
newton  meters  due  to  thermal  shape  distortions. 


In  the  above  calculation  a  mass  value0 


of  54.  08  kg  was  used. 


220 


Determination  of  the  PAGEOS  Figure 

The  light  curve  of  15  September  1967,  showing  the  largest  and  most 
periodic  brightness  fluctuations  yet  recorded,  was  selected  for  detailed 
reduction.  Figure  1 1 -A  shows  the  entire  light  curve  at  approximately 
the  same  scale  as  in  Figures  2  through  7.  The  obvious  periodicity, 
shown  in  the  expanded  section  (Figure  11  -B),  results  almost  certainly 
from  the  large  specular  brightness  component.  The  specular  highlight 
which  contributes  to  the  specular  component  is  only  an  area  subtending  a 
0.  25-degree  arc  at  the  center  of  the  sphere,  i.  e.  ,  an  area  some  7  cm  in 
diameter.  This  implies  that  the  brightness  values  represented  by  the 
section  AA  (one  rotation)  were  produced  by  the  virtual  image  of  the 
sun  as  it  mapped  out  a  th  i  band  nearly  circumscribing  the  satellite,  i.  e. 
a  trace  some  4-8  cm  wide  and  96  meters  in  path  length  which  nearly  joins 
back  on  itself  as  the  satellite  rotates  under  the  solar  image.  Figure 
ll-C  is  a  polar  plot  of  corresponding  radii  of  curvature,  assuming  con¬ 
stant,  pre-launch  reflectivity  values.  A  first-order  smooth  joining  of 
these  curvatures  produces  the  apparent  cross  section  of  Figure  11  -D. 

At  this  scale,  the  small  curvature  fluctuations  between  the  deep  minima 
are  not  discernible,  although  the  measured  data  definitely  show  a 
preponderance  of  small  surface  anomalies  on  one  side  as  compared 
with  the  other.  The  prolate  spheroidal  shape  that  this  cross  section 
implies,  holds  only  if  the  observe  band  is  close  to  being  the  maximum 
circumference  of  the  actual  figure.  Cross  sections  of  less  peripheral 
measurement  (whose  planes  do  not  pass  through  the  figure  center)  imply 
greater  localized  curvature  distortions  to  be  consistent  with  observed 
radii  of  curvature. 

There  are,  of  course,  other  possible  interpretations  of  the  brightness 
fluctuations  recorded  on  September  15th.  However,  those  that  come  to 
mind  seem  less  probable  than  the  prolate  spheroid  model.  If  the  satellite 
is  spherical,  the  deep  brightness  minima  might  be  accounted  for  by  loss 
of  specularly-reflective  surface  area  from  regions  on  the  sphere  sub¬ 
tending  nearly  60  degrees  at  the  center  of  the  sphere  and  opposite  each 
other,  or  by  loss  of  reflectivity  due  to  skin  wrinkling.  It  seems  unlikely 
that  the  first  of  these  mechanisms  for  reduction  of  specular  scattering 
would  obtain  locally  rather  than  generally  over  the  entire  sphere  since 
it  involves  removal  or  degradation  of  the  aluminized  layer  on  a  sub- 
macroscopic  scale.  The  second  of  these  is  certainly  admissible  if  one 
postulates  two  opposed  spherical  caps  with  progressively  wrinkled 
zones  as  the  virtual  solar  image  approaches  the  poles  of  the  caps.  It 
is  felt  beyond  the  scope  of  this  paper  to  analyze  the  nature  of  the  wrinkling 
to  yield  a  smoothly  increasing  change  from  a  scattering  distribution  of 
a  purely  specular  form  to  one  which  is  quasi-diffuse  in  form.  The 
difficult  concept  is  the  statistical  distribution  of  the  macroscopic  specular 
elements  which  must  redirect  the  light  in  a  diffuse  distribution  such  that 
another  zone  of  wrinkles  does  not  reinforce  the  central  virtual  image. 

The  end  caps  also  must  subtend  an  included  angle  of  some  60  degrees, 
i.e.  they  are  16  meters  in  diameter  much  larger  than  the  one-meter 
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diameter  polar  caps  of  the  PAGEOS  construction  which  are  the  areas  of 
unique  character  on  the  vehicle  where  stresses  might  develop  at  opposite 
points  on  the  nominal  sphere. 

This  wrinkle  mechanism  is  much  more  complex  and  constrictive  in 
its  character  than  the  assumption  that  the  vehicle  has  distorted  to  a 
spheroid,  but  it  cannot  be  ruled  out.  On  the  other  hand  the  spheroid 
could  result  from  stresses  produced  by  satellite  spin.  Of  non- spherical 
shape  possibilities,  the  prolate  spheroid  is  the  most  nearly  spherical 
consistent  with  measured  brightness  fluctuations.  Whatever  shape  and/ 
or  surface  condition  is  responsible  foi  producing  the  deep  brightness 
minima,  it  persisted  for  at  least  the  twenty  days  between  the  25  August 
and  15  September  (1967)  records,  but  may  have  existed  considerably 
longer . 

Conclusions 

From  the  generally  constant  brightness  on  the  early  orbits  and  the 
good  agreement  of  the  measured  brightness  to  that  expected  for  a  100- 
foot  sphere,  successful  deployment  of  PAGEOS  and  its  full  inflation 
were  achieved.  The  brightness  irregularities  evident  only  six  hours  after 
deployment  indicate  definite  departure  of  the  deployed  vehicle  from  a 
smooth  sphere.  It  should  be  noted  that  the  slow  0.  25  magnitude  fluctuations 
observed  on  Rev.  3  are  consistent  with  the  conclusions  of  the  Smithsonian 
Astrophysical  Observatory1  4  that  the  vehicle  inflation  was  perfect,  based 
on  apparent  constancy  of  the  trailed  image  at  01*128rn26s,  (20  seconds  after 
complete  inflation).  The  photographic  record  is  simply  not  capable  of 
showing  such  small  modulations  in  brightness  of  a  trailed  image.  The  data 
presented  here  suggest  initial  rotation  of  PAGEOS  about  its  center  of  mass, 
and  show  what  appears  to  be  a  rate  of  change  pattern.  The  gradual  evolu¬ 
tion  in  late  June  and  early  July  1966  of  the  relatively  small  initial  bright¬ 
ness  discontinuities  into  the  major  ones  of  mid- July  indicates  that  these 
deformations  of  the  balloon  material  were  due  to  evolution  of  stresses  in 
the  skin  rather  than  due  to  the  sublimed  inflation  material  as  suggested  by 
others.7  However,  the  discontinuities  of  form  were  most  drastic  some 
two  to  ten  days  after  the  satellite  grazed  the  earth's  shadow  for  the  first 
time  and  thus  skin  temperatures  would  have  dropped  to  as  low  as  -100°C 
and  allowed  either  the  benzoic  acid  or  anthraquinone  inflation  gases,  if 
any  remained,  to  solidify  and  accrete  on  the  skin.  Reference  6  contains 
an  estimate  that  the  inflation  gases  should  essentially  be  lost  by  revolu¬ 
tion  80  (4  July  1966),  based  on  a  projected  thermal  history.  However, 
such  thermal  shocking  of  the  plastic  as  would  have  occurred  in  this 
temperature  cycling  would  also  have  intensified  mechanical  stress  redistri¬ 
bution.  This  point  probably  deserves  more  study. 

From  the  relatively  large  amplitude  of  the  atmospheric  scintillation, 
as  well  as  from  the  lack  of  a  phase  angle  effect,  one  can  conclude  that 
the  vehicle  appears  as  a  bright  stellar  point  and  hence  is  a  highly  specular 
convex  form,  like  Echo  I,  and  unlike  Echo  II  or  the  4-meter  Explorer 
balloon  satellites  or  even  the  9-meter  USAF  evaporative  grid-sphere 
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(1966-63A),  all  of  which  exhibit  diffuse  scattering  characteristics  and 
show  lower  levels  of  scintillation  as  would  be  expected  of  objects 
subtending  angular  extents  of  0.  5  to  5  arc  seconds. 

The  one  practical  complication  which  the  brightness  variations  of 
the  PAGEOS  vehicle  might  produce  is  a  difficulty  in  continuous  detection 
of  the  vehicle  near  its  maximum  slant  range  for  ballistic  cameras  of 
limited  aperture,  such  as  the  BC-4  and  MOTS  cameras.  This  is  almost 
certain  to  be  true  since  the  minimum  brightness  events  will  correspond 
to  approximately  stellar  magnitude  +  4. 8  at  the  slant  ranges  of  10,  000 
km  which  are  expected  later  in  its  life  (Ref.  6).  Alternately,  the 
occurrence  of  bright  maxima  might  produce  camera  images  which  are 
larger  than  the  40-60  microns  desired  for  best  determination  of  the 
image  location.  The  computation  of  nominal  exposure  time  and  camera 
aperture  to  be  used  at  a  given  site  on  each  transit  will  prove  to  be  less 
satisfactory  than  would  have  been  possible  without  the  brightness  varia¬ 
tions  . 

It  would  appear  that  for  the  purpose  of  remote  satellites  intended 
to  have  a  brightness  which  is  both  large  and  nearly  invariant  (except  for 
predictable  phase  angle  variations)  a  diffusely-scattering  inflatable 
sphere  would  be  the  best  design.  This  type  of  vehicle  would  also  suffer 
the  least  from  atmospheric  scintillation.  The  image  displacement  from 
the  center  of  mass  due  to  the  phase  effect  might  prove  less  troublesome 
than  the  brightness  variations  of  random  rather  than  predictable  character. 
The  maximum  image  displacement  of  PAGEOS  is  only  some  1.  0  arc 
seconds  at  a  2400  km  slant  range  and  the  image  size  only  some  2  arc 
seconds,  comparable  if  not  smaller  than  the  seeing  disc  for  all  but  the 
best  of  observing  sites.  However,  if  a  collapsible  satellite  is  required 
with  a  specular  surface  to  produce  a  nearly-constant  point-like  source, 
the  use  of  some  of  the  new  polyurethane  "memory-stressed"  foams  might 
be  used  to  establish  sphericity  to  a  predictable  tolerance. 

It  has  been  shown  that  continuous  monitoring  by  broad-band  photo¬ 
electric  photometry  with  due  care  to  obtain  continuity  of  record  and 
accuracy  of  calibration,  is  a  particularly  useful  data  source  in  providing 
information  concerning  PAGEOS'  size,  form,  surface  characteristics, 
and  rotation  about  its  center  of  mass.  Other  workers  have  shown  that 
additional  information  concerning  surface  characteristics  is  available 
by  use  of  narrow-band  and  polarimetric  sensing,  7,13  but  unless  these 
ancillary  data  are  sensed  simultaneously  (‘.hrough  the  use  of  beam 
splitters  or  additional  telescopes),  continuity  of  the  often  complicated 
brightness  fluctuation  patterns  is  lost  by  the  interruptions  in  the  record. 
This  is  even  more  important  if  the  vehicle  is  rapidly  rotating  and/  or 
non -spherically  symmetrical  as  in  the  case  of  most  payloads  and  all 
upper  stage  rockets. 
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